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1 Introduction. 

The realization of Bose Einstein condensates (BEC) and quantum degenerate Fermi gases with 
cold atoms have been highlights of quantum physics during the last decade. Cold atoms in the 
tens of nanokelvin range are routinely obtained via combined laser- and evaporative-cooling tech- 
niques-l^ For high-enough densities (> 10^^ cm^), the atomic de Broglie wavelength becomes 
larger than the typical interparticle distance and thus quantum statistics governs the many-body 
dynamics of these systems. Characteristic features of the physics of cold atomic gases are the mi- 
croscopic knowledge of the many -body Hamiltonians which are realized in the experiments and the 
possibility of controlhng and tuning system parameters via external fields. External field control 
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of contact inter-particle interactions can be achieved, for example, by varying the scattering length 
via Feshbach resonances,^ while trapping of ultracold gases is obtained with magnetic, electric and 
optical fields.'^ In particular, optical lattices, which are artificial crystals made of light obtained via 
the interference of optical laser beams, can realize perfect arrays of hundreds of thousands of mi- 
crotraps,!^ allowing for the confinement of quantum gases to one-dimensional (ID), 2D and 3D 
geometries and even the manipulation of individual particles.'^ This control over interactions and 
confinement is the key for the experimental realization of fundamental quantum phases and phase 
transitions as illustrated by the BEC-BCS crossover in atomic Fermi gases,'^ and the Berezinskii- 
Kosterlitz-Thouless transition"^ for cold bosonic atoms confined to 2D. 

Breakthroughs in the experimental realization of BEC and degenerate Fermi gases of atoms 
with a comparatively large magnetic dipole moment, such as ^'^Cvj^^^ ^^^Ei^^and ^^^Dy atoms^^^^ 
(dipole moment 6i1b^1}1b and 10/ig, respectively, with /i^ Bohr's magneton), and the recent as- 
tounding progress in experiments with ultracold polar molecules^^^^ have now stimulated great 
interest in the properties of low temperature systems with dominant dipolar interactions (see re- 
views Refs.'^^^^ for discussions of various aspects of the problem). The latter have a long-range 
and anisotropic character, and their relative strength compared to, e.g., short-range interactions 
can be often controlled by tuning external fields, or else by adjusting the strength and geometry of 
confining trapping potentials. For example, in experiments with polarized atoms, magnetic dipo- 
lar interactions can be made to overcome short-range interactions by tuning the effective i'-wave 
scattering length to zero using Feshbach resonances.^^^^ This has already led to the observation 
of fundamental phenomena at the mean-field level, such as, the anisotropic deformation during 
expansion and the directional stability'^^'^ of dipolar BECs. Heteronuclear polar molecules in 
a low vibrational and rotational state, on the other hand, can have large permanent dipole mo- 
ments along the intemuclear axis with strength ranging between one tenth and ten Debye (1 De- 
bye c:^ 3.335 X 10^-^^ C -m). In the presence of an external electric field (with a typical value of 
10-^ — 10^ V/cm) mixing rotational excitations, the molecules can be oriented in the laboratory 
frame and the induced dipole moment can approach its asymptotic value, corresponding to the 
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permanent dipole moment. This effect can be used to tune the strength of the dipole-dipole in- 
teraction.!^ Additional microwave fields allow for advanced tailoring of the interactions between 
the molecules, where even the shape of interaction potentials can be tuned with external fields, in 
addition to the strength. This tunability of interactions forms the basis for the realization of novel 
quantum phenomena in these systems, in the strongly interacting limit. 

As a result of this progress, in recent years dipolar gases have become the subject of intensive 
theoretical efforts, and there is now an extensive body of literature predicting novel properties for 
these sy stems .'^^'^ It is the purpose of this review to provide a summary of these recent theoretical 
studies with a focus on the many-body quantum properties, to demonstrate the connections and dif- 
ferences between dipolar gaseous systems and traditional condensed-matter systems, and to stress 
the inherent interdisciplinary nature of these studies. This work covers spatially homogeneous as 
well as trapped systems, and includes the analysis of the properties of dipolar gases in both the 
mean-field (dipolar Bose-Einstein condensates and superfluid BCS pairing transition) and in the 
strongly correlated (dipolar gases in optical lattices and low-dimensional geometries) regimes. 

We tried our best to include all relevant works of this exciting, ever expanding field. We 
apologize in advance if some papers (hopefully, not many) are not appearing below. 

2 The dipole-dipole interaction 

For polarized dipolar particles, interparticle interactions include both a short-range Van der Waals 
(vdW) part and a long-range dipole-dipole one. The latter is dominant at large interparticle sep- 
arations and assuming a polarization along the z-axis as in Fig. [TJa) the interparticle interaction 
reads 

Vdrf(r) = 3(l-3cos2 0). (1) 

Here d is the electric dipole moment (for magnetic dipoles d^ should be replaced with jU^, with /i 
the magnetic dipole moment), r is the vector connecting two dipolar particles, and 9 is the angle 
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Figure 1: (a) Geometry for the interaction of two aligned dipoles. (b) Tunability of the dipole- 
dipole interaction by using a time-varying aligning field. The angle between the dipole orienta- 
tion and the z-axis determines the strength and the sign of the effective interaction. 

between r and the dipole orientation (the z-axis). The potential ^^^^(r) is both long-range and 
anisotropic, that is, partially repulsive and partially attractive. As discussed in the sections below, 
these features have important consequences for the scattering properties in the ultracold gas, for 
the stability of the system as well as for a variety of its properties. 

2.1 Scattering of two dipoles 

The long-range character (~ r^^) of the dipole-dipole interaction results in all partial waves con- 
tributing to the scattering at low energies, and not only, e.g., the ^•-wave, as is often the case for 
short-range interactions. In fact, for dipole-dipole interactions the phase shift 5i in a scattering 
channel with angular momentum I behaves as 5/ ~ for Z > and small k (see, e.g., Refs.'^ 



The effect of the anisotropy of the interaction is instead that the angular momentum is not 
conserved during scattering: for bosons and fermions the dipole-dipole interaction mixes all even 
and odd angular momenta scattering channels, respectively. Due to the coupling between the 
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various scattering channels, the potential V^d then generates a short-range contribution to the total 
effective potential in the 5- wave channel (/ = 0). This has the general effect to reduce the strength 
of the short-range part of the interaction. 

Thus, for two bosonic dipolar particles (even angular momenta) the scattering at low energies 
is determined by both the long-range and the short-range parts of the interaction. This is in contrast 
to the low energy scattering of two fermionic dipoles (odd angular momenta), which is universal 
in the sense that it is determined only by the long-range dipolar part of the interaction, and is 
insensitive to the short-range details. 

For a dilute weakly interacting gas the above results allow a parametrization of the realistic 
interparticle interaction between two particles of mass m in terms of the following pseudopoten- 
tiall3oLllI](see also Refs.E^l and^Sl) 

y(r) =^5(r) + ^(l-3cos2 0), (2) 

with 

Anfp-aid) 

g = — (3) 

m 

parametrizing the short-range part of the interaction. We note that the long-range part of the pseu- 
dopotential y (r) is identical to the long-range part of the original potential and the scattering length 
a{d) controlling the short-range part depends on the dipole moment. This dependence is impor- 
tant'^ when one changes the dipole moment, using external, e.g. electric, fields, as explained 
below. 



The strength of the dipole-dipole interaction can be characterized by the quantity 

ad = 

fi 



7 



which has the dimension of length and can be considered as a characteristic range of the dipole- 
dipole interaction, or dipolar length. This length determines the low energy limit of the scattering 
amplitudes, and, in this sense, Ud is analogous to the scattering length for the dipole-dipole interac- 
tion. For chromium atoms with a comparatively large magnetic moment of 6/ifi (equivalent dipole 
moment d — 0.056 Debye) we have Ud ~ 2.4 nm. For most polar molecules the electric dipole 
moment ranges in between 0.1 and 1 Debye, while Ud ranges from 1 to 10^ nm. For example, the 
dipole moment of fermionic ammonia molecules ^^NDs is J = 1.5 Debye with = 712 nm, while 
for H^^C^^N it increases to d = 2.98 Debye and Ud = 3620 nm. This latter value of the effective 
scattering length is an order of magnitude larger than, for example, the one for the intercomponent 
interaction in the widely discussed case of a two-species fermionic gas of ^Li, where au = — 1 14 
nm. Thus, the strength of the dipole-dipole interaction between polar molecules can be not only 
comparable with but even much larger than the strength of the short-range interatomic interaction. 



2.2 l\inability of the dipole-dipole interaction 



One spectacular feature of the dipole-dipole interaction is its tunability. In Sect. 2.2.1 we first 



review methods for tuning the strength and sign of dipolar interactions with an eye to cold atoms. 



and then in Sect. 2.2.2 we discuss tunability for the specific case of polar molecules, where both 



the strength as well as the shape of interactions can be engineered. 
2.2.1 TunabUity of interactions in cold atoms 

In Ref.!^ a technique has been developed to tune the strength as well as the sign of dipolar in- 
teractions in atomic systems with a finite permanent magnetic dipole moment. This technique 
uses a combination of a static (e.g., magnetic) field along the z-axis and a fast rotating field in the 
perpendicular xy-p\a.m such that the resulting time dependent dipole moment is [see Fig. ^h)] 

d{t) = d {cj-cos^ + [exCos{Q.t) +eySm{Q.t)] sin^j. 
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Here Q. is the rotating frequency of the field and the angle ^,O<0<;r/2is determined by the 
ratio of the amplitudes of the static and rotating fields. The above expression implies that the 
dipoles follow the time-dependent external field adiabatically. This in turn sets an upper limit on 
the values of the rotating frequency Q., which should be (much) smaller than the level splitting in 
the field. However, if the frequency Q. is much larger that the typical frequencies of the particle 
motion, over the period In/Q. the particles feel an average interaction 

(Vrf(r)) = ^(l-3cos20)a(,^). 

The latter differs from the interaction for aligned dipoles, Eq. ([T]), by a factor a{^) = (3cos^ ^ — 
l)/2, which can be changed continuously from 1 to —1/2 by varying the angle 0. Thus this 
method allows to "reverse" the sign of the dipole-dipole interaction and even cancel it completely 
for ^ = a cos 1 / -\/3 = 54.7, similar to NMR techniques. We note that an analogous technique 
can be also applied for the electric dipole moments of, e.g., polar molecules. We will review 
applications of this method below. 

2.2.2 Effective Hamiltonians for polar molecules 

In the following we will be often interested in manipulating interactions for polar molecules in 
the strongly interacting regime. In particular, we will aim at modifying not only the strength but 
also the shape of interaction potentials, as a basis to investigate new condensed matter phenom- 
ena. This usually entails a combination of the following two steps: (i) manipulating the internal 
(electronic, vibrational, rotational, ...) structure of the molecules, and thus their mutual interac- 
tions, using external static (DC) electric and microwave (AC) fields, and (ii) confining molecules 
to a lower-dimensional geometry, using, e.g., optical potentials, as exemplified in Fig. [2j Under 
appropriate conditions, the resulting effective interactions can be made purely repulsive at large 
distances (e.g., at characteristic distances of lOnm or more), as in the 2D example of Fig. [3]^a). 
On one hand, this has the effect to suppress possible inelastic collisions and chemical reactions 
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occurring at short-range (i.e., at characteristic distances of Uc < Inm), and on the other hand it 
allows to study interesting condensed matter phenomena originating from the often-unusual form 
of the two-body (or many -body-) interaction potentials. In the next few subsections we review 
techniques for the engineering of the interaction potentials which will be used in the many-body 
context in Sect. [6]below. 

Our starting point is the Hamiltonian for a gas of cold heteronuclear molecules prepared in 
their electronic and vibrational ground-state, 



Here the first term in the single particle Hamiltonian corresponds to the kinetic energy of the 
molecules, while Vtrap(ri) represents a trapping potential, as provided, for example, by an optical 
lattice, or an electric or magnetic trap. The tcrmH^ describes the internal low energy excitations 
of the molecule, which for a molecule with a closed electronic shell ^E(v = 0) (e.g. SrO, RbCs 
or LiCs) correspond to the rotational degree of freedom of the molecular axis. This term is well 
described by a rigid rotor H^^ = = 5J? with B the rotational constant (in the few to tens of 
GHz regime) and J/ the dimensionless angular momentum. The rotational eigenstates for a 

quantization axis z, and with eigenenergies 57(7-1-1) can be coupled by a static (DC) or microwave 
(AC) field E via the electric dipole moment d/, which is typically of the order of a few Debye. 

In the absence of electric fields, the molecules prepared in a ground rotational state 7 = have 
no net dipole moment, and interact via a van-der-Waals attraction Vydw ~ —C^jr^, reminiscent of 
the interactions of cold alkali metal atoms in the electronic ground- states. Electric fields admix 
excited rotational states and induce static or oscillating dipoles, which interact via strong dipole- 
dipole interactions Y^d with the characteristic Xjr' dependence given in Eq. ([T]). For example, a 
static DC field couples the spherically symmetric rotational ground state of the molecule to excited 
rotational states with different parity, thus creating a non-zero average dipole moment. The field 




+ Vtrap(r,) - d,-E(0 -f £ ydd(r, - r,). 



(4) 
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strength therefore determines the degree of polarization and the magnitude of the dipole moment. 
As a result, the effective dipole-dipole interaction may be tuned by the competition between an ori- 
enting, e.g., DC electric field and the quantum (or thermal) rotation of the molecule. This method 
effectively works for the values of the field up to the saturation limit, at which the molecule is 
completely polarized (typically 10^ V/cm). 



The many body dynamics of cold polar molecules is thus governed by an interplay between 
dressing and manipulating the rotational states with DC and AC fields, and strong dipole-dipole 
interactions. In condensed matter physics one is often interested in effective theories for the 
low-energy dynamics of the many-body system, after the high-energy degrees of freedom have 
been traced out. The connection between the full molecular A^-particle Hamiltonian (|4]) includ- 
ing rotational excitations and dressing fields, and an effective low-energy theory can be made 
using the following Born-Oppenheimer approximation: The diagonalization of the Hamiltonian 



Hbo = E; 



+ El<jKid(r;- for frozen spatial positions {r,} of the molecules 
yields a set of energy eigenvalues Vjp ({r/}), which can be interpreted as the effective interaction 
potential in the single-channel many-body Hamiltonian 



A' 

3D 



i=l I 



g; + V„,(r,) 



+ viY(M)- (5) 



The term V^f {{fi}) represents an effective A'^-body interaction, which can be expanded as a sum 
of two-body and many-body interactions 

({r,}) = t^'"' (r,-r,) + £ (r,-,r„r,) + . . . , (6) 

i< j i< j<k 

where in most cases only two-body interactions are considered. The dependence of ^^.^({r/}) 
on the electric fields E provides the basis for the engineering of the many body interactions, as 
described below. 
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We note that the attractive part of the interaction potential can induce instabilities in a dipolar 
gas at the few body level as well as at the many-body level (this latter case will be discussed 
in Sect. |3] below). For example, for several experimentally relevant mixed aUcali-metal diatomic 
species such as KRb, LiNa, LiK, LiRb, and there exist chemically reactive channels that 

are energetically favorable, leading to particle recombination and two-body losses in the gas. The 
rate of chemical reactions can be strongly enhanced by dipole-dipole interactions which can attract 
molecules in a head-to-tail configuration [e.g., = in Fig. (TJ^)] to distances on the order of 
typical chemical interaction distances, ac ^ 1 nm.l^^^ One aim of interaction engineering is to 
control these interactions in order to stabilize the gas against particle losses. This will enable the 
study of complex condensed matter phenomena in these systems. 

2.2.3 Stabilization of dipolar interactions in 2D 

The simplest example of stabilization of dipolar interactions against inelastic collisions is sketched 
in Fig.[2]and consists of a system of cold polar molecules in the presence of a polarizing DC electric 
field oriented in the z-direction, and of a strong harmonic transverse confinement Vtrap = 



with frequency CI and characteristic length ah = ^Jh/{m£L). The latter is provided, e.g., by an 
optical lattice along z. 

Figure [3|b) shows a countour plot of the interaction potential V{p,z) for two dipoles in this 
quasi-2D geometry, where 

Here p = {p,^,z) represents the distance between the two molecules in cylindrical coordinates, 
and r = |p|. The first term is the isotropic vdW potential, assuming the molecules are in their 
rotational ground state, with a vdW length a = (2;r/r(^)^) {ijiCf^/h^^ ^/^j58159| second term 
is the anisotropic dipolar potential, with induced dipole moment d and dipolar length a^- 

Figure |3];b) illustrates essential features of reduced dimensional collisions: for finite d, the 
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Figure 2: System setup: Polar molecules are trapped in the (x, j)-plane by an optical lattice made of 
two counter-propagating laser beams with wavevectors ±kL = i^LCz. (blue arrows). The dipoles 
are aligned in the z-direction by a DC electric field Edc = ^^DC^z (red arrow). An AC microwave 
field is indicated (green arrow). Inset: Definition of polar (i?) and azimuthal (cp) angles for the 
relative orientation of the inter-molecular collision axis ri2 with respect to a space-fixed frame, 
with axis along z. [Adapted from Ref.!^ 

repulsive dipole-dipole interaction overcomes the attractive van-der-Waals potential in the (z = 0)- 
plane at distances r > cic, realizing a repulsive in-plane potential barrier (blue dark region). In 
addition, the harmonic potential confines the particles's motion in the z direction. The combination 
of the dipole-dipole interaction and of the harmonic confinement thus yields a three-dimensional 
potential barrier separating the long-distance, where interactions are repulsive, from the short- 
distance one, where interactions are attractive and inelastic processes can occur. If the collision 
energy is smaller than the height of the barrier at the saddle point (white circles), the particles' 
motion is confined to the long-distance region, where particles scatter elastically. Particle losses 

(re) 

are due to tunneling through the potential barrier at a rate Kj . Within a semiclassical (instanton) 

(re) 

approximation valid for ^ max{a/,,a}, the tunneling rate Kj ' is well approximated by the 

13 



(a) 



(b) 



MI^[9) [arb. units] 





p [arb. units ] 



Figure 3: (a) Effective repulsive interaction potentials in 2D. Solid line: Dipolar potential 
yJP(p) = D/p^ induced by a DC electric field in the configuration of Fig. 2] [Edc > and 
Eac(0 = ]. Dash-dotted line: "Step-like" potential induced by a single AC microwave field 
and a weak DC field [see also Fig. |5|. Dashed line: Attractive potential induced by the combi- 
nation of several AC fields and a weak DC field. The potentials V'JP(p) and the separation p are 
given in arbitrary units, (b) Potential V{p,z)/ {hQ) versus p/tth and z/a/, for |m| = 1 for ^OR^^j^b 
in a t2 = 2;r(50kHz) trap, where ah = 56.4 nm, d = 6.25 nm'^and d = 0.2 D. White circles: sad- 
dle points. Dashed lines: semiclassical trajectory for the collision of two molecules. The dashed 
half-circles show d and a^. [Panels (a) and (b) adapted from Refs-f^and,!^ respectively.] 



exponential form 



(7) 



The constant c has been recently computed numerically by Julienne, Hanna and Idziaszekl^ to 
be c ~ 3.03, while the "attempt rate" ft)p^^ for the scattering of two isolated dipoles reads cOp oc 
{hK^a^/ ju), independent of particles' statistics. Here K is the momentum for a collision with rela- 
tive kinetic energy = tP'K^ / (2/i), with Uk = 2k /k the DeBroglie wavelength. For particles in a 



crystalline configuration (see Sect. 6.1 below), (Op will be proportional to the frequency of phonon 



oscillations around the mean particle positions (Dp ~ d^ / (jla^) , with a the mean inter-particle 
distance. The expression Eq. (|7]) shows that coUisional losses may be strongly suppressed for any 
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molecular species for a large enough dipole moment or strength of transverse confinement. 



In ultracold collisions one often has the following separation of length scales: ^ <3/j ^ 
d ^ Gc, and can be tuned by, e.g, increasing the external DC field. Figure Qa) and (b) show 
numerical results for reactive and elastic collision rates of bosonic and fermionic KRb molecules, 
respectively, and for several strengths of transverse confinement. Here a^, ah and d are on the 
order of hundreds of nm, tens of nm, and less than 10 nm, respectively. Because of the moderate 
dm&x = 0.5D of KRb molecules, here Ud < an and the semiclassical regime of large aa of Eq. Q 
is not reached. Nevertheless, in stark contrast to collisions in 3D,'^ the figure shows that the 
ratio between elastic and inelastic collision rates increases rapidly with d, signaling an increased 
stability with increasing d. For bosons, the exact numerical results (thick lines) approach rapidly 
the semiclassical instanton limit (thin lines) with increasing d. The behavior of the inelastic rates 
for fermions is explained in detail in Refs.'^^'^ 

Recent landmark experimental results from the JILA group with fermionic KRb molecules 
show a strong suppression of inelastic collisions and increase of elastic ones with d, in excellent 
agreement with the predictions of Fig. |3] This opens the way to the study of strongly correlated 
phenomena in these systems. 

2.2.4 Advanced interaction designing: Blue-shielding 

By combining DC and AC fields to dress the manifold of rotational energy levels it is possible 
to design effective interaction potentials (r,— r^) with (essentially) any shape as a function 
of distance. For example, the addition of a single linearly-polarized AC field to the configuration 
of Sect. |2.2.3| leads to the realization of the 2D "step-like" potential of Fig. [3]^a) (black dashed- 



dotted line), where the character of the repulsive potential varies considerably in a small region of 
space. The derivation of this effective 2D interaction is sketched in Fig. [sj^a-c).'^^^ The (weak) 
DC-field splits the first-excited rotational (7 = l)-manifold of each molecule by an amount h5, 
while a linearly polarized AC-field with Rabi frequency Q. is blue-detuned from the {\g) — \e))- 
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(b) fermions (j=1) 
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Figure 4: Quasi-2D elastic (el) and reaction (re) rate constants Kj for (a) identical KRb bosons 
(j = 0) and (b) identical KRb fermions (j = 1) at a collision energy of £ = fcs(240nK) = /^(SkHz) 
for three different trap frequencies Q./{27r) = 25 kHz, 50 kHz, and 150 kHz. Heavy lines: CC. 
Light lines: UBA or instanton. Horizontal lines show the unitarity limits. Vertical arrows show 
where rjj = 100. [From Ref.!^ 



transition by HA, see Fig. [5[a). Because of h5 and the choice of polarization, for distances p ^ 
{d^/hSy^^ the relevant single-particle states for the two-body interaction reduce to the states |^) 
and \e) of each molecule. Figure [Sj^b) shows that the dipole-dipole interaction splits the excited 
state manifold of the two-body rotational spectrum, making the detuning A position-dependent. 
As a consequence, the combined energies of the bare groundstate of the two-particle spectrum 
and of a microwave photon become degenerate with the energy of a (symmetric) excited state at 
a characteristic resonant (Condon) point pc = (J^/M)^/-^, which is represented by an arrow in 
Fig.[5]^b). At this Condon point, an avoided crossing occurs in afield-dressed picture, and the new 
(dressed) groundstate potential inherits the character of the bare ground and excited potentials for 
distances p ^ pc and p <ti pc, respectively. Fig. [5|d) shows that the dressed groundstate potential 
(which has the largest energy) is almost flat for p ^ pc and it is strongly repulsive as l/p^ for 
p ^ pc, which corresponds to the realization of the step-like potential of Fig. ^a). We remark 
that, due to the choice of polarization, this strong repulsion is present only in the plane z = 0, 
while for z 7^ the groundstate potential can become attractive. The optical confinement along z 
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Figure 5: Design of the step-like potential of Fig. |3[a) (black dashed-dotted line): (a) Rotational 
spectrum of a molecule in a weak DC field. The DC field splits the (7 = 1) -manifold by an 
amount h5. The linearly -polarized microwave transition with detuning A and Rabi frequency Q. 
is shown as an arrow, (b) BO-potentials for the internal states for Q = (bare potentials), where 
\g;e)± = {\g;e) ± \e;g))/\/2. The resonant Condon point pc is indicated by an arrow, (c) AC- 
field-dressed BO-potentials. The dressed groundstate potential has the largest energy. 



of Sect. 2.2.3 is therefore necessary to ensure the stability of the system. 

The interactions in the presence of a single AC field are described in detail in Ref.,!^ where 
it is shown that in the absence of external confinement this case is analogous to the (3D) optical 
blue- shielding developed in the context of ultracold collisions of neutral atoms ,1^^^^ however with 
the advantage of the long lifetime of the excited rotational states of the molecules, as opposed to 
the electronic states of cold atoms. The strong inelastic losses observed in 3D collisions with cold 
atomsl^^^^ can be avoided via a judicious choice of the field's polarization, eventually combined 
with a tight confinement to ensure a 2D geometry (as e.g. in Fig. [2] above). For example, in Ref.^^ it 
is shown that in the presence of a DC field and of a circularly polarized AC field the attractive time- 
averaged interaction due to the rotating (AC-induced) dipole moments of the molecules allows for 
the cancelation of the total dipole-dipole interaction. The residual interactions remaining after 
this cancelation are purely repulsive 3D interactions with a characteristic van-der-Waals behavior 
(r) ~ {d^/ hA) I r^ . This 3D repulsion provides for a shielding of the inner part of the interaction 
potential and thus it will strongly suppress inelastic collisions in experiments. 

Recent works'®!^ have considered the microwave spectra of alkali-metal dimers including hy- 
perfine interactions. It is an important open question to determine the effects that the presence 
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of internal states such as, e.g. hyperfine states, have on the broad class of shielding techniques 
described above. 



3 Weakly interacting dipolar Bose gas 
3.1 BEC in a spatially homogeneous gas. 

Let us discuss now the influence of the dipole-dipole interaction on the properties of a homoge- 
neous single-component dipolar Bose gas This can be most conveniently done in the language 
of second quantization. For this purpose we introduce particle creation and annihilation field op- 
erators \j/^{r) and i/A(r) satisfying standard bo sonic commutation relation 



i^(r),v/^(r') = VA(r)v/t(rO- V/^(r')VA(r) = 5(r-r') 



0. 



The corresponding second quantized Hamiltonian of the system then reads 



H 



-—V^-ll 
2m 



1 



+ ^ / JrJr>t(r)i^At(r')y(r-r')i^A(r')i^A(r), 



(8) 



where m is the mass of the particles, V{r) is the interparticle interaction, and the chemical potential 
jU fixes the average density n of the gas. We consider the case when the system is away from any 
"shape" resonances,'^,'22l2il ^nd, therefore, replace the original interparticle interaction with the 



'This and the next Sections are substantially revised and updated version of the corresponding part of Ref.^^^l 
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pseudopotential O). Assuming that the system is dilute, na ^ 1, we can write the Hamiltonian as 



H = / dY^f\r) 



+ 



2m 2 



^ I drdr'xif^ (r) (r')y^ (r - r') xif{r') v/(r) , 



(9) 



where is given by Eq. ([Tj) and g = Anh^a/m [as compared to Eq. we omit the dependence of 
the scattering length a on the dipole moment d]. Note that the scattering length has to be positive, 
a > 0, to avoid an absolute instability due to local collapses.!^ 

As we will see below, an important parameter that determines the properties of the system 
described by the Hamiltonian ([9]) is 



£dd 



4nd^ 



md 



1 ad 



3 g 3h^a 3 a 



(10) 



It measures the strength of the dipole-dipole interaction relative to the short-range repulsion. In the 
case £dd < 1 , the short-range part of the interparticle interaction is dominant while the dipole-dipole 
interaction results in only small corrections. For a positive scattering length a, the system is stable 
and exhibits BEC at low temperatures. This case corresponds to earlier experiments'^^^ with Cr 
BEC {£dd ~ 0.16^. It was found that the corrections due to magnetic dipole-dipole interaction 
between ^^Cr atoms are of the order of 10%. 

For the opposite case £dd > 1, the anisotropic dipole-dipole interaction plays the dominant 
role resulting in instability of a spatially homogeneous system,'^,'^.'^This instability can be seen 
in the dispersion relation £'(p) between the energy E and the momentum p of excitations in the 
Bose-condensed gas, which can be easily obtained within the standard Bogoliubov approach: 



^(P) - 




(11) 



2m 



p^+4mng[l + £dd{icos^e-l)]. 
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Here 6 is the angle between the excitation momentum p and the direction of dipoles, and V(p) = 
g+ (4;r<i^/3)(3cos^0 — 1) is the Fourier transform of V{r) = g5{r) +Vd{r). For e^^ > 1, the 
excitation energies ^(p) at small p and 6 close to ;r/2 become imaginary signalling the instability 
(collapse). This instability of a spatially homogeneous dipolar Bose gas with dominant dipole- 
dipole interaction is a result of a partially attractive nature of the dipole-dipole interaction. 

3.2 BEC in a trapped gas. 

The above consideration shows that the behavior of a spatially homogeneous Bose gas with a strong 
dipole-dipole interaction is similar to that of a Bose gas with an attractive short-range interaction 
characterized by a negative scattering length a < 0. In the latter case, however, the collapse of 
the gas can be prevented by confining the gas in a trap provided the number of particles A^^ in the 
gas is smaller than some critical value Nc, N < Nc (see, e.g.,^^). This is due to the finite energy 
difference between the ground and the first excited states in a confined gas. For a small number 
of particle this creates an effective energy barrier preventing the collapse and, therefore, results 
in a metastable condensate. The same arguments are also applicable to a dipolar BEC in a trap, 
see Refs.'^and,'^with one very important difference: The sign and the value of the dipole-dipole 
interaction energy in a trapped dipolar BEC depends on by the trapping geometry and, therefore, 
the stability diagram contains the trap anisotropy as a crucial parameter. 

3.2.1 Ground state. 

The Hamiltonian for a trapped dipolar Bose gas reads 




,2 



H 




(12) 



+ 
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where 



(13) 



is the trapping potential and we again use the pseudopotential (|2]) for the short-range part of the in- 
terparticle interaction assuming that the system is away from "shape" resonances. For the trapping 
potential we consider the experimentally most common case of an axially symmetric harmonic 
trap characterized by the axial (O^ and radial (Op trap frequencies. The aspect ratio of the trap / 
is defined through the ratio of the frequencies: / = y^OJp/ojv = Iz/lp, where 1^ = ^/hjmco^ and 
Ip = ^Jh/m(Op are the axial and radial sizes of the ground state wave function in the harmonic 
oscillator potential ( [T3| ), respectively. For Z < 1 one has a pancake-form (oblate) trap, while the op- 
posite case / > 1 corresponds to a cigar-form (prolate) trap. Taking into account the anisotropy of 
the dipole-dipole interaction, one can easily see that the aspect ratio / should play a very important 
role in the behavior of the system. 

The standard mean-field approximation corresponds to taking the many-body wave function in 
the form of a product of single-particle wave functions: 

^(ri,...,r^;0 = nV^i(r;'0- (14) 

The condensate is then described by the condensate wave function ^f{r,t) = \/N\if\{Y,t) normal- 
ized to the total number of particles, / dr \ i//(r, t)\ = N, and governed by the time-dependent 
Gross-Pitaevskii (GP) equation 



-^V^-li + U,rir)+gMrj)\^ + 



,2 f , ,l-3cos20 I , 

+ / dr ^ v/(r ,f 

J Ir-r'l 



V/(r,f). (15) 



The validity of this approach was tested in Refs. EUandEH by using many-body diffusion Monte- 
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Carlo calculations with the conclusion that a GP equation with the pseudopotential (|2]) provides a 
correct description of the gas in the dilute limit na^ <^ 1 . Note that, being the product of single- 



particle wave functions, the many -body wave function ( 14 1 does not take into account interparticle 
correlations at short distances due to their interaction, which takes place at interparticle distances 



|r — r'l < a^/ = mcfi/h^. This change of the wave function is taken into account in Eq. (15) by 



the contact part of the pseudopotential ([2]) [the fourth term in the right-hand-side in Eq. ( [15] )] but 
ignored in the last term of Eq. ( [T5] ) because the main contribution to the integral comes from large 
interparticle distances (of order the spatial size of the condensate). 



Let us first consider stationary solutions of Eq. ( 15 ), for which V^(r,f) = i//b(r) and V^o(r) obeys 
the stationary GP equation 



2 / , / I - 3cos^0 I , .2 
|r — r'l 



V/o(r) = /iV^o(r), 



(16) 



where = Numerical analysis of Eqs. ( 15 ) and ( 16 ) was performed in Ref s . ^BQ1^II11J 74 | 76 L[77] 



on the basis of numerical solutions of the non-linear Schrodinger equation ( 16 ) together with varia- 
tional considerations with the Gaussian ansatz for the condensate wave function. In Refs. 
the problem was treated using diffusive Monte-Carlo calculations, while the authors of Ref.!^ 
apply the Thomas-Fermi approximation that neglects the kinetic energy and allows to obtain ana- 
lytical results. 

We begin the discussion of the results with the case of a dominant dipole-dipole interaction. 



£dd ^ 1» such that the third term in the left-hand- side of Eq. (16) can be neglected. This case 
demonstrates already all important features of the behavior of dipolar condensates. The general 
case will be briefly discussed at the end of this section. 

Let us introduce the mean-field dipole-dipole interaction energy per particle 



17 1 /" J J /I / m2 1 -ScOS^e I ; |2 



(17) 



r — r 
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which together with the trap frequencies co^ and cOp are important energy scales of the problem. 
One can easily see that the value of the chemical potential ji and the behavior of the dipolar 
condensate are determined by the aspect ratio of the trap /, the quantity V/hcOp, and the parameter 
o = Na^/lp. Notice also that the anisotropy of the dipole-dipole interaction results in squeezing 
the cloud in the radial direction and stretches it in the axial one (along the direction of dipoles) 
in order to low the interaction energy. For this reason the aspect ratio of the cloud L = L^/Lpis 
always larger than the aspect ratio / of the trap. Here and Lp are the axial and the radial sizes of 
the cloud, respectively. 

We now summarize the results of the stability analysis of the dipolar condensate with £jj ^ 1 
(Eq. ( 16 1 with g = 0),l^,'^^E9l(see also Ref.l^for the stability analysis in a general harmonic trap). 
The mean-field dipole-dipole interaction is always attractive, V < 0, for a cigar shaped trap / > 1 
causing instability (collapse) of the gas if the particle number exceeds a critical value Nc- This 
critical value depends only on the trap aspect ratio /. It was found that the shape of the cloud with 
N close to Nc is approximately Gaussian with the aspect ratio L 2.1 for a spherical trap (/ = 1), 
and L ^ 3.0 for an elongated trap with / ^ 1. 

For a pancake shaped trap with / < 1, the situation is more subtle. In this case there exists a 
critical trap aspect ratio ^ 0.43, which splits the pancake shaped traps into soft pancake traps 
(/^ < Z < 1) and hard pancake traps (Z < Z^). For soft pancake traps one has again a critical number 
of particles Nc such that the condensates with > Nc are unstable. For close to Nc and Z — t- l^, the 
aspect ratio of the cloud Lc approaches the aspect ratio of the trap, Lc — ?■ Z*. Note that in this case 
the collapse occurs even in a pancake shaped cloud with positive mean dipole-dipole interaction V 



due to the behavior of the lowest quadrupole and monopole excitations (see Section 3.2.2 1. 

For hard pancake traps, it was argued in Refs .l^andl^that the dipolar condensate is stable for 
any because the dipole-dipole interaction energy V is always positive. On the other hand, by 
using more advanced numerical analysis and larger set of possible trial condensate wave functions, 
the authors of Ref.!^ found that the dipolar condensate in a hard pancake trap is also unstable for 
sufficiently large number of particles. Similar conclusions were drawn in Ref.^^ It was found that 
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the critical values of the parameter a for the instability to occur are orders of magnitude larger 
than in soft pancake and cigar shaped traps. In addition, the regions in parameter space were 
discovered where the maximum density of the condensate is not in the center of the cloud such 
that the condensate has a biconcave shape. (Analogous behavior of the condensate in a general 
three-dimensional harmonic trap were found in Ref.jl-^see also^ and.l^ These regions exist also 
in the presence of a small contact interaction with \a\ < O.laj, but their exact position and size 
depend on a. It is important to mention that condensates with normal and biconcave shapes behave 
differently when the instability boundary is crossed. The condensate with a normal shape develops 
a modulation of the condensate density in the radial direction, so-called "radial roton" instability 



similar to the roton instability for the infinite-pancake trap (/ — oa)i2M ggg Section 3.2.3 On the 
other hand, it is the density modulations in the angular coordinate that lead to the collapse of 
biconcave condensates - a kind of "angular roton" instability in the trap. In the latter case one has 
spontaneously broken cylindrical symmetry. 

The behavior of the trapped dipolar condensate can be simply captured by means of a Gaussian 
variational ansatz for the condensate wave function V^(r): 



where the equilibrium radial size Lp and the cloud aspect ratio L can be found by minimizing the 
energy. Note that in order to describe biconcave shaped condensates, one has to consider (see 



Refy^ SL linear combination of two wave functions: the first one is a Gaussian (18 1 and the second 
one is the same Gaussian multiplied by H2{x/Lp)+H2{y/Lp), where H2 is the Hermite polynomial 
of the second order. 

For large values of the parameters Na / li, where / = p or z, are large (but still na^ ^ 1 ), one can 
use the Thomas-Fermi approximation to find the chemical potential and the shape of the cloud.l^ 
This case corresponds to the small the kinetic energy, as compared to other energies, and, therefore. 
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we can neglect the corresponding term with derivatives in Eq. ([16]). The GP equation then becomes 



'^{co^pP^ + (oh^) + ^g\xir{r,t)\^ + 



,2 f , ;i - Scos^e I ,2 



r — r 



V/o(r) = AiV^o(r)• 



(19) 



The solution of this equation reads 



VA^(r) =n(r) =no 1 



4 Rl 



with the chemical potential 



where hq is the density of the condensate in the center of the trap and 



F{L) 



1 GiL) - 1 



3 L2-1 ■ 



(20) 



The energy of the condensate is 



1 . . f \5 



(21) 



and the radii of the condensate in the radial and axial directions Rp and are 



Rr 



Anmay^L 



9 F{L) 
2L2-1 



1/5 



R, = LR 



P' 



(22) 
(23) 



and the corresponding aspect ratio of the cloud L can be found from the equation 



' 2/4 ^ i L2 - 1 



+ (e.w-l)( 1-^ 1 =0. 



(24) 
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Note, that the above equation coincides with the equation on the aspect ratio for the Gaussian 
variational ansatz ( 18 ) when the kinetic energy contribution is neglected, as shown in Ref.l^It was 



also found that the Thomas-Fermi approximation agrees well with numerical results when used to 
analyze the stability of the condensate. However, the critical number of particles cannot be found 



in the Thomas-Fermi approximation because both terms in the expression (21 ) for the energy have 



the same dependence N'^/^ on the number of particles A'^ after taking into account the expressions 



( |22) ) and ([23]) for Rp and R^. 

Let us now briefly discuss the stability of a dipolar condensate in the general case with g ^0. 
It is obvious that for an attractive short-range interaction with ^ < the condensate can only be 
(meta) stable for a small number of particles. For a repulsive short-range interaction with ^ > and 
weak dipole-dipole interaction < e^j < 1, the condensate is always stable. For e^^ > 1 the dipolar 
condensate can be only metastable for number of particles smaller than a critical value, A'^ < A^^^ 
which depends on e^^ and the trap aspect ratio /. This means that the (metastable)condensate 
solution provides only a local minimum of the energy, while the global minimum presumably 
corresponds to a collapsed state with L — )■ oo or, for Z < 1, a kind of density modulated state. 

3.2.2 Collective excitations and instability. 

We have already mentioned that collective excitations play an important role in the stability anal- 
ysis of a dipolar condensate. They also determine the dynamics of the gas and, therefore, are of 
experimental interest. 

For a trapped dipolar condensate, the analysis of excitations is usually performed on the basis 
of the Bogoliubov-de Gennes equations which can be obtained by linearizing the time-dependent 



GP equation ( 15 1 around the stationary solution V^o(r). This can be achieved by writing a solution 



of Eq. ( 15 1 in the form 

V/(r,?) = \ifo{r)+£[u{r)exp{-im)+v*{r)exp{-icot)] 
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where the second term describes small (e ^ 1) oscillations of the condensate around V^(r) with 



(complex) amplitudes u{r) and v(r). To the first order in e the linearization of Eq. ( 15 1 gives the 
Bogoliubov-de Gennes equations 



h(Ou{r) = [-—V^-ii + Utr{r) + 
Zm 

+ ij JrV(r-r')|v/o(r')|^]"(r) + 
+ |jrV(r-r')|vAo(rO|'v(r), 



(25) 



-hQ)v{r) = [-— V2 - ^ + UM + 
Im 

+2 J JrV(r-r')|v^o(r')|^]v(r) + 
+ jdr'V{r-r')\Yo{r')\\{r), 



(26) 



where V{r — r') is given by Eq. The solution of these linear equations provides the eigenfunc- 
tions {uj{r),Vj{r)) with the amplitudes Uj{r) and vj(r) obeying the normalization condition 



dr[u*{r)uj{r)-v*{r)vj{r)] = d, 



and the corresponding eigenfrequencies COj of the collective modes. The Bogoliubov-de Gennes 



equations (25) and (26) can also be obtained by diagonalizing the Hamiltonian (13) in the Bo- 
goliubov approximation, which corresponds to splitting the field operator i^(r) into its mean-field 
value i//b(r) and the fluctuating quantum part expressed in terms of annihilation and creation oper- 
ators dCj and aj of bosonic quasiparticles (quanta of excitations): 



VA(r) = i^o(r)+£K(r)d; + v*(r)a;] 



The normalization condition for the amplitudes Uj{r) and Vj{r) ensures the bosonic nature of the 
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excitations: The operators dt; and dc] obey the canonical Bose commutation relations. 

Nonlocality of the dipole-dipole interaction results in an integrodifferential character of the 



Bogoliubov-de Gennes equations (25) and (|26|), making it hard to analyze them both analytically 



and numerically. A simpler way is to study the spectrum of small perturbations around the ground 



state solution of the time-dependent GP equation ( |15| ) (see Ref.'^ for this approach to atomic 
condensates). Using this approach in combination with the Gaussian variational ansatz,!^^!^ or 
the Thomas-Fermi approximation,^^ it is possible to obtain analytic results for several low energy 
excitation modes. 

As an illustration, let us consider a Gaussian variational wave function 

V/(r,0= (27) 
= A{t) n expj-i^l— ^ + /r]a,(0 + /T]2j8,(0 

ri=x,y,z { ^^^riV) 

The variational parameters here are the complex amplitude A, the widths Rrj, the coordinates of 
the center of the cloud t]o, and the quantities a,j and /3,j related to the slope and the curvature, 
respectively. The normalization of the wave function to the total number of particles A'^ provides 
the constraint 

= n^/^ \A{t) l^R^RyR^ = const. (28) 
To find the equations governing the variational parameters, we notice that the time-dependent 



GP equation ( [15] ) are equivalent to the Euler-Lagrange equations for the action 

dtL (29) 
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with the Lagrangian 



L = 



1m 



- f/tr(r)] drdv' I V/(r) I V(r - r') | v/(r') 



We therefore can obtain an effective Lagrangian Lgff that depends on the variational parameters by 



inserting Eq. ([28]) into Eq. (30) and integrating over space coordinates. We obtain 



-hN(p + 



+ ■ 



Ri 



+ d Jrexp 



Bcos^e 



(30) 



where (p is the phase of A [the modulus of A was excluded by using Eq. (28 1] and we set t]o(0 = ^ 
and OCri{t) = for simplicity (this corresponds to ignoring the so-called sloshing motion of the 
condensate). The standard Euler-Lagrange variational procedure 



d dL 



dL 



dt \ dqj I dqj 







with qj =(Rri, Prj) provides equations of motion for the parameters Rj^, and /3,j : 



niRrj 
2hRr 



and 



R11 + 



dR 



-U(R)^0. 



(31) 
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Figure 6: Low-energy excitation modes of a dipolar condensate. 

The above equation describes the motion of a particle with a unit mass in the potential 

f fi^ (0^R\\ N 
^^^^ ^ \ [im^ ^ ~ ) + 4V2ny^R,RyR, 

Therefore, the frequencies of small amplitude oscillations around the stationary solution can be 
read from the second derivatives of the potential U (R) at its minimum. In this way one can obtain 
the frequencies for the first three compressional excitation modes. In Refs.l^and, ^ these frequen- 
cies and the corresponding shapes of the cloud oscillations were found for a cylindrical symmetric 
trap, Rx = Ry = Rp, Rz = LRp, see Fig. |6] In the considered cylindrical geometry with dipoles 
are oriented along the z-axis, the projection M of the angular momentum on the z-axis is a good 
quantum number that can characterize the mode: One has M = for modes 2 and 3 and M = I 
for mode 1. The modes 2 and 3 as often called breathing and quadrupole modes, respectively, 
and we will follow this convention here. (In the Thomas-Fermi approximation, one can find an- 
alytical expressions for these modes, see Ref.l^ Important is that with increasing the strength of 
the dipole-dipole interaction, the quadrupole mode 3 demonstrates the tendency towards instabil- 
ity, and becomes unstable when e^^ reaches some critical value. This character of instability via 
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softening of the mode 2 is similar to that in a Bose gas with a short-range attractive interaction 
(a < 0). 

The situation for a dipolar gas with dominant dipole interactions is more complicated,'^,'^.'^ 
It was found (see Ref.!^ that the instability of collective modes of a dipolar EEC reminds that 
of a gas with an attractive short-range interaction only if the trap aspect ratio is larger than the 
critical one, / ^ (numerically was found I > 1.29): The lowest frequency "breathing" mode 
2 becomes unstable when the parameter o = Na^/lp — t- cTc- The variational approach discussed 
above provides the scaling behavior of its frequency £02 near the critical point (see Ref.'^: (02 ~ 
(c7c — a)^, with /3 = 1/4, which is very close to the experimental value /3 ~ 0.2 for Chromium 
BEC.El 

For intermediate values of I above (0.75 <l < 1.29), the mode which drives the instability 
(the lowest frequency mode) is a superposition of breathing and quadrupole modes with the expo- 
nent /3 still close to 1/4. The mode has the breathing symmetry (mode 2) for o far below Oc, while 
it changes and becomes quadrupole-like (mode 3) as o approaches the critical value Oc- 

For / close to (Z < 0.75) the lowest frequency is the quadrupole mode 3. The frequency of 
this mode Gh, tends to zero as o approaches the critical value, Gh, ~ {Oc — o)^ , with the exponent 
/3 ~ 1/4 if / is not too close to l^. When / approaches one has Oc °°, and /3 > 1/4. Finally, 
when / = the frequency of the lowest frequency quadrupole mode £03 can be zero only for 
(J = oo.lMl(Note that this result cannot be reproduced within the Gaussian variational ansatz, which 
in general does not provide reliable results close to the instability, see Ref.l^ 

Collective modes for the case / < /* were analyzed in Ref.'^ on the basis of the Bogoliubov- 
de Gennes equations (25) and ( [26] ). The two possible type of solutions for the stable condensate 
were already mentioned above: A pancake (normal) shaped condensate (the maximum condensate 
density is in the center of the trap), and a biconcave shaped condensate (the maximum condensate 
density is at some distance from the center of the trap). It was found that in the case of a pancake 
condensate, the mode which drives the instability has zero projection of angular momentum on the 
z-axis, M = 0, and consists of a radial nodal pattern. The number of the nodal surfaces increases 
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with decreasing / (flattening of the condensate). This "radial roton" mode in a confined gas can 
be viewed as an analog of the roton mode in an infinite-pancake trap from Ref.,'^ see below. In 
a biconcave condensate near the instability, the lowest frequency mode has non-zero projection 
of the angular momentum on the z-axis, M ^ 0. This mode is an "angular roton" in the trap: 
For a biconcave-shaped condensate, the maximum density is along the ring, and an angular roton 
corresponds to density modulation along this ring. The instability in this case corresponds to the 
collapse of the condensate due to buckling of the density in the angular coordinate, and, therefore, 
breaks the cyhndrical symmetry spontaneously (see Ref.l^for more details). 

3.2.3 Roton instability of a quasi 2D dipolar condensate. 

Let us now discuss the effects of the long-range and anisotropic character of dipole-dipole forces 
in the physically simpler case of an infinite pancake shaped trap, with the dipoles perpendicular to 
the trap plane.!^ It was found that a condensate with a large density hq can be dynamically stable 
only when a sufficiently strong short-range repulsive interaction is present. Otherwise, excitations 
with the certain in-plane momenta q become unstable when the condensate density hq exceeds 
the critical value Hc- Interestingly, the excitation spectrum of a stable condensate with the density 
riQ < ric has a roton-maxon form similar to that in the superfluid helium (see also Ref.'^ for the 
quasi-2D version of this problem). 

The time-dependent GP equation for the condensate wave function y/(r, t) of dipolar particles 
harmonically confined in the direction of the dipoles (z-axis) reads 



2m 



? ? ? 

V^-^ + —^+g\xf/{r,t)\^ + 



+ j dr'Vd{r~r')\^f{r',t) 



V/(r,0, (33) 



where ft)^ is the confining frequency. Let us assume the ground state to be uniform in the in-plane 
directions such that the ground state wave function V^o(z) is independent of the in-plane coordinate 
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p = (x^y). We can then integrate over p' in the dipole-dipole term of Eq. (33 1 with the result 



¥o{z)=0, 



(34) 



where = Snd^/S. This one-dimensional equation is similar to the GP equation for short-range 
interactions. The simplest case corresponds to g+gj > 0, where the chemical potential jU is always 
positive. Let us consider the case ji ^ hcOz (large condensate density) such that we can use the 
Thomas-Fermi approximation to find the condensate density profile in the confined direction: 



noiz) = Voiz)= noil -z^/L") 



where nQ = ii/{g + ga) is the condensate maximum density and L = (2/i / m(o|) is the Thomas- 
Fermi size. 



Eq. ( 33 1 can now be linearized around the ground state solution i//o (z) to obtain the Bogoliubov- 



de Gennes equations for the excitations. These equation are Eqs. (25) and (26) with ?7tr(r) = 
m(o|z^/2 and V^o(r) = M^o{z). Having translational symmetry in the in-plane directions, we can 
characterized the solutions of these equations by the momentum q of the in-plane free motion. In 
addition to q, we also have an integer quantum number j > related to the motion in the z-direction 
such that the amplitudes {M(r),v(r)} have the form {M(z),v(z)}exp(zqp). After introducing the 
new functions /± = m ± v, the Bogoliubov-de Gennes equations readf^ 



(Of- — ^kin/+: 

(0f+ = H^J_+H,,,[f. 



(35) 
(36) 



where 



2m 



-— y + <? H 
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is the kinetic energy operator and 



H,M-] = ^{g + gd)wli.z)f-{z)+ (37) 

oo 

- 2^gd¥o{z) J dz'\iroiz)exp{-q\z-z\)f-{z) 



is the interaction operator. The solution of the above equations provides excitation frequencies 
(Oj{q) which depend on both j and q. The most relevant for the stability analysis is the lowest 
frequency branch COo{q) for which the confined motion is not excited in the limit q ^ 0. 

Because of nonlocality of the dipole-dipole interaction, an effective coupling [the second term 



in the right-hand- side of Eq. (38)] becomes momentum dependent. For small in-plane momenta 



qL<til, excitations of the lowest branch are essentially two-dimensional with a repulsive effective 
coupling, and their spectrum has been found in Ref .'^ These excitations are phonons propagating 
in the jcy-plane with the sound velocity c^: 



COo{q) = Csq, Cs = {111 /3m 



,1/2 



In the opposite limit of large in-plane momenta qL^ I, the excitations are three-dimensional 
and the interaction term is then reduced to 



Eqs. (35) and (36) are then equivalent to the Bogoliubov-de Gennes equations for the excita- 
tions in a condensate with a short-range interaction with the strength 2g — gj. This interaction 
is repulsive if the parameter /3 = g/g^i > 1/2, and all excitation frequencies in this case are real 
and positive for any in-plane momentum q and condensate density hq. In the other case /3 < 1/2, 
the interaction is attractive resulting in dynamical instability of a condensate with regard to high 
momentum excitations at a sufficiently large density. 



The analysis in the Thomas-Fermi regime of the system of equations (35 ) and (36 ) in the most 
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Figure 7: Dispersion law COo{q) for (a) j8 = 1 /2, jl/hcOy = 343; (b) j8 = 0.53, n/hco^ = 46 (upper 
curve) and j8 = 0.47, /i/^COz = 54 (lower curve). The solid curves show numerical results. (Taken 
from Ref 

interesting case qL I and /3 close to the critical value 1/2 was performed in Ref.'^ It was 
found that at the critical value /3 = 1/2 the momentum dependence of the excitation frequencies is 
characterized by a plateau [see Fig.|7]^a)], and the j-th branch reads 



coj (q) = + h^(0^[l+ jU + 3) /2] , » 1 , 



where = h^q^/2m. 



For j8 7^ 1/2, the lowest branch of the spectrum is 



where the condition /i£^|2/3 — 1|/(1 +/3) ^ fP'ol was assumed. Eq. (38) provides us with two 
types of behavior of the lowest-frequency mode (OQ{q). It is either monotonously increases with q 
[see Fig. [TJ^b)] when /3 > 1/2, or has a minimum if /3 < 1/2. Being combined it with the fact that 
(Oo(^) grows with ^ for ^ 1 , the existence of this minimum results in a roton-maxon character 
of the spectrum as a whole [see Fig. |7]^b)]. This type of the excitation spectrum in an infinite 
pancake trap can be understood as follows: For small in-plane momenta q <^ excitations have 
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two-dimensional character and are phonons because dipoles being oriented perpendicular to the 
plane of the trap, repel each. On the other hand, excitations with large momenta q ^ have 
three-dimensional character and, hence, the repulsion between them is reduced. The excitation 
frequency therefore decreases with an increase of ^,reaches a minimum, and starts to increase 
again as the excitations continuously enter the single-particle regime. 



The roton minimum for j8 close to 1/2 found from Eq. (38 1 is located at ^ = {16il5/15h(0^y/^/h 
where 5 = 1/2 — /3, and L = ^Jhjmxa^ is the harmonic oscillator length for the confined motion, 
and correspond to the excitation frequency 



«Omin= V^^W^ -(8/i5/15)2. 

This minimum becomes deeper with increasing the density (chemical potential) or 5, and reaches 
zero at ^ = \/2/lz for iid/fico^ = 15/8. Excitations for larger values of }i5 /fico^ have imaginary 
frequencies for <5r ~ and, therefore, the condensate becomes unstable. 



Eqs. (35) and (36) for various values of j8 and ix/fico^ were solved numerically in Ref.'^The 
results for the excitation spectrum in the Thomas-Fermi regime are shown in Fig.|7} demonstrating 
a good agreement between numerical and analytical approaches. 

For non- Thomas-Fermi condensates, the stability does not require as strong short-range repul- 
sive interaction as in the Thomas-Fermi regime because of a large kinetic energy in the confined 
direction. The spectrum of excitations in this case also has a roton-maxon character, although the 
appearance of the roton minimum and the instability take place at smaller values of /3, see Ref.'^ 
for details. 

Up to now, a roton-maxon dispersion was observed only in liquid He with strong interparticle 
interactions. Dipolar condensate provides the first example of a weakly interacting system with 
a roton-maxon excitation spectrum. This spectrum can be controlled and manipulated by chang- 
ing the density, the strength of the confinement, and the short-range interaction, starting from the 
Bogoliubov-type spectrum, then creating the roton minimum, and finally reach the point of insta- 
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bility. 

It is important to point out that the existence of the roton minimum with (^r 7^ at a given 
/3 < 1/2 for jl/hco^ just below the point of instability is likely to indicate the existence of a new 
ground state presumably with a periodic density modulation. This is in contrast to the instability of 
condensates with attractive short-range interaction, which is driven by unstable long wavelength 
excitations resulting in local collapses. In this case the chemical potential is negative and not 
bounded from below such that no new ground state exists. In the Section [6] we will show that the 
excitation spectrum of a two-dimensional dipolar gas in the strongly interacting regime n2a^ > 1 
also has the roton minimum, and the system undergoes a liquid to solid quantum phase transition. 

4 Weakly interacting dipolar Fermi gas 

In this Section we discuss fermionic dipolar gases in the weakly interacting regime. Most of the 
discussion will be devoted to a single-component (polarized) dipolar gas with only brief mention- 
ing some results available for two- and more component dipolar systems. 

The crucial differences in the behavior of many-body fermionic systems as compared to bosonic 
ones are related to the Pauli principle: identical fermions are not allowed to be in the same quan- 
tum state. As a result, the many-body wave function of a single component Fermi gas should 
be antisymmetric with respect to permutations of the positions of any two particles. In the second 
quantization, this requires that the field operators i/A(r) and xfr^ obey the canonical anticommutation 
relations 

{ i?r(r) , i^/t (rO } = V>(r) v/^ (r') + y/t (r') v/(r) = 5 (r - r') , 
{VA(r),VA(rO} = {v/^(r),VAt(rO}=0. 

As a direct consequence, the wave function of a relative motion of two identical fermions is allowed 
to have components with only odd values, L = 1,3, . . ., of the angular momentum, and vanishes 
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when the interparticle distance tends to zero. Therefore, the low-energy scattering of two identical 
fermions is insensitive to the short-range part of their interaction and is solely determined by the 
long-range dipole-dipole part Vj. As a result, for a single-component polarized dipolar Fermi gas 



Another consequence of the Pauli principle is that the state of a many-body system of fermions 
at a low temperature T is completely different from that for bosons. The average number of ideal 
fermions in a quantum state i with the energy £, is given by the Fermi-Dirac distribution 



where /i is the chemical potential, which depends on T and, as usual, ensures the fixed total 
number of particles A^^ = The ground state {T = 0) therefore corresponds to all quantum 

state with e, < /i(r = 0) = being completely completely occupied [n(e,) = 1], while the states 
with e, > £f are are being empty [n(e,) = 0]. The energy £f is called the Fermi energy and sets the 
typical energy scale a many-body system of fermions. 

The ground state of an ideal homogeneous Fermi gas with £p = /2m corresponds to the 
so-called Fermi sphere: All quantum states with momenta p below the Fermi momentum pf — 
\flm£p are occupied and the states with p > pf are empty. The states with momentum p = pf 
form a surface in the momentum space called the Fermi surface, which separates the filled and 
empty states. Semiclassical state counting provides the relation between the Fermi momentum pf 



we can omit the contact term in Eq. ([2]), and the corresponding Hamiltonian then reads 



H 




(39) 



where is given by Eq. ([T]) and i7tr(r) is the trapping potential (if present). 
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and the density n of a single-component homogeneous gas: 

For a trapped Fermi gas we can establish the similar relation but between the local Fermi 
momentum Pf{^) and the local density of the gas n(r), 

where 



PF{r) = y/2m[^-U^,{r)], (42) 

provided the chemical potential /i is much larger than the level spacing in the trapping poten- 
tial f/tr(r)- This condition corresponds to a large number of particle in the trap, most of them 
occupying high energy states of the trapping potential. The wave functions of these states are qua- 



siclassical (see, for example,"^, and the calculation of the gas density results in Eq. 41 which is 
the essence of the local-density (Thomas-Fermi) approximation. This approximation is legitimate 
when the trapping potential changes slowly over the distances of the order of the average inter- 



particle separation h/pf. For an ideal Fermi gas in a harmonic potential expression ( |42[ ) 

gives 



PF(r) = ^Im \\i - [/tr(r)] 

\ 



/ z?" 

V ^TFx '^TFy ^TFz 

where Pf{0) = \/lmjx is the Fermi momentum in the center of the trap and Rtfu = A/2/i/m(o2 is 
the Thomas-Fermi size of the gas cloud in the a-direction. The density of the gas in this approxi- 
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mation, according to Eq. ([4T]), reads 



/ 2 2 2 \^/^ 

nMr)=no[l-^, ^ , (43) 

\ ^TFx ^TFy ^TFz J 

where hq = (2m/i)-^/^/ (6;r^^^) is the density in the trap center. The calculation of the total number 
of particle with the use of the above density distribution relates the chemical potential ji to the total 
number of particle and the parameters of the trap: 



II = fia){6N) 



1/3 



where w = (cOxCOy.ft)^)^/-^. 

For understanding the properties of the fermionic systems it is important to keep in mind that 
the ground state in the form of a filled Fermi sphere stores a large amount of kinetic energy. This 
guaranties applicability of the perturbation theory for dilute dipolar systems with pfUd/fi ^ 1. 
Another consequence is the improved stability of fermionic dipolar gases, as compared to the 
bosonic ones, against collapse due to the attractive part of the dipole-dipole interaction. This can 
be understood as follows: The energy per volume for a homogenous dipolar Fermi gas with density 
n effectively attractive two-body interaction can be written as 

^(") 3 2 2 

— ^ = -efH-Adn, 

where the first term is the kinetic energy of the filled Fermi sphere and the second term is the 
interaction energy with some numerical coefficient A of the order unity. The first term scales as 



^2/3 [see £q (40)] and provides an energy barrier between states with small n and positive energy 
and collapsing states with n 00 and negative energy. Therefore, one expects the stability against 
collapse when the system is dilute: n^^^^ ^ aj or, equivalently, pfUd/h ^ 1, and instability in 
the dense system with n^^^^ < when the interaction energy becomes comparable or larger that 
the kinetic energy. Applying this arguments to a trapped single-component dipolar Fermi gas, one 
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expects to have a stable gas for pp{0)ad/h < 1 or N^/^a^/l < 1, where / is the oscillator length 
and we use Eq. (43) to obtain ~ [pf{0)l/h]^. We provide more details on the issue of stability 
later. 



4.1 Effects of dipole-dipole interactions. 

When considering effects of interparticle interactions in Fermi systems, one has to keep in mind 
two possible scenarios depending on whether the properties of the system (the ground state and 
excitations) change continuously or abruptly when interactions are switched on. In the first case 
an interacting system is called normal Fermi liquid (in other words, belongs to Fermi liquid uni- 
versality class) and has properties that are very much similar to those of an ideal Fermi gas. Of 
course, the interaction leads to appearance of new features (collective modes, for example), which 
are absent in a non-interacting gas, but for many applications the system can be considered as 
an ideal gas of fermionic non-interacting quasiparticles. For weak interparticle interactions, the 
properties of the interacting system can be obtained with the help of perturbation theory starting 
from the non-interacting Fermi gas. In the second scenario, the ground state and excitations of 
interacting system are qualitatively different from those of the non-interacting Fermi gas, and a 
system is in non-Fermi liquid universality class. This scenario is usually associated with breaking 
of some symmetries of an ideal gas: phase (or gauge) symmetry in a superfluid Fermi liquid or 
translational symmetry in a charge-density wave or crystal state. The new ground state cannot 
be continuously connected with the filled Fermi sphere (ground state of a non-interacting Fermi 
gas), and, therefore, one has to go beyond simple perturbative expansions to describe those states. 
It is important to mention that one does not necessarily need a strong interaction for the second 
scenario. For example, even an infinitesimally small attractive interaction results in a superfluid 
ground state. The smallness of the interaction in this case manifests itself in low (much smaller 
that Tf) critical temperature Tc - the temperature above which the superfluid properties disappear 
and the system returns to normal Fermi liquid. In contrast, the charge-density wave state requires 
strong interaction, and this state disappears at temperature comparable or larger that Tf when the 
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gas is essentially classical. 

As we will discuss below, depending on an experimental setup, both scenarios are possible in 
a polarized dipolar gas: A 3D polarized dipolar gas is in the superfluid state for low temperatures, 
T < Tc <^ Tp, and in the normal state (Fermi-liquid) for T > Tc. The state of a monolayer of 
polarized dipoles depends on temperature and relative angle between the dipole moments and 
the motion plane of molecules. For the perpendicular orientation of dipoles, the gas is in the 
normal state, but, starting from some critical tilting angle, becomes a superfluid at small enough 
temperatures T < Tc. In both cases, the increase of the strength of the dipole-dipole interaction 
leads to the instability of the homogeneous state resulting to a collapse or formation of density- 
wave state with broken translational symmetry. 

4.2 Normal (anisotropic) Fermi liquid state. 

We begin with discussion of a normal Fermi liquid state of a dipolar fermi gas, which is a generic 
state for a fermionic dipolar gas at finite (Tc<T < Tp) temperatures, as well as for a purely repul- 
sive (in a monolayer, for example) dipolar gas, in a weakly interacting regime kpud < 1. Following 
the original idea of Landau, an interacting normal Fermi system (Fermi liquid) can be described in 
terms of fermionic quasiparticles, which can be viewed as particles together with disturbances they 
produce in the system due to interactions with another particles (particles surrounded by particle- 
hole excitations) - dressed particles. In the ground state, the quasiparticles occupy all states with 
energies smaller or equal than the chemical potential ep forming a filled Fermi sphere (in a 
spatially uniform dipolar gas this corresponds to a deformed Fermi sphere in momentum space due 
to anisotropy of the dipole-dipole interaction, see below). Excited states are obtained by moving 
some quasiparticle from occupied states below jU to empty ones above - creation particle-hole exci- 
tations. The advantage of this description is that weakly excited states correspond to small number 
of particle-hole excitations near the Fermi surface and, hence, can be described using the dilute gas 
approximation. Note that, although we are talking about filled quasiparticle states inside the Fermi 
sphere, quasiparticles in the Fermi liquid are well-defined only in the vicinity of the Fermi surface 
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where their energies e(p) are much larger than the inverse of their Ufe-times Tp due to decay via 
creation of particle-hole pairs. (In a weakly interacting gas the quasiparticles are well-defined for 
all momenta.) This is because the presence of occupied states below the Fermi energy strongly 
reduces the phase space volume for such processes, and, as a result, the life-time Tp of quasipar- 
ticle near the fermi surface in a Fermi liquid are much larger than the corresponding time Tc in 
a classical gas with the same interparticle interactions and density, Tp ~ [sp / e{Tgt)Yxc Tc. But 
those are quasiparticles we actually need to describe low-energy excitations of the Fermi system 
and its behavior at low temperatures and under weak external perturbations. 

The change of the quasiparticle distribution 5n^ (we assume here a spatially homogeneous gas) 
results in the change of the energy of the system 

dE = £[e(p) + ju]5np + £/(p,p')5«p5v: (44) 
P P,P' 

where e(p) = 5£'/5«p|g„p=o is the quasiparticles energy (counted from the chemical potential ju) 
and the second term describes the interaction between quasiparticles with /(p, p') = 5^E / 5n^5n^! \ Snp=o 
being the Landau /-function, which plays a crucial role in the Fermi-liquid theory, and can be ei- 
ther calculated perturbatively (if the gas is weakly interacting) or measured experimentally. Note 
that the function / describes the change of the quasiparticles energy under the change of quasipar- 
ticle distribution as a result of their interaction, 

5e(p) = L/(p>p')5v' 
p' 

which gives rise to Fermi-liquid corrections and make possible collective motion of quasiparticles 
(collective modes) even when collisions between quasiparticles can be neglected. 

For states close to the Fermi surface (the boundary between occupied and empty states), the 
quasiparticle energy has the form 

e{p)^^{p-PF), 
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where pf is the Fermi momentum specifying the Fermi surface in momentum space, and m* is 
the effective mass. The Fermi momentum jc/- is related to the density in the same way as in the 
ideal gas, Eq. ( [40] ), reflecting the fact that numbers of particles and quasiparticles are equal, while 
the effective mass m* can be expressed in terms of the /-function (see, for example Ref.®|). The 
compressibility K = rT^dil/dn, where = dE / dn is the chemical potential, is another impor- 
tant quantity, which can also be expressed in terms of /-function. For a stable system one must 
have K > 0. Therefore, the knowledge of the compressibility as a function of system parameters 
provides us with stability conditions of the system against collapse. The stability of the system 
against possible deformations 5np of the Fermi surface around its equilibrium form (Pomeranchuk 
criterion,'^i'2Q|]) ^^j^ j-,g obtained from the requirement that the change of the energy caused by this 



deformation, Eq. (44), is positive. In this way one can detect instabilities different from collapse, 
related to the nonuniform change of the Fermi surface. 

The /-function determines also collective modes in the Fermi liquid (Landau zero sound), 
which correspond to a coUisionless coherent dynamics of particle-hole excitations. The simplest 
way to describe zero sound is to use a semiclassical (or Wigner) quasiparticle distribution function 
n(r,p,?), which is the Fourier transform of a single -particle density matrix with respect to the 
relative coordinate, 

n(r,p,0 = J ^^r'^v/^(r + r72,f)v/(r-r72,0)exp(-/pr7^), 

and describes the local momentum distribution of particles at position r. In the ground state of 
a spatially homogeneous system, no(r,p,?) = 0{pf — p) corresponds to a filled Fermi sphere. 
For a thermal equilibrium state, the step function 0{pF — p) has to be replaced with the Fermi- 
Dirac distribution, n7(r,p,?) = [exp(e(p)/r) + Ij^^Time evolution of non-equilibrium distribu- 
tions n(r, p,?) = neq(p) + 5n(r,p,f) are described by the quasiparticle kinetic equation 

, d de d de d ^ , , „ 
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where e = e(p) +Y.p' fi'Pi'P')^'^i^i'PJ) the collision integral, which normally appears on the 
right-hand side, is set to zero assuming low temperatures, as discussed above. The solutions of 
this equation of the form 5n(r,p,?) ~ x(p)exp[/(kr— cot)] <^ neq(p) with CO = ck and c > are 
called Landau zero sound and described coherent motion of particle-hole pairs - propagation of 
a deformation of the Fermi surface. Generically, the solutions of this kind exist when /(p,pO is 
positive (for more details and exact criterion see, for example, Refs."^. Note that the condition 
c > vf separates the zero- sound from the continuum of particle-hole excitations and ensures its 
long life-time. In the opposite case the energy of zero-sound waves would be inside the continuum 
of particle-hole excitations and, hence, the waves would rapidly decay into incoherent particle-hole 
excitations (Landau damping). 

4.2.1 Anisotropic Fermi surface and single-particle excitations 

Due to anisotropy of the dipole-dipole interaction, the Fermi surface in a dipolar gas is not a sphere 
any more and the modulus of the Fermi momentum depends on the direction. The effective mass 
becomes a tensor that can be defined from the relation between the Fermi momentum p/r and 
fermi velocity \f = (9e(p) / (9p|p=p^, pfi = m*jVfj. This can easily be seen by using the following 
variational ansatZjI^^HS 

n(p) = e[p^F- UpI + p]) - o'pI] (46) 

and find the variational parameter a by minimizing the total energy of the system with the interac- 
tion energy calculated in the Hartree-Foch approximation: 
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where V is the volume of the system and only exchange (Fock) terai contribute to the dipole-dipole 
interaction energy because the direct Hartree contribution 



Ed = drdr'n{r)VAir-r')n{r') 



where n{r) is the gas density, vanishes in a homogeneous gas as a result of angular integrations. 
It was found that /3 < 1 so that the Fermi surface is deformed into a spheroid stretched along the 
direction of dipoles (prolate spheroid). These findings were supported by microscopic calculations 
in the spirit of Landau liquid theory in Refs.l^ and.!^ The quasiparticle energy calculated from 



Eq. (47 1 reads 

f dp' 



^"') = l^-/(dlp''^'>'->'>"'''-''' 



which corresponds to the following Landau /-function for a spatially homogeneous gas 

/°-(p,pO = -yd(P-p') (48) 

with only exchange contribution. Note that the condition of spatial homogeneity of the gas is 
essential for validity of Eq. ( [48] ). This is because the Fourier component of the dipole-dipole 
interaction Vd(q) in non-analytic for q — (the limit depends on the direction q approaches zero). 
As a result, the direct (Hartree) contribution vanishes only in the spatial homogeneity gas, in which 
one has Vd{q) averaged over the direction of q, which is zero. In an inhomogeneous gas, this is not 
the case and one also has the contribution of the direct dipole-dipole interaction, see, for example. 



Eq. (52) describing spatially inhomogeneous variations of the quasiparticle distribution. 

Setting e(p) to zero gives the position of the Fermi surface in momentum space p/^ = npF{n), 
where n is a (radial) unit vector (direction) in momentum space. The chemical potential jl then 
has to be defined selfconsistently from assuming a fixed gas density n = kp/67r^: 

(■pF{n) 
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For weak interaction one 



/>f (n) = hkp 



1 , 

I + —adkf{?> co^^ On - I] 
9n 



where 6^ is the angle between n and the z-axis. This gives j8 = 1 — la^kf /9k. The energy and the 
chemical potential are 



£ _ 3 n\l 



1 



81;r' 



-{adkf 



and 



2m 



1 



28 



;{adkf 



405 ;r2 

After expanding the quasiparticle energy in the vicinity of the Fermi surface, e(p) ~ V/r(p — p/r), 
one finds that the tensor of the effective mass has only longitudinal m2(n) and and transverse 
my(n) components: 

\F 



n . . . +ee- 



where is the polar angle unit vector and 



m 



m2(n) 
m 



1 — —a^kf (3 cos^ 0n — 1 ) 



1 



my(n) 3;r 



6;r 



adkf sin29n- 



Calculations for moderate strengths of the dipole-dipole interactions and finite temperatures were 
performed in Refs., 1221,1241961 including the trapped case (Ref.'^, as well as 2D (monolayer) and ID 
(tube) gases, and two-component dipolar gas (Ref-I^. We mention here only some details for a 
monolayer and refer to these references for more details. 

In a 2D dipolar gas (monolayer), when the chemical potential ji is much smaller than the 
frequency (O^ of the transverse confining potential in the z-direction, ji <^ co^, the transverse wave 
function of particles is limited to the ground state wave function ^o{z) of the harmonic oscillator, 
such that i/A(r) = \i/{p)^o{z), where p = (x^y) = (p cos (pp,p sin (pp) is the in-plane vector. The 
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corresponding effective 2D dipole-dipole interaction for the in-plane motion 



dzdz!Uzfv^{p,z-7!)Uzf 

^ {2P2(cos0)^(1/2,O;pV2/2)+ 



1 



X 



P2(cos0) — -sin 0cos2(j!)p 



^(1/2,-1;pV2ZJ)} 



X 



(49) 



where ^{a,b;z) is the confluent hypergeometric function and 9 is the angle between the direction 
of dipoles (in the (;c,z) -plane) and the motion (x,y) -plane, has the following Fourier transform 



d 



PJ' 



(50) 



where w{x) = xexp{x^)er{c{x) with erfc(x) being the error function and 0p is the angle between p 
and the x-axis. For p ^ Pf ^ h has 



(P) ~ -TT— p[2P2(cos0) -sin^0cos2(pp] 



(51) 



which is linear in p. (Strictly speaking, expression (51 1 contains also a constant which depends on 



the regularization of the Fourier integral at the origin. This constant corresponds to a short-range 
inerparticle interaction and, hence, has no physical effect in a single component Fermi gas because 
all its contributions should vanish upon proper antisymmetrization. We therefore set this constant 
to zero.) Within the Hartree-Fock approximation one then obtains (assuming adkp ^1) 



Pf{t^) = ^kf 



8 2 

1 -I- —^adhp sin 9 cos 2(j!)p 



for the position p = np/r (n) of the Fermi surface, where now n = (cos (pn, sin ^n) is the unit vector 
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of direction in the -plane, and 



1 + [P2 (cos 9)-— sin^ 9 cos 2^p] 



mlDii^) 3;r 10 

-adkp sm 0cos2(j!)p 



for the longitudinal (along n) and transverse (perpendicular to n) components of the effective 
mass, respectively, that for 0=0 (dipoles are perpendicular to the plane and, therefore, the 

system has rotational symmetry around the z-axis) the deformation of the Fermi surface disappears 
and Yp = PF/fn* with m* /m^ 1 —Aadkf/3n. 

4.2.2 Collective modes (Landau zero sound) 



Collective modes in a dipolar gas can be studied on the basis of the kinetic equation, Eq. (45 ). For 
small deviation of the quasiparticle distribution from equilibrium of the form 

5/7(r,p,?) ~ x(p)exp[z(kr-a)?)] <neq(p), 

the kinetic equation reduces to the following equation on the unknown function x(p) 



[a)-kVp£(p)]x(p) = 

"^'^(-'^/(W^C^-"'"''''"'''- '''' 

where Vp = (9/(9p and the wave vector is assumed to be much smaller than kf = pf/h,k<^kf. In 
the Hartree-Fock approximation, the /-function is /(k,p,p') = Vdlk) — Vd(p — p'), see comments 
below Eq. (48) and Ref s., 1251.123 gq. (52) was analyzed numerically for a 3D gas in Refs.l^S andl^^, 



and in 2D gas in Refs.,'2a.E21(Note that in a 2D gas - monolayer - the first term in the /-function can 
be omitted following the arguments from the end of the previous section.) The results for the sound 
velocity in a 3D gas together with the propagation limit due to particle-hole continuum are shown 
in Fig. [8] as a function of the propagation angle (the angle between k and the z-axis). They show 
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Figure 8: Speed of sound (solid curve) in a 3D dipolar Fermi gas as a function of the angle of 
propagation 0^ relative to the direction of polarization, for dipolar interaction strengths adkf = 1 
(left panel) and aakp — 3 (right panel). The speed is measured in units of v^- The dashed curve 
represents a lower bound on the speed of any undamped mode. (Taken from Ref.'^ 

that the collective mode can propagate only in the certain cone of directions around the direction of 
dipoles polarization. This is counterintuitive to some extend because this is the direction in which 
two dipoles attract each other; on the other hand in the perpendicular direction, in which dipoles 
repel each other and, hence, one would expect the existence of the collective mode, no zero sound 
is possible. This can be understood by noting that the zero-sound propagation is dominated by the 
exchange contribution, not direct one, and, therefore, the intuition based on the direct interaction 
does not work, see Ref.l23for more discussions of this issue. The sound velocity depends strongly 
on both the propagation direction 0^ and interactions strength ajkf: It increases monotonically 
with 0k for ajkf < 1 and becomes a non-monotonic in 0k for a^kf > 1, see Fig. [8} 

In a dipolar monolayer (quasi-2D gas), the situation is even more intriguing because the exis- 
tence of zero-sound and the value of the sound velocity strongly depend on the propagation direc- 
tion ( (^k is the angle between k and the x-axis), on the tilting angle 6, and on the strength of the 
interaction. In this case, there is no collective modes if the tilting angle is smaller that some critical 
value that depends on the interaction strength, see Fig. |9j This is again counterintuitive because 
the direct interaction for small tilting angles is purely repulsive (dipoles are almost perpendicular 
to the motion plane), and one would expect stable collective zero-sound modes. However, similar 
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Figure 9: The speed of zero sound in a 2D dipolar Fermi gas as a function of the tilting angle Q 
for adkp = 1 and the propagation angles ^ = (upper panel) and ^ = n/2 (lower panel). The 
shaded region correspond to strong damping. (Taken from Ref.l^ 
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to the 3D case, the /-function contains only exchange contribution, and this explains such peculiar 
behavior of collective modes. (The collective modes without exchange contribution were consid- 
ered in Ref.l^^ Note also that with increasing the tilting angles from the critical one, the directions 
in plane, in which one has propagating zero-sound, changes from those around the projection dy 
of the dipole moment d on the plane to those around the direction perpendicular to dii (when the 



polarization of dipoles approaches the plain), see Fig. 10 



Note, however, that the above results for collective modes were obtained with the /-function 
in the lowest-order Hartree-Fock approximation, /(k,p,p') = Vd(k) — Vd(p — p')- Taking higher 
order terms into account can change the situation: As it was shown in Ref.^^ for the case of a 
2D dipolar Fermi gas polarized perpendicular to the motion plane, 0=0, the inclusion of second 
order contributions to the /-function results is the appearance of a stable zero-sound mode with 
the velocity s = Vf[l + 2{adkf)^]. 



4.3 BCS pairing in a homogeneous single-component dipolar Fermi gas. 

The partial attractiveness of the dipole-dipole interaction opens the possibility for BCS pairing in 
a fermionic many-body dipolar system at sufficiently low temperatures. As we will see in this 
section, the pairing in dipolar systems has generically an unconventional character (different from 
a singlet isotropic 5-wave pairing as in a two-component fermionic system with an isotropic attrac- 
tive interaction), and a superfluid state has many peculiar properties that are different from those of 
conventional superconductors. Indeed, the 5-wave (together with other even angular momentum) 
two-particle interaction channel is forbidden in a single-component Fermi gas by the Pauli prin- 
ciple. On the other hand, the angular part of the matrix element for the dipole-dipole interaction 
between the states with the angular momentum L = 1 (p-wave channel) and its projection on the 
z-axis M = is negative (i.e. corresponding to an attractive interaction): 

(L= 1,M = O|1-3cos20|L= 1,M = 0) = -y < 0, (53) 
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Figure 10: The speed of zero sound in a 2D dipolar Fermi gas as a function of the propagation 
angle ^ for adkf = 1 and 6 = n/4 (upper panel) and 6 = n/2 (lower panel). The shaded region 
correspond to strong damping. (Taken from Ref.!^ 
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and, therefore, can lead to BCS pairing. (The matrix elements between the states with M = ±1 are 
positive.) It easy to see that this pairing should be anisotropic, reaching its maximum amplitude 
in the direction of dipolar polarization when two dipoles attract each other, and being zero in the 
perpendicular directions corresponding to repulsive dipole-dipole interaction. As we will see, the 
dominant contribution has />-wave symmetry. 

The Cooper pairing in a polarized single-component dipolar Fermi gas has been discussed in 
Refs.'^and'^ within the BCS approach with the restriction to purely p-wave pairing. An exact 
value of the critical temperature and the angular dependence of the order parameter for a dilute gas 
were found in Ref.li^ 

After omitting the contribution of the short-range part of the interparticle interaction, as dis- 
cussed above, the Hamiltonian of a homogeneous single-component polarized dipolar Fermi gas 
reads 



2m 



H = Jdr\if\r) 

+ ydrdr'xif\r)xif\r')Vd{r-r')xjf{r')f{r). (54) 

We considered the property of the system with this Hamiltonian in the dilute limit na^ ^ 1 
and at temperatures T much smaller than the chemical potential /i (or the Fermi energy Cp), 
T II = £f ^ {h^ /2m){6Tt^nf-/^ , relevant for Cooper pairing. In this case one can neglect 
the corrections 5/1 to the chemical potential due to the dipole-dipole interaction because 5jU ~ 
d^n ~ £F{na\y/'^ < Ep. 

The BCS pairing corresponds to a nonzero value of the order parameter 

A(ri -r2) = Vd(ri-r2) {w{r\)W2)) , 

which can be viewed as a wave function of Cooper pairs, because of anticommutativity of the 
fermionic field operators, A(ri — ¥2) changes sign under the exchange of particles ri < — )■ r2 form- 
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ing a pair. As a consequence, the order parameter in momentum space 



^(p) = J drexp{-ipr/h)A{r) 



is also antisymmetric, A(— p) = — A(p). Because of the anisotropy of the dipole-dipole interaction, 
the angular momentum L of the relative motion of two particles is not a conserved quantum number, 
but its projection on the z-axis (in the considered geometry) does. We can therefore write A(p) in 
the form 

A(p)= Az,(p)Y^o(p), (55) 

oddL 

where Yim{p) the spherical harmonics and p is the unit vector in the direction of the momen- 
tum p. We keep in the sum only odd angular momentum L following the discussion above and set 



M = in every term. This is because M is a conserved and, following Eq. (53 1, only for M = 
one has an attractive interaction. 

A nonzero order parameter and, therefore, the superfluid properties in the system appear for 
temperatures below some temperature Tc which is the critical temperature of the superfluid transi- 
tion. This critical temperature and the order parameter A for temperatures below Tc can be found 
from the gap equation,ti2^II03] jj^g momentum representation and assume the order param- 

eter to be real a real function of momentum p) 



where ^(p) = a/A2(p) -f- {p^/2m — ji)^ is the energy of single-particle excitations in the superfluid 
gas. The effective interparticle interaction is described by the function y(p,p') = Vd{p — 'p') + 
5y (p, p'). Here Vd(q) is the Fourier transform of the bare dipole-dipole interaction potential Vd(r) : 
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Vd(q) = yj2(3cos2(0q)-1), (57) 

with 0q being the angle between the momentum q and the z-axis, and 5y(p,p') corresponds to 
corrections to the bare interparticle interaction resulted from many -body effects. The effective 
interaction y(p,p') describes all scattering processes in the system which transform a pair of par- 
ticles with momenta p' and — p' into a pair with momenta p and —p. The leading process here is 
the direct scattering of the two particles on each other [the term V^;(p — p')]^ while the many-body 
corrections 5y(p,p') describe processes of higher order in Vj. The leading terms in 5V are sec- 
ond order in (see Ref. ^^^) and correspond to scattering processes in which the two colliding 
particles interact with each other indirectly with an involvement of particle-hole excitations which 
they create in the system (see more details in Ref.liOU). These processes are important even in the 
weakly interacting case (although they are of the second order in the small parameter) because they 



result in the pre-exponential factor in the expression for the critical temperature, see Eq. (59 1. 



The gap equation (56) can be simplified for temperatures just below Tc because for such tem- 



peratures the order parameter A(p) is small and, hence, the right-hand-sides of Eq.(56) can be 
expanded in powers of A(p). The resulting equation 



A(P) 



dp' 

ilnh) 



;^(P,P') 



^K{p')^\lp') 



(58) 



where K{p) = tanh(i^ /2T)/2^ and ^ = /2m — jl, is equivalent to the Ginzburg-Landau equation 
for a spatially homogeneous order parameter. . 

Note that Eq. ( [58] ) always has the trivial solution A = which corresponds to a normal phase of 
the Fermi gas. The Cooper pairing is associated with a nontrivial solution of the gap equation ( |58| ) 
which exists for temperatures T <Tc. To find the critical temperature Tc, it is sufficient to keep only 



the linear term in the square brackets in the right-hand-side of Eq.(58 ) because A — t- for T — )■ 
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The corresponding linearized gap equation allows also finding the momentum dependence of the 
order parameter. The absolute temperature dependent value of A is determined by the nonlinear 



term in the right-hand- side of Eq.( 58 ). 



The result for the critical temperature reads (see Ref.^^for details) 

r.= 1.44e.exp(-^). (59) 



For temperatures T close to Tc, the anisotropic order parameter A(p) on the Fermi surface, p = pf, 
has the form 



A{pFP)^2.5T,J\-^(j>oiP), (60) 



Tc 



where 



(^o(p) = V2sin (^^cos0pj (61) 

with Gp being the angle between the vector p and the z-axis. 

For momenta away from the Fermi surface the order parameter is 

A(P)^-^ / ^Vd{p-PFP')A{PFP'). 

The dependence of the order parameter A(p) on the modulus p of the momentum p for several 
values of the angle Op is shown in Fig. [TT| This momentum dependence of the order parameter 
in a dipolar Fermi gas is in contrast to that for pairing (both s- and jc-wave) due to a short-range 
interparticle interaction, in which case the order parameter is a constant for momenta p ^ h/ro, 
where ro is the range of the interparticle interaction, and rapidly decays p> fi/vQ. 

The anisotropy of the order parameter in the momentum space described by the function 



(^o(n) = -\/2sin[(;r/2)cos0], see Eq. (60), provides another difference from the conventional s- 
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Figure 1 1 : The order parameter A(p, 0) [in units of A 0)] as a function of the momentum p (in 
units of pf ) for various values of the polar angle d. (Taken from Ref.li^ 
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wave pairing, say in a two-component Fermi gas with a short-range intercomponent attractive 
interaction. As a result of this anisotropy, the gap |A(j9irp)| in the spectrum of single -particle ex- 
citations in a dipolar superfluid gas depends on the direction of momentum p: The gap reaches its 
maximum in the directions parallel to the direction of the dipoles (Op = 0, n), while it vanishes in 
the direction perpendicular to the dipoles {Op = n/2). Similar anisotropic is expected in the prop- 
erties of collective excitations, and, as a result, in the response of the dipolar superfluid dipolar 
Fermi gas. 

The vanishing of the single-particle gap at Op = n/l (and arbitrary azimuthal angle (p, i.e. on 
the line on the Fermi surface p = pf) results in the dependence of the specific heat of the gas 
at low temperatures T ^ Aq ~ Tc (the contribution of collective excitations is proportional to T^). 
Note that for the conventional i'-wave pairing, the low-temperature specific heat is determined by 
the contribution of collective modes (oc T^), while the contribution of single -particle excitations is 
exponentially suppressed. 

The above mentioned properties of a superfluid dipolar Fermi gas is similar to those of the polar 
phase of superfluid liquid ^He. This phase, however, cannot be realized in experiments because 
it has higher energy than experimentally observed A and B phases (see, e.g. Ref.^^^. Note that 
several heavy-fermion compounds in a superconducting state (for a review of superconducting 
phases of heavy-fermion compounds see, e.g. Refs.,'^,'^ and^^^ also have lines of zeros of 
the order parameter on the Fermi surface and, as a consequence, the T^-dependence of the low- 
temperature specific heat (see, e.g. Ref.^^^. 

4.4 BCS pairing in a trapped single-component dipolar Fermi gas. 

Similar to the bosonic case, the trap geometry strongly influences the BCS pairing in a polarized 
dipolar Fermi gas as a result of the dipole-dipole interaction. It is natural to expect that cigar- 
shaped traps are more favorable for pairing then pancake- shaped ones because the dipole-dipole 
interaction is on average attractive in the former case and repulsive in the latter one. As a result, 
the critical temperatures in cigar-shaped traps should be higher. This question was addressed in 
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Refs.tnSlandP 

In a trapped gas, the BCS order parameter 



A(ri,r2) = yd(ri-r2) (VA(ri)vir(r2)) 



(62) 



depends on both coordinates ri and r2 and not only on their difference ri— r2 as in the spatially 
homogeneous case, because the translational symmetry is broken by the trapping potential. To find 
the critical temperature Tftrap in the trap, it is sufficient to consider the linearized gap equation: 



Mri,r2) = -Ki(ri-r2) J Jr3Jr4iS:(ri,r2;r3,r4)A(r3,r4) 



(63) 



with the kernel 



i^f^ ^ V tanh(.^v,/2r) +tanh((^v,/2r) 



Vi,V2 



where i^y = £v — M and 0v(r) are the eigenenergies (shifted by the chemical potential jl) and the 
eigenfunctions of the single -particle Schrodinger equation in the trap 



2m 



A + t/tr(r) Uv(r) = ev<|)v(r) 



where f/tr(r) is given by Eq. (13). Note that the gap equation (63) does not contain the mean 



field Hartree-Fock corrections. They lead to unimportant change of parameters of the Hamiltonian 



(54) and, therefore are not relevant for pairing. The many-body contributions 5V to the effective 



interparticle interaction which are also absent in Eq. (63) [compare with Eq. (56)], will be taken 
into account later. 

The results of the analysis of Eq. (63) in Ref.U^ are the following. When (Op ~ (O^ ^ 



(shallow nearly spherical trap), where Tc is the critical temperature of the BCS transition in a 
spatially homogeneous gas with the density n equals to the central density hq in the trap, it is 
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convenient to perform the Fourier transformation with respect to the relative coordinate r = r i — r2, 



A(R,p) = J Jrexp(-/pr/^)A(R + r/2,R-r/2), 

where R = (ri + r2)/2 is the coordinate of the center of mass. The characteristic scale of the p- 
dependence of A(R,p) is of the order of the Fermi momentum pf, while for the 7?-dependence it 
is of the order of the size of the cloud Rtf, which is much larger than the typical size = pf/mTc 
of pairing correlations (coherence length), Rtf ^ <^o- We therefor can write the order parameter 
on the local Fermi surface as A[R,p = ii;?^(R)] = A(R)(po(n), where the function A(R) obeys the 
equation (see Ref.l^for details) 



7C(3) / PF 



48;r2 



^ + 2S^)[^(^)-^^7^^^^)' '''' 



with fx = fy = 1 — 3/;r^, = The solution of this equation which is formally equivalent 

to the Schrodinger equation for a three-dimensional anisotropic harmonic oscillator, for the lowest 
eigenvalue gives the following expression 



Tr-Z 



ctrap 



In 



Tc 



' ctrap 



7C(3) neF 



24nj2 



(On 



+ 



(O7 



1 + 



(65) 



for the change of the critical temperature due to the presence of the trapping potential. According to 
this expression, the critical temperature in the trap is always smaller than that in the homogeneous 
gas. In the considered case we have (Oa/Tc <^ 1 but kCf /2An(f- > 1 (weakly interacting gas). 
Taking into account that IC, (3)/(48;r2) ^ 0.018 we see that the difference between T^trap 

and Tr is 

small if K£f /lAnd^ is not very large. 

The order parameter just below Tc is given by the corresponding eigenfunction and has the 
Gaussian form 

A(R)ocexp(- £ 4/2/L), 

a=x,y,z 
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Figure 12: The function F versus the inverse trap aspect ratio /. (Taken from Ref.lii^ 

where 



, Pf ICOa 




48;r2 V 24nd^) 



-n 1/4 



mCOa \l Tc _ 

is the characteristic size in the a-th direction. If again the quantity nep /lAnd^ is not very large, 
we have /Ace ^ ^tf^ where Rjp = pf/mcOa is the Thomas -Fermi radius of the trapped gas cloud 
in the a-th direction. This means that the pairing takes place only in the central part of the gas. 

For a large number of particles the gas is in the Thomas-Fermi regime [Eq. ([43])], and we have 
the following relation 

A^ = jU(no)V6ft)ztt)p, 

between the total number of particle A'^ in the gas and the density hq in the center, where /i(no) = 
{671^?? nof-)^ /Im is the chemical potential. For a fixed and hq, this gives 



fctrap ■ 



CO 



1 + 



lAnd^ 



(66) 



where (o = {(0^0)^)^/^ and I = -J (Op/ CO;, is the trap aspect ratio. The function 
F{1) = ^l^{3)/48n^[2^1 - + ^1 +6/;r2/-4/3] ^j^own in Fig. 
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The minimum 
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of F{1) provides the optimal trap aspect ration I* = 1.23 (cigar shaped trap, as expected) for the 
highest critical temperature in the trapped gas with the fixed total number of particles and the 
density in the center hq [the critical temperature Tc of the homogeneous gas as a function of no 
is given by Eq. ([59])]. The optimal value is a result of the competition between the anisotropic 
dipole-dipole interparticle interaction and the finite-size effects: The former favours larger / (cigar 
shaped traps) while the latter, due to the zero boundary condition on the order parameter, acts on 
pairing destructively imposing an upper limit on I. 



Numerical solution of the linearized gap equation (63 1 for the case jl ^ (U^, (Op ~ Tc shows the 
existence of a critical trap aspect ratio Ic for a given interaction strength F = 36n{0)d^ / nil and a 
number of particles A^. This critical aspect ratio corresponds to zero critical temperature, r^rap = 0, 
such that no pairing is possible in traps with I < Ic- Alternatively, for fixed values of the trap 
aspect ratio / and the number of particles A^, BCS pairing takes place if the interaction parameter 
r is large enough, F > Ff. The existence of the critical values Ic and Tc can be understood as 
follows: Due to the fact that the order parameter changes sign when the direction of the z-axis is 
reversed, single-particle states involved in forming the order parameter should have different (by 
an odd integer) quantum numbers n^. As a result, these states have different energies (the minimum 
difference is ft),) and can be paired only if the energy gain due to pairing (which is of the order 
of the critical temperature Tc) exceeds this difference. In the limiting case of an infinite pancake 
trap with the confinement only in the z-direction, this results in the appearance of a critical trap 
frequency (o^c = l.STc^^ with no pairing possible in traps with G)^ > ft)zc- 

The dependence of the critical interaction strength Fc on the trap aspect ratio / for different A'^ 
is shown in Fig. [T3| As expected, the critical interaction strength Ff decreases with increasing /. 
On the other hand, with increasing the number of particles A^, the interaction parameter becomes 
larger, F ~ N^^^, and the critical aspect ratio Ic decreases. 



The order parameter Ao(R) for a cigar shaped trap with / = 2.2 (see Fig. 14) exhibits a non- 
monotonic behavior with the distance from the trap center. This is to be compared with a monotonic 
behavior of the BCS order parameter in a two component Fermi gas with a short-range attractive 
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Figure 13: Critical lines Tc versus the aspect ratio / for different numbers of particles A'^. The BCS 
pairing takes place above the depicted curves. (Taken from Ref.Ui^ 
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Figure 14: The order parameter for the aspect ratio / = 2.2 (cigar shaped trap). Th e solid line 
shows A(z, p = 0) and the dotted line corresponds to A{z = 0,p). (Taken from Ref.^^^ 

interaction'HIL'I^ under similar conditions. (It should be noted, however, that an oscillating and 
highly non-monotonic behaviors of the order parameter in a two component Fermi gas was ob- 
tained in Ref.^in the regime of an intershell pairing <ti CO.) 

4.5 BCS pairing in a two component dipolar Fermi gas 

Adding the second component in a dipolar Fermi gas opens the possibility for a singlet interspecies 
pairing (equivalent to the ^•-wave pairing in a two component atomic Fermi gas) in addition to the 
triplet intraspecies pairing considered in Sec. |4.3[ and competition between them. Assuming 
chemical stability of a mixture of two species of polarized dipolar particles with equal masses, 
concentrations, and dipole moments (a mixture of fermionic polar molecules with two different 
hyperfine states, for example), this problem was considered in Ref .l^^ The corresponding Hamil- 
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tonian reads 

H = '£j JrVAt(r) -£v2-^ ^^(r) + l £ | JrJr>t(r)v/;(rOy„«'(r-r')V^„'(rOv^«(r), 

where a = ± denotes two different species, the intraspecies interaction Vaa(r) is the dipole-dipole 
one, Vace(r) = Vdi^), while the interspecies interaction Vaa'i^) for Gc' ^ a contains a short-range 
part Vs{r) = {ATZas/m)5{r) in addition to the dipole-dipole interaction, ^^^'(r) = V'^/(r) -(-y^(r), 
parametrized by the 5-wave scattering amplitude Us. 
The order parameter of singlet interspecies pairing 

A.(ri -rz) = y_+(ri-r2) (v/_(ri)vir+(r2)) 

in momentum space is now a sum of all partial waves with even angular momentum L and zero 
azimuthal quantum number M, 

A,(p)= £ A,z,(p)Y^o(p). (67) 

evenL 



The corresponding gap equation is similar to Eq. (56) with y(p,p') = ^^(p — p') + ATZas/m-\- 
5V(p,pO- 

hi the BCS approach, when one neglects many-body contributions including the deformation 
of the Fermi-sphere in momentum space, the critical temperature of the singlet superfluid transition 
is (see Ktf^ 



and the corresponding order parameter parameter on the Fermi surface reads [compare with Eq. ( 6 1 1] 



4- (pfP) ~ cos [cos dp a/3/|A,.|] , (69) 
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where \Xs\ is the largest positive root of the equation 



The order parameter ( 69 ) is now a symmetric function under p — )■ — p as it should be for a singlet 
pairing. 

With the deformation of the Fermi surface taken into account, the problem was solved numer- 
ically (see Ref.f^ for details). It turns out that without contact interaction (as = 0) the critical 
temperature of the triplet intraspecies pairing is always higher and, therefore, the system under- 



goes the transition into the triplet intraspecies BCS state, see Sec. 4.3 The ground state in this case 
is a mixture of two intraspecies triplet superfluids. However, switching on an attractive contact in- 
teraction (as < 0) one can increase the critical temperature of the singlet pairing and make it larger 
than for the triplet one such that the ground state corresponds to an interspecies singlet superfluid. 
The critical value of askf as a function of D = nd^/Sf and the corresponding phase diagram are 



shown in Fig. 15 



4.6 BCS pairing in a dipolar monolayer 

Let us now consider the possibility for superfluid pairing in a polarized dipolar monolayer - a 
polarized single-component dipolar gas confined to a (quasi)2D geometry by a harmonic trapping 
potential V{z) = m(0^z^/2, assuming hco^ ^ £f, where £f = pp /2m is the Fermi energy of a Fermi 
gas with the 2D density njo = Pp/4-7th^. An important parameter of the problem is the angle 
between the z-axis (normal to the 2D motion plane) and the direction of the dipole polarization, 
see Fig. [16} The effective 2D interaction between dipoles is given by Eqs. ( [49] ), ( [50[ ), and pT) . 



To see the possibility of BCS pairing, one has to look at the component in the ;?-wave channel. 
Straightforward but lengthy calculations'^^ show that it becomes negative when sin0 > 2/3, and, 
therefore, the system becomes unstable against BCS pairing. Note that the critical value 6c = 
arcsin2/3 = 0.73 (42°) is larger than the value arcsin 1 /v/3 = 0.62 (35°) of the angle above 
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Figure 15: The dependence of the critical scattering length a* on the strength o f the dipole interac- 
tion D = nd^/£f, and the resulting superfluid phase diagram. (Taken from Ref.^^^ 
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Figure 16: Fermionic monolayer: Fermionic dipoles are confined to the xy-plane. The dipoles are 
ahgned in the ;cz-plane and form an angle G with the z-axis. 
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which the dipole-dipole interaction has attractive directions in the -plane. 

Due to the anisotropy of the effective interaction (for > 0), the azimuthal quantum number 
is no longer a conserving quantity and, therefore, the superfluid order parameter contains all odd 
harmonics in the azimuthal angle cp. In momentum space, 

oo 

A(k) = £A„(fc)cos[(2n-l)(pk] 

with the components A„(fc) coupled to each other through the gap equation. With only n = I 
component taken into account, the problem was solved in Ref.,^^^ in which is was also found 
that the Fermi surface deformation due to anisotropy of the effective interaction does not play 
any significant role in the pairing problem and only slightly decreases the critical angle Oc when 
the strength of the interaction increases. This result was confirmed by the analysis of Ref.,'^in 
which all components of A(k) were taken into account. It appears that the critical angle remains 
practically the same, and higher components of A(k) with n > 1 are visible only at angles very 
close to Be- For larger values of 6, one has A„>i/Ai < 10^^ such that the pairing has indeed the 
j9-wave character, A(k) ^ Ai (k) cos (p^. 

The critical temperature for the BCS pairing readsl^^ 



2e^ 

— e^exp 

7t 



3n 1 



AkfUd (9/4) sin^e-l 



Fi{e)F2{kFlz, 0), 0>0c = arcsin2/3, (70) 



where 



and 



Fi(0) 



0.52 - 2.47 sin^ + 2.83 sin^ 
0. 18- 0.81 sin2 + 0.91 sin^e 



F,(,,e)^,/(^), f^e^^ 0-25-i.i3sin^e + i.2Ssiny 

^' ' ' ' ' 0.18 -0.81 sin^e +0.91 sin^e 



The two functions Fi{0) and F2{kFlz,0) in Eq. (70) describe deviations from the simplest BCS 



approach (with no account of the Fermi surface deformation): The exponent in Eq. (70) is simply 



the p-wave component of the effective dipole-dipole interaction on the undeformed Fermi surface. 
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Figure 17: The critical temperature as a function of the tilting angle 6 for adkp = 2.5 and kfl ^ 0.2. 
These values correspond to a gas of fermionic ^^ND3 molecules with the density n2D = lO^cm^^ 
and C0z = 2nx lOOKHz. (Taken from Ref.-^^) 

These deviations result from many-body corrections to the interparticle interaction and Fermi- 
liquid effects, as well as from the second order contribution to the two-body quasi-2D scattering 
amplitude. The latter includes virtual transitions to intermediate states which are not necessarily 
limited to the ground state of the transverse confinement. The virtual transitions to excited states 
of the transverse confinement, together with many-body contributions (those include corrections 
to the interparticle interaction, Fermi surface deformation, and Fermi-liquid effects - effective 
mass), contribute to the function Fi{6). The virtual transitions to the states with ground state 
motion in the transverse confining potential contribute to both functions Fi{0) and F2{kfl^,9) 
(see Ref.'^ for details). The latter transitions provide the second order Bom term to purely 2D 
scattering amplitude on the effective potential ( [49] ), which has the so-called anomalous scattering 
contribution ~ fc^lnfc (see more discussions below) due to long-range power decay of the potential. 



This anomalous contribution gives rise to the function F2{kflz, 6). Eq. (70) predicts a rapid growth 
of Tc with increasing of the angle 6 from the critical value Oc = arcsin2/3, to the values of the 



order of tens of nK for realistic experimental parameters; see Fig. 17 



We obtain our results for the superfluid critical temperature using the mean-field approach. 
However, as it is well-known, this approach in two dimensions is only applicable at zero tempera- 
ture, while at finite temperature the long-range order is destroyed by phase fluctuations and, there- 
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fore, the mean-field order parameter is zero. In this case, the transition into the superfluid phase 
follows the Berezinskii-Kosterlitz-Thouless (BKT) scenario,'^!!,^^ and.^^ In the weak coupling 
limit, however, as it was pointed out by Miyake,'^the difference between the critical temperature 
calculated within the mean-field approach Tc and the critical temperature of the BKT transition 
7b KT can be estimated as Tc — Tbkt Iv^ and, therefore, small as compared to Tc- As a result, 
our mean-field calculations provide a reliable answer for the critical temperature in the considered 
weak coupling regime a^kp < 1 . 

Another interesting possibility to create a topologically nontrivial superfluid state in a mono- 
layer was considered in Refs.'^^*^'and,'^^^'in which the authors consider a monolayer with fermionic 
dipoles oriented perpendicular to the plane {0 = 0) and use the RF-dressing technique discussed 
to create an effective (time- averaged) attractive potential Veff (P — > °°) ~ —din/ P'^- At 



in Sec. 



2.2.1 



short distances, the potential Veff (p) has a repulsive core that prevents low-energy particle from ap- 
proaching each other and, therefore, suppresses inelastic collisions resulting in losses. The leading 
p-wave 2D scattering amplitude for identical fermions was found^ to be 



h{k) = [%i{kp)V,ff{p)xifl+\p)2npdp ^ -l-kr, + ~{kr,)Hn{Ckr,), (71) 
Jo 3 m 2 m 

where Ji(z) is the Bessel function, ^fj^^\p) is the radial wave function of the p-wave relative 
motion, the constant C is determined by the behavior of the potential at short distances, and the 
length scale = mdl^^/fP' depends on the details of the RF-dressing (for BCS pairing has to 
be positive, for more details see in Ref.^^^. The first term in Eq. (71 1 corresponds to the p-wave 
Bom amplitude for the potential Vefflp) = —{f^/m)r^p^^ while the second term contains both 
the anomalous scattering due to the dipole-dipole tail of the interparticle interaction and the short- 
range contribution. 

For kr^ <^ 1, the first term in Eq. ( [7T| ) is negative and, therefore, leads to p-wave superfluid 
transition. The most stable low-temperature ;?-wave superfluid phase in 2D has + ipy symmetry, 
Ak = A(fc) exp(z(Pk), because this is the only jc-wave superfluid phase with a non-zero energy gap 
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on the entire Fermi sphere. The numerical solution of the gap equationj^^^^EH] shows that A{k) 
raises linearly for k<kf and approaches a constant ne^^Tc for k>kf. The critical temperature Tc 
reads 



37t 



where 

K ~ 0.16exp -— — 
V 64 

with the numerical coefficient A determined by short-distance behavior of the interparticle interac- 
tion. Note that the value of the critical temperature is very sensitive to the short-range part of the 
effective potential Veff(p) and, by modifying it, can be varied within a few orders of magnitude, 
reaching the values of the order of tens of nK for realistic experimental parameters that correspond 
to the life-time of the system of the order of seconds; see Ref. for details and discussions. 

The resulting + ipy superfluid pairing spontaneously breaks time -reversal invariance (the 
degenerate time-reversal partner is the px — ipy state). This phase belongs to the class of the 
so-called topological superconductors and can exist in one of two topologically distinct phases, 
depending on the sign of the chemical potential jU.'^^The phase with /i < is topologically trivial 
(may be continuously deformed to the vacuum state), while the phase with /i > is topologically 
non-trivial ( cannot be continuously deformed to the vacuum) and has several very interesting 
properties. One of the most interesting of them is that the vortices in this superfluid carry localized 
zero-energy states, described by a Majorana fermion. These Majorana states obey non-Abelian ex- 
change statistics j^^^l^^^l and can possibly be used for topologically protected quantum information 
processing.^i^ In the considered case of superfluid state of dipoles in a monolayer, the chemical 
potential is positive, /i > 0, and the resulting superfluid phase is topologically non-trivial. 



4.7 BCS pairing in a bilayer dipolar system 

The single-component fermionic bilayer dipolar system (see Fig. [18]) provides an example of a 
relatively simple many -body system in which an entire range of nontrivial many -body phenomena 
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Figure 18: The setup of the dipolar bilayer system: Two layers with the thickness Iq of a ID optical 
lattice formed by two counterpropagating laser waves with wavevectors and — k/, are filled with 
dipoles oriented perpendicular to the layers. The interlayer distance lisn/ki. An interlayer Cooper 
pair or molecule 
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are solely tied to the dipole-dipole interparticle interaction with its unique properties: long-range 
and anisotropy. The long-range character provides an interparticle interaction in single-component 
Fermi gases inside each layer that otherwise would remain essentially noninteracting. For the 
considered setup, this intralayer interaction is always repulsive and gives rise to the crystalline 



phase for a large density of particles (see Sec. 2.2.3 1. More important, the long-range dipole-dipole 



interaction couples particles from different layers allowing them interact in the 5-wave channel 
that is dominant at low energies, allowing formation of bound states and BCS pairing,'^.'^^ The 
Hamiltonian of the system reads 

H= £ I JrVA^(r)|-|-A+im£0,V„-M'}v^«(r) + ^/ JrJr>^(r)y/t(r/)y^(r _r')v>/3(r')V^a(r),(72) 

where a = ± is the layer index, z± = z± //2, with I being the layer separation, VAce(r) with r = 
(p,z) is the field operator for fermionic dipolar particles (p = xtx + y^y) on the corresponding 
layer a, A = Ap + / is the Laplace operator, (O^ is the confining frequency in each layer such 
that 1^ = ^Jh/mG)y, and /i' is the chemical potential. The last term describes the intra- (a = /3) 
and interlayer (a 7^ /3) dipole-dipole interparticle interactions. Assuming a strong confinement, 
hG)v ^ lJi',T, where T is the temperature, we can write ^a{^) = Wa{p)'l>Q{za) and, therefore. 



reduce the Hamiltonian (72 1 to 



= L^/^P¥^a(p)|-^Ap-Ai|v/a(p) 

+ [dpdp'xirUpMip')Vapip-p')Mp')Mp)^ (73) 

for a two-component fermionic field Xjrcdp), Ct = ±, with shifted chemical potential 11 = 11' — 
h(0z/2. The intracomponent (intralayer) interaction Vcea(p) coincides with V^^{p) in Eq. (49 1 for 
0=0, and the intercomponent (interlayer) interaction is 
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The considered system is characterized by three characteristic lengths: the dipolar length = 
md-^/h^, the interlayer separation /, and the mean interparticle separation inside each layer ~ kp^ 
with kp = V^Jtn being the Fermi wave vector for a 2D single-component fermionic gas with the 
density n. Therefore, the physics of the system is completely determined by two dimensionless 
parameters which are independent ratios of the above lengths. The first parameter g = aa/l (the 
ratio of the dipolar length and the interlayer separation) is a measure of the interlayer interaction 
strength relevant for pairing. In experiments with polar molecules, the values of the dipolar length 
ad is of the order of 10^ 10'^ nm: for a "^^K^^Rb with currently available fa 0.3D one has Ud ^ 
170 nm (with fti 600 nm for the maximum value d w 0.566D), and for ^Li^^^Cs with the tunable 
dipole moment from d = 0.35D to d = 1.3D in an external electric field ~ IkV/cm the value of 
ad varies from fa 260 nm to fa 3500 nm. For the interlayer separation / = 500 nm these values 
of ad corresponds to g < 10. The second parameter kpl measures the interlayer separation in units 
of the mean interparticle distance in each layer. This parameter can also be both smaller (dilute 
regime) and of the order or larger (dense regime) than unity for densities n = 10^ -j- 10^cm~^ 
(for, example, for / = 500 nm one has kpl = 1 for n 3 • lO^cm"^). The two parameters g and 
kpl determine the regime of interlayer scattering at typical energies of particles (~ Fermi energy 
Ep = h^kp/2m), and their product, gkpl = adkp, as usual, controls the perturbative expansion in 
the system and, therefore, many-body effects. 

The interlayer interaction has a very specific form resulting from the anisotropy of the inter- 
action: Two particles from different layers attract each other at short and repel each other at large 
distances, respectively, as a result of different mutual orientations of their relative coordinate and 
of their dipole moments. A peculiar property of V2Z)(p) is 

j dpV2D{p) = 0. 
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This means that its Fourier transform 



vanishes for small q. 



V2d{(D = I dpV2D{p)e-"'P = gqle-"^', (75) 

m 



V2D{q ^ 0) ^ gql ^ 0. 

m 



The potential well at short distances is strong enough to support at least one bound state for any 
strength of interlayer coupling. For a weak coupling between layers (g <^ 1), the bound state is 
extremely shallow and has an exponentially large size^^(see alsJ^^: 



4^2 
Eb^ — ^exp 



8 128 2521 
'? + T5^" 450""^^+ 



(76) 



for the binding energy and 



Rt = ^Jh^/mEb ~ /exp(4/^^) > / 

for the size, respectively. However, in the intermediate and strong coupling cases (g ^ 1) the size 
of the deepest bound state becomes comparable with the interlayer separation: 



Eh = {h^/ml^)2g{\-^g 

and 

Rb-mg)-"\ 

respectively. 



The specific properties of the interlayer potential (75): It decays exponentially for large mo- 
menta l^^, while it is proportional to kfor k <^ l^^ , lead to different regimes of scattering and, 
therefore, of the BCS pairing, depending on the relation between g and kfl. This can be conve- 
niently formulated in terms of the vertex function r(£',k,k'), where the arguments E, k,and k' are 
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independent of each other. This function satisfies the following integral equation^^ 

r(£,k,kO = y2D(k-k')+ / -^y2Z)(k-q) x \ .^ r(£,q,kO. iii) 

J {2ny E ~h q^/m + iO 

The 2D scattering amplitude fk{(p), where ^ is the angle between k and k', corresponds to (m/h'^)r{E,\iL, k') 
with E = fp-k^/m = h^k'^/m. The solution of this equation.l^^^H^ reads: 

For g < kpl < 1, the leading contribution to scattering is given by the first Bom term 

lit??- , 

r(£,k,k')^ ^ k-k' /; (78) 

m ' ' 

For exp(— 1/^^) <^kl < g < 1, the scattering is dominated by the second order Born contribu- 
tion 

r(£,k,kO^-— ^, (79) 

m 4 

which is momentum and energy independent and, hence, is equivalent to a pseudopotential Vb(p) = 
-i2nh^/m){gy4)5{p); 

For exp(— 1/^^) ^ kl <^ g < 1, higher order contributions become important and one has to 
sum leading contributions from the entire Born series. The result of this summation is 

r(£,k,kO-— -— 4y— , (80) 

m \n{Eb/E)+in 



where E^ is the energy of the bound state from Eq. (76 1. This expression recovers the standard 
energy dependence of the 2D low-energy scattering and has a pole at E = —Ej,, as it should be. 
The real part of the scattering amplitude, being zero at E = Eh, changes from negative to positive 
values for E > Ej, and E < Ej,, respectively. Note that within the lowest order terms, a unique 
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expression for the scattering amplitude can be written as 



,2 r 



2 



r(£,k,kO 



^|k-k'|Z 



\n{Eb/E) + i7t 



m 



As one can see, the interlayer scattering amplitude is negative in the 5-wave channel [for the 
case exp( — 1/^^) ^kl <^ g < 1 this requires E ep ^ Ej,, which is realistic in the limit ^ < 1]. 
This means that at sufficiently low temperatures, the bilayer fermionic dipolar system undergoes a 
BCS pairing transition into a superfluid state with interlayer 5-wave Cooper pairs, characterized by 
an order parameter A(p) ~ {V-i'P)¥+i ~P) ) with x(fa (p) being the field operator in the momentum 
space, which is independent of the azimuthal angle cp, A(p) = A{p). 

The analysis of the corresponding gap equation was performed in Refs.'^ and^^ (the latter 
includes many -body effects), and we present here only the results for the critical temperature Tc in 
the experimentally most interesting case kfl ~ 1 and ^ < 1 such that gkfl = a^kf < 1 (results for 
other cases scan be found in Ref.:^^ 




n 1 



X exp 



Agkplyikpl) 1 - {A/n)gkFly{kFl)a{kFl) 



T:{g,kfl), 



7t 



where 




L„(z) and I„(z) being the modified Struve and Bessel functions, respectively. 




with 
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Figure 19: The function f{x). (Taken from Ref. 



and the function f{x) is shown in Fig. 19 The dependence of the function x{g,kfl) on kfl for 



several values of g is shown in Fig. 20 We see that the critical temperature decreases very rapidly 
for kfl > 1 due to the fast decay of the scattering amplitude. The optimal value of kpl is around 
0.5 with the critical temperature reaching values of the order of 0.1/i for ^ 0.9 that corresponds 
to gkfl^OASK 1. 

This BCS state with interlayer Cooper pairs occurs in the weak (interlayer) coupling regime 
when the size of the bound state is larger than the interparticle separation (in other words, the 
Fermi energy is larger than the binding energy). With increasing interlayer coupling, the BCS state 
smoothly transforms into a EEC state of tightly bound interlayer molecules when the interparticle 
separation is larger than the size of the bound state; see, for example, Refs.'^ and.^^ Of course, 
the EEC regime and EEC-ECS crossover are possible only when the mean interparticle separation 
in each layer is larger than the distance between the layers. 

Let us now discuss possible physical realizations of the interlayer pairing. In the experiments 
with polar molecules, the values of the dipolar length aj are of the order of 10^ ^ 10^ nm: for 
a "^^K^^Rb with currently available d ~ 0.3 D one has ~ IVOnm (with aj ~ 600nm for the 
maximum value d ~ 0.566D), and for ^Li^^'^Cs with a tunable dipole moment from d = 0.35D 
to J = 1.3 D (in an external electric field ~ IkV/cm) the value of varies from aj ^ 260 nm 
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Figure 20: The function t(^,jc) for g = 0.7 (solid line), g = 0.8 (short-dashed line), and g = 0.9 



(long-dashed line). (Taken from Ref. ^) 

to ^ 3500 nm. For the interlayer separation / of the order of few hundreds nanometers, the 
corresponding values of the parameter g can be both smaller and larger than unity (g < 10). 

The values of the parameter kfl are also within this range for densities n = 10^ ^ 10^ cm^^ (for 
example, one has kpl = 1 for / = 500nm and n ^ 3 ■ 10^ cm^^). Note, however, that the optimal 



values of this parameter are around kfl ~ 0.5 (see Fig. 20) and, hence, the optimum value of the 
interlayer separation is related to the density, which, in turn, should be large enough to provide 
a substantial value for the Fermi energy. For'^^K^^Rb molecules at the density n ^ 4 - 10^ cm^^ 
in each layer one has £f ^ lOOnK and kf =. Therefore, the interlayer separation / should be 
relatively small, / < 150nm, to meet the optimal conditions. For I = 150nm one then has g ^ 1.1 
(with current d ~ 0.3D), kfl ^ 1, and Tc ~ 0. 1£f ~ lOnK. Note that strictly speaking these values 
of parameters g and kfl do not correspond to the weak coupling regime considered in this paper, 
rather to the intermediate regime of the BCS-BEC crossover. However, based on the experience 
with the BEC-BCS crossover in two-component atomic fermionic mixtures, in which the critical 
temperature continues to grow when approaching the crossover region from the BCS side, we 
could expect that the above value of the critical temperature provides a good estimate for the onset 
of the superfluidity in the intermediate coupling regime. 
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4.8 Stability of fermionic dipolar systems 

As we have already mentioned, a polarized homogeneous fermionic dipolar gas becomes unstable 
for strong dipole-dipole interaction {a^kf > 1). Unfortunately, for this strongly interacting regime 
the usage of the Hartree-Fock approximation could not be rigorously justified. It is commonly 
accepted, however, that the Hartree-Fock method is able to provide a correct qualitative picture 
even in this regime, although its quantitative results should be taken with care. 

The simplest way to get a quantitative insight into the instability region for a 3D dipolar gas 
is to use the variational ansatz (46) to calculate the compressibility. It turns outj ^^l^^l that the 



compressibility if a homogeneous gas becomes negative for aakp > 9.5 signalling the instability 
of the gas leading to a collapse. 



For a trapped gas with particles,!^ the ansatz (46) for the momentum distribution has to be 



generalized to the ansatz for the Wigner distribution function 



n{r,p) = e[pF{r)^-Upl+pj)-OC^P, 



where ppi^Y = pj^ ~ ^'^^m^il^i^^ + ^ is the position-dependent square of the Fermi 

momentum, A and /3 are variational parameters, pf = (48//)^/^A^/^/?//co, and la, = ^ffijmxi) with 
CO = (cOpCO^)^/"^. Minimization of the energy with respect to the variational parameters a, /3, and 
A gives equilibrium density and local momentum distribution for a given value of the trap aspect 



ratio / = a/ (Up /ft),. The calculation of the compressibility shows that it becomes negative for any 
value of the trap aspect ratio I provided the value of the parameter N^^^cid / la, which measures the 
strength of interparticle interaction, is sufficiently large. The dependence of the critical value of 



CQ on the trap aspect ratio / is shown in Fig. |2l| Note the found in Ref.li^ critical aspect 
ratio Ic = 0.19, below which the trapped gas is stable for any value of N^/^ad/lco, is an artefact of 
the approximation used in this paper resulted from neglecting the deformation of local momentum 
distribution, see discussion in Ref.!^ 

In a dipolar monolayer, the calculations of the compressibility in Ref.'^l^ shows that the col- 
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Figure 21: Critical value of N^/^a^kf as a function of the trap aspect ratio I. The solid line 
represents the full variational calculation, while the dashed line is obtained by forcing a = 1. 
(Taken from Ref.^) 

lapse of the system takes place for adkp ^ 1 for = ;r/2 (dipoles are oriented parallel to the plane) 
but for rapidly increasing value of a^kp when the angle between the dipole polarization and the 
plane increases. More careful analysis of the both long wavelength (collapse) and finite wavelength 
(density wave) instabilities of a dipolar monolayer can be performed by looking at the stability of 
the collective modes: A collective mode becomes unstable when its frequency tends to zero. The 
corresponding equation valid for a general wave vector k of the collective mode reads 



^co + e(p') -e(p' + ^k) 



(81) 



where r(k,p,p') = V^'^{hk) - V^'^ [p - p' + hk) for our case. [Note that Eq. (52) follows from 



(81 1 after taking the limit k ^0 and using the relation x(p) = [co — kVpe(p)]')^(p).] As a result, 
an instability occurs when the equation 



^P' rrt « ^,^ n(pO-n(p^ + ;ik) , 



{2nh) 



e(pO-e(p' + ^k) 



(82) 



has a non-trivial solution for some value of k. Note that for f (k,p,p') = V^^{hk) (when only the 
direct interaction is taken into account and the exchange one is neglected), ydip) is p-independent 
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and Eq. ([8T]) reduces to 



l-Vi^{hk)n{(O,k)=0, (83) 



where 

dp' n{p') - n{p' + hk) 



U{(0,k) 



{2%hY hco + e(pO - e(p' + hk) 



is the 2D polarization operator. Eq. ( 83 1 is used to study long wavelength {k — t- 0) plasmon os- 
cillations in electrically charged systems (see, for example, ^■^^). Although keeping only direct 
interaction in the long wave-length limit is legitimate for Coulomb systems (because of divergence 
of the Coulomb interaction, while the exchange one is finite due to non-zero momentum transfer, 
|p — p'l pp)^ this approximation gives physically incorrect results in a Fermi system with a finite 
Fourier component of the interparticle interaction for small momentum transfer (like in the consid- 
ered case of a dipolar monolayer). In this case, the direct and the exchange contributions are of the 
same order and keeping only the former results in unphysical results. Actually, for a short range in- 
terparticle interaction (with a momentum-independent Fourier component), the two contributions 
cancel each other, as it should be in a single-component Fermi gas. Similar considerations are also 
applied to the analysis of instabilities in a dipolar systems on the basis of Eq. ( [82] ): Keeping the 
exchange contribution in this equation is essential in order to obtain correct results consistent with 
fermionic statistics of particles. 



For — 7- 0, Eq. ( 82 ) is equivalent to the Pomeranchuk criterion,'^^^ formulated in the frame- 



work of Landau Fermi-liquid. Numerical solution of this equation (for k = shows that the 
instability of the system for Q > 3;r/8 corresponds to mostly isotropic with some addition of the 
quadrupole (~ cos2(p) deformation of the Fermi surface, i.e. to the collapse, that takes place 



for a^kf > 1 (see Fig. 22). For 9 < 3;r/8, the leading instability corresponds to the p-wave 
(~ sincp) deformation of the Fermi surface. However, this instability is unobservable because 
for these values of the tilting angle 9 the system undergoes a density-wave instability at smaller 



values of adkp'^^, This instability corresponds to a non-trivial solution of Eq. ( 82 ) with finite k 
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Figure 22: Phase diagram of the 2D dipolar Fermi gas at T = 0. For < a^kp < 1.5 and small 
tilting angles 9 the system is a normal Fermi liquid (NFL). Thje transition to the supefluid state 
(SF) occurs at the critical angle 9c = 0.72. At moderately strong interactions, the system either 
collapses or undergoes the transition into the density-wave phase (DW). (Taken from Ref.*^^) 

and takes place for ajkf > 1.5. For < 9 < 3n/S, the corresponding vector k is along the j-axis 
(9k = ±7r/2) and has the modulus that is twice larger than the Fermi wave- vector in the j-direction, 
k = h^^2pf{n/2)ey. For 0=0 (isotropic case), the instability vector k has no preferable direc- 
tions and the system is believed to become unstable against formation of a crystalline state with 
the triangular lattice. 

5 Dipolar multilayer systems 

Let us now briefly discuss known results on dipolar multilayer systems. The stability against 
formation of inhomogeneous (density wave) phases in fermionic dipolar multilayer system was 
discussed in Ref.^^^for a particular choice of the tilting angle 9 = arccos 1 / y/3 and in Refl^^for 
an arbitrary 9. The analysis of these papers show that inclusion of exchange interactions tends to 
stabilize the homogeneous state resulting in higher values of critical dipolar interaction strength 
as compared to the simple random-phase- approximation (RR\) approach. On the other hand, for 
multiple layers this critical dipolar interaction strength decreases with the number of layers. 

Another interesting feature of multilayer systems of dipoles is the formation of many-body 
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bound states in the form of a chain (or filament) made of one dipole in each layei^^^^^^^ (It was 
argued in Ref.li^ that bound states involving two molecules from the same layer do not exist.) 
The binding energy of such chains increases with the number of involved molecules (or layers). 
As a result, the ground state contains chains of maximum length, while at finite temperatures, the 
competition between entropy that favors shorter chains and energy preferring longer ones results in 
a non-monotonic dependence of the distribution on the length of the chains. For bosonic dipoles, 
quantum fluids of such self-assembled chains (dipolar chains fluid) and superfluidity of dipolar 
chains was considered in Ref.l^^ For fermionic dipoles, the situation is even more interesting 
because chains with even number of dipoles are bosons, while with odd number of dipoles are 
fermions. In this case, even at zero temperature, there is an interplay between the Fermi statistics 
in the form of a Pauli principle giving rise to a finite kinetic energy of a filled Fermi sphere, and the 
binding energy, see Ref.,'^ where a Bose-Fermi mixture of self-assembled noninteracting. chains 
was considered for the simple case of a three-layer system of fermionic polar molecules oriented 
perpendicular to the layers. For a more general case for both bosonic and fermionic dipoles, which 
also includes interactions between chains, see Ref.^^ 

The superfluidity in fermionic dipolar multilayer systems was addressed in Ref.'^ The inter- 
layer character of Cooper pairs in this case leads to the competition for pairing among adjacent 
layers resulting in a dimerized superfluid state as the ground state, in which the system can be 
viewed as a stack of bilayers with interlayer pairing correlations inside each bilayer and no such 
correlations between layers belonging to different bilayers. This state is characterized by a quasi- 
long-range superfluid order in every bilayer. At some finite critical temperature, this phase un- 
dergoes a phase transition into a dimerized "pseudogap" phase with only short-range superfluid 
correlations. These correlations disappear above the second critical temperature, and the system is 
in the normal phase (see details and proposals for experimental detections of the phases and phase 
transitions in Ref.^^^. 
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6 Strongly interacting dipolar gas 

Strong correlations are at the core of a number of fundamental phenomena in many-body physics, 
ranging from the formation of self-assembled ionic crystals to exotic phases such as high-Tc su- 
perconductivity and spin liquids. The regime of strong correlations between particles in a gas is 
generally obtained when the strength of the inter-particle interactions becomes comparable to, or 
larger than, the average kinetic energy. There are two main avenues to achieve this regime of strong 
correlations: (i) the first is to decrease the kinetic energy by placing particles on a lattice - which 
in the case of a dilute gas has the effect of increasing the effective mass - and the second (ii) is to 
increase the relative strength of interactions. In gases of dipolar particles strong correlations can 
be achieved either way. In the following we first discuss the phase diagram for dipoles confined to 



2D with tunable interactions (Sect. 6.1 ), and then we review several works on exotic many-body 



phases for interacting dipoles trapped in optical lattice. Sect. 6.2 



6.1 Two-dimensional dipoles: phase diagram. 

The conceptually simplest example, although remarkably rich from a physics point of view, is a 
system of cold polar molecules in a DC electric field under strong transverse confinement. The 
setup is illustrated in Fig. [2]^a). A weak DC field along the z-direction induces a dipole moment d 
in the ground state of each molecule. These molecules interact via the effective dipole-dipole in- 
teraction VjP(r) = D{r^ — 3z^) / according to their induced dipoles, with D = d^. For molecules 
confined to the ;c, j-plane perpendicular to the electric field this interaction is purely repulsive. For 
molecules displaced by z > r/\/3 the interaction becomes attractive, resulting in few -body and 
many-body instabilities. As discussed in Sect. 2.2[ these instabilities can be suppressed by a suffi- 



ciently strong 2D confinement with the potential Vtrap(Z;) along z, due to, for example, an optical 
force induced by an off-resonant light field.!^ 
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Figure 23: Tentative phase diagram for bosonic dipoles in 2D [see setup in Fig. [2| in the T — 
plane: crystalline phase for interactions > Tqm and temperatures below the classical melting tem- 
perature Tm (dashed line) The superfluid phase appears below the upper bound T < nfp'n/lm 
(dotted iine).l^^^l^ The crossover to the unstable regime for small repulsion and finite transverse 
confinement Q for polar molecules is indicated (hatched region). Hatched regions for ~ rqm and 
at the crystal-normal phase transitions correspond to the existence of possible exotic phases [see 
the Open Questions section in the text]. [Adapted from Ref.'^ 
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The 2D dynamics in this pancake configuration is described by the Hamiltonian 



which is obtained by integrating out the fast z-motion. Equation ( 84 ) is the sum of the 2D kinetic 
energy in the x,y-p\ane and the repulsive 2D dipolar interaction 

Kff (P)=^Vp', (85) 

with p-j = {xj —Xi.yj —yi) a vector in the ;c,y-plane [solid line in Fig. [3|a)]. The distinguishing 
feature of the system described by the Hamiltonian ([84]) is that tuning the induced dipole moment 
d drives the system from a weakly interacting gas (a 2D superfluid in the case of bosons or a 2D 
Fermi liquid'^^, to a crystalline phase in the limit of strong repulsive dipole-dipole interactions. 
This transition and the crystalline phase have no analog in the atomic bose gases with short range 
interactions modelled by a pseudopotential of a given scattering length. 

A crystalline phase corresponds to the limit of strong repulsion where particles undergo small 
oscillations around their equilibrium positions, which is a result of the balance between the repul- 
sive long-range dipole-dipole forces and an additional (weak) confining potential in the x, j-plane. 
The relevant parameter is 

Epot (f/a^ (fm ad 
i^kin n jma'^ n a a. 

which is the ratio of the interaction energy and the kinetic energy at the mean interparticle distance 
a. This parameter is tunable as a function of d from small to large. A crystal forms for ^ 1, 
when interactions dominate, which for a dipolar crystal is the limit of large densities. This density 
dependence is different from that in Wigner crystals with \/r- Coulomb interactions, as realized 
e.g. with laser cooled trapped ions.^^In the latter case = (e^/a) jt?' jma^ ~ a and the crystal 
forms at low densities. In addition, the charge e is a fixed quantitiy, while d can be varied as a 



88 



function of the DC field. 



Figure 23 shows a schematic phase diagram for a dipolar gas of bosonic molecules in 2D as 
a function of and temperature T . In the limit of weak interactions < 1, the ground state is 
a superfluid (SF) with a finite (quasi-)condensate. The SF is characterized by a superfluid frac- 
tion ps{T), which depends on temperature T, with p5(r = 0) = 1. A Berezinskii-Kosterlitz- 
Thouless transition^I^ElII from the superfluid to a normal fluid occurs at a finite temperature 
Tkt = TtpsfP'n/lm, as expected in two dimensions. Recent numerical results in Ref. ^"^"^ obtained 
with an exact Path-Integral Monte-Carlo technique (PIMC)'^^have shown that the superfluid frac- 
tion Ps{T), and thus Tkt, has a non-monotonic behavior as a function of the interaction strength rj, 
reaching a maximum of about psiT^j) = 0.9 at ~ 1. 

In the opposite limit of strong interactions rj^l the polar molecules are in a crystalline 
phase for temperatures T < T,n with T,„ ^ 0.0S9D/a^, see Ref .1^^ The configuration with min- 
imal energy is a triangular lattice with excitations given by acoustic phonons, with character- 
istic Debye frequency TicOq / {fP' / ma^) ~ 7.9^/rj. The intermediate strongly interacting regime 
with > 1 has been investigated using several numerical, especially quantum Monte-Carlo, tech- 
niques in Refs.fSSHHHIll! Using an exact PIMC technique, a quantum melting transition from the 
crystalline to the superfluid phases has been determined to occur at a critical interaction strength 
rQM = 18±4,.lS4] 



Confined geometries: The addition of an in-plane parabolic confinement of frequency 0)\\ as 
realized in experiments with magnetic or optical traps introduces a term £,mft)|?p?/2 in Eq. ([84]). 



The strength of interactions is now characterized by the dimensionless ratio T = with 



lh/m(0\\ the harmonic oscillator length. The physics of classical mesoscopic crystals with 
T ^ 1 and dominant thermal fluctuations has been first discussed in the context of excitonic ma- 
terials. ^^^^^^^ It was found that small dipolar clusters with N <AQ particles confined to parabolic 
potentials may not crystallize in a triangular lattice, but rather arrange in shell-like structures. 
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Figure 24: (color online) (a-d) Monte Carlo snapshots of the density of particles in all mesoscopic 
phases forN =13 dipoles, as a function of the effective mass T. (a) superfluid; (b) supersolid; (c) 
ring-like crystals; (d) classical crystal. [From Ref.li^ 
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Finite-r melting of these structures would usually proceed through two separate (non-sharp) tran- 
sitions, corresponding to the loss of rotational and radial order. 

Motivated by dipolar molecular and atomic gases, the focus has now shifted towards the low- 
temperature regime of dominant quantum fluctuations J ^^^ ^ ^^^ ^^^^^^ For bosons, it was recently 
determined^^^^^^that the quantum melting of mesoscopic crystals (t 1) into a superfluid (t < 1) 
is a sharp crossover involving two intermediate phases: these are ring-crystals, with vanishing su- 
perfluid fraction, and mesoscopic superfluids with a modulated (e.g., non-homogeneous) density. 



Snapshots of quantum Monte-Carlo simulations for N — 13 particles are shown in Fig. 24 for all 
of these quantum phases. 

Having determined the low-temperature phase-diagram both in the homogeneous situation 
and in confined geometries, the remaining question is whether these phases, and in particular 
the crystalline phase emerging at strong dipole-dipole interactions, are in fact accessible with 
polar molecules. As discussed in Sec. 2.2[ stable 2D configurations for the molecules exist in 



the parameter region where the combination of strong dipole-dipole interactions and transverse 
(optical) trapping confines the particles' motion to the large distance region with a > l±, with 
l± ~ {12d'^/mQ?y/^ the position of the saddle points in Fig. jsj^b) (white circles). For a given 
induced dipole d the ground-state of an ensemble of polar molecules is thus a crystal for mean 
interparticle distances l± ^ a < flmax^ where amax = d^m/h^rqM corresponds to the distance at 
which the crystal melts into a superfluid. For SrO (RbCs) molecules with permanent dipole mo- 
ment d = 8.9D {d = 1.25D), amin ~ 200nm(100nm), while amax can be several jim. Since for 
large enough interactions the melting temperature can be of the order of several jUK, the self- 
assembled crystalline phase should be accessible for reasonable experimental parameters using 
cold polar molecules. 

Open questions: Remarkably, there are still important open questions concerning the phase 
diagram of two-dimensional dipoles, and in particular the order of the phase transitions between 
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the solid and liquid phases at zero and finite temperature. While as often happens^ numerical re- 
sults are usually consistent with direct first order transitions, exotic intermediate phases may occur 
in either case: (i) at finite T ~ r,„ an intermediate hexatic phase characterized by a short-range 
positional and a quasi-long-range orientational (sixfold) order may exist between the solid and the 
isotropic liquid phases. Evidence for this hexatic phase, originally proposed by Kosterlitz, Thou- 
less, Halperin, Nelson, and Young,^ ^ ! ^ ^ ^ ^ I has been recently found numerically in Ref.^^Cii) At 
low-temperature theoretical results suggest the presence of an intermediate microemulsion bubble 
phase between the superfluid and the solid. First introduced in the context of 2D electron gases 
in Si MOSFETs by Spivak and Kivelson,'^^^' and never observed so far (however, see below), a 
microemulsion should prevent a first order transition in two dimensions. The observation of either 
phase with cold dipolar gases would constitute a breakthrough for condensed matter theory. 

Strong correlations can also occur in a weakly interacting dipolar gases subject to rotation,'^^^''^^^ 
which, as extensively reviewed in Ref.,!^ represents a key element to engineer effective magnetic 
fields in ultracold atomic and molecular gases. Strongly correlated phases similar to the Laughlin 
quantum Hall states have been proposed in Ref.'^for dipolar bosons and for in Ref.l^^for dipolar 
fermions, while the transition from the Laughlin liquid state to a dipolar crystal state was addressed 
in Ref.f^ Recent worlJ^ has now provided quantitative estimates for the realization of Abelian 
and non- Abelian gauge fields with polar molecules. 

6.2 Optical lattices 

Hubbard Hamiltonians are model Hamiltonians describing the low-energy physics of interacting 
fermionic and bosonic particles in a lattice.l^^They have the general tight-binding form 

" = -I.jHahcy+ I -^n.anj,^,. (87) 
i,j,(y ij,(y,(y' 
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Here Z7,.c7 (b] ^y) are the destruction (creation) operators for a particle at site / in the internal state 
o, Jfj describes coherent hopping of a particle from site / to site j (typically the nearest neighbor), 
and Ufj'^' describes the onsite (/ = j) or offsite (/ 7^ j) two-body interactions between particles, 
with Hi cy = bj^bi a- Hubbard models have a long history in condensed matter physics, where they 
have been used as tight-binding approximations of strongly correlated systems. For example, for 
a system of electrons in a crystal hopping from the orbital of a given atom to that of its nearest 
neighbor, o represents the electron spin. A (fermionic) Hubbard model comprising electrons in a 
2D lattice with interspecies onsite interactions is thought to be responsible for the high-temperature 
superconductivity observed in cuprates.^^ 



In recent years, Hubbard models have been shown to provide excellent microscopic descrip- 
tions of the low-energy physics of interacting bosonic and fermionic atoms trapped at the bottom 
of an optical lattice. '^^^'^ Since the interactions between cold atoms are short-ranged, in these 



systems Hubbard Hamiltonian typically have onsite interactions only [U^ ■' in Eq. ([87])]. This is 
readily shown for the simple case of single- species (o = o') bosonic atoms with contact interac- 
tions, such as ^^-^Cs atoms prepared in their absolute internal (hyperfine, fine, ...) ground-state, and 
trapped in the lowest band of an optical lattice. In the limit in which all energies involved in the 
system dynamics are small compared to excitation energies to the second band and neglecting the 
often-small overlap beyond nearest neighboring densities, the microscopic many-body Hamilto- 
nian reduces to one of the form of Eq. ( [87] ) with nearest- neighbor hopping energy J = Jf ^j^i and 
on-site interactions U = U^'^^''^ with 

J =-/<(r)(-|^ + VoW)w,+ i(r)d3r, 

U =4^/|w,(r)|4d3r. (88) 

Here, {w,(r)} is a complete set of single-particle basis functions, known as Wannier functions, 
which are linear combinations of exact solutions of the Schrodinger equation in the periodic op- 
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tical potential Vo{r) = Y.a=i3^o.a^^^^{kaXa) (known as Bloch functions), and are localized at 
individual sites j (here we focus on the lowest lattice band only). The optical lattice has a depth 
Vo,a proportional to the intensity of the confining laser beams, with wavevector ka- In atomic 
systems, the Hamiltonian parameters in Eq. ([88]) can then be accurately controlled independently 
using external (optical, magnetic) fields: by increasing the intensity of the lattice laser light, J de- 
creases exponentially, while U can be broadly tuned, e.g., by varying the scattering length using 
magnetic Feshbach resonances .l^l 

The resulting Bose-Hubbard Hamiltonian (BHH) has been extensively studied in condensed 
matter physics.'^'^ When the atom number is commensurate with the number of lattice sites, the 
BHH predicts a zero-temperature phase transition from a superfluid (SF) phase to a Mott insulator 
(MI) with an increasing ratio of the on site interaction U (due to repulsion of atoms) to the tunneling 
matrix element J. In the MI phase the density (occupation number per site) is pinned at an integer 
value n = 1, 2, and the excitation spectrum shows a gap of the order of U, corresponding to 
particle-hole excitations. When the density is not integer, the low-energy phase is superfluid for 



all strengths of the ratio J/U. The associated phase diagram is sketched in Fig. 25 as a function of 
the chemical potential /i and the ratio J/U, the lobes denoting MI regions of constant density. 

Spectacular experiments with ultracold atoms have lead to the first observation and charac- 
terization of this superfluid/Mott-insulator quantum phase transition for bosonic atoms by 
looking at the interference of the expanded cloud, the measurement of the gapped excitations and 
the (lack of) conductivity in the Mott phase. Further experimental work with fermions may re- 
solve the phase diagram of the fermionic Hubbard model in 2D by performing an analog quantum 
simulation of Eq. ([87]) with two-species cold fermions .liSQHi82| 



6.2.1 Dipoles on a 2D lattice monolayer 

The long-range and anisotropic character of dipole-dipole interactions add to the Bose-Hubbard 
model new possibilities to observe quantum phases of fundamental interest. The simplest case is 
that of an ensemble of single-species dipoles (a = o') which are all polarized perpendicular to the 
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Figure 25: (color online) (a) Sketch of the phase diagram of the Bose-Hubbard model at T = 
in the plane /i/i7 vs J /U . The regions MI denote Mott Insulator phases, with an integer average 
onsite density n. [Adapted from Ref.'^^ 
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Figure 26: (color online) (a) Phase diagram of Hamiltonian ([89]) at T = 0. Lobes: Mott solids 
(densities indicated); SS: supersolid phase; SF: superfluid phase. DS: devil's staircase. Panels 
(b-d): sketches of the groundstate configuration for the Mott solids in panel (a), with density 
n = 1/2, 1/3 and 1/4, respectively. [Adapted from Ref.ti^ 

lattice plane, resulting in isotropic in-plane interactions. This adds to the BH Hamiltonian terms 
of the kind Zi<j Ul^'^nfnJ, with f/^^ = V/rfj. 

Extended Hubbard models have been extensively studied in literature. It has been predicted 
that in 2D lattices the presence of finite range interactions gives rise to novel quantum phases, like 
the charge-density wave (checkerboard), which is an insulating phase with modulated density, and 
the supersolid (SS) phase, with coexistence of superfluidity and of a periodic spatial modulation 
of the density, different from the one of the lattice. This latter phase has particularly interesting 
and storied history. First proposed in the context of the search for the groundstate of helium, it has 
been the subject of extended theoretical and experimental investigations for almost 40 years. While 
recent experiments may have spotted it in bulk solid helium, its very own existence in free-space, 
that is in the absence of an underlying periodic potential, is a matter of active research and debate. 
With few exact theoretical tools available in the strongly-interacting regime, quantum Monte-Carlo 
(QMC) methods have so far established SS behavior in free-space to be based on defects and 
disorder mechanisms, such as the presence of superfluid dislocations and grain boundaries. 
Crucially, any commensurate bulk solid (including "^He) should be insulating,^^ notwithstanding 
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recent theoretical proposals.'^ 

Several theoretical studies have demonstrated SS behavior in tight-binding lattice models.^ ^^* ^ ^'^ ^ i ssf tTgH 
Model systems with nearest-neighbor (NN) or next-neighrest-neighbor (NNN) interactions have 
been generally considered. The variety of theoretical models and techniques which have been used 
has resulted in a zoo of predictions. A general conclusion is that SS behavior seems to be favored 
by finite-range interactions as well as finite on-site interactions, U^'^ > J^jj_i- Quantum Monte- 
Carlo methods have determined the SS to occur for the following models: i) hard-core bosons 
(infinitely-large U^'-^) on a triangular lattice with NN interactions only, for densities comprised 
between 1/3 < n < 2/3; ii) hard-core bosons on a square lattice with NN and NNN interactions 
for n < 0.25 and 0.25 < n < 0.5 between a "star" and a "stripe" solid at half filling; Hi) soft-core 
bosons (finite value of the ratio U^"^ / Jf-^^) on a square lattice with NN and n > 0.5. Phase sepa- 
ration, characterized by a negative compressibility, has been predicted to occur in several models, 
for example in the latter case (iii) with NN interactions and n < 0.5. However, its origin can be 
traced back to the finite-range character of the interactions (NN) considered in that specific model, 
and thus phase- separation may be expected to disappear in the case of the infinite-range interac- 
tions considered below. 

Dipolar atomic and molecular gases trapped in optical lattices can provide physical systems 
where the dynamics is microscopically described by extended-Hubbard Hamiltonians with long- 
range, anisotropic, interactions,'^ with Hamiltonian parameters are tunable with external fields. 
In particular, polar molecules in optical lattices can provide for strong offsite dipolar interactions, 
of the order of hundreds of kHz, decaying with distance as 1/r^. Due to these strong interactions, 
two molecules cannot hop onto the same site, and thus the particles may be treated as effectively 
"hard-core" (soft-core particles may be realized with dipolar magnetic atoms). For a lattice of 2D 
polarized hard-core dipoles, the microscopic extended Bose-Hubbard Hamiltonian in the presence 
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Figure 27: (color online) n vs. /i. (a): Solids and SS for a system with linear size L = 12 and 
J/V — 0.05. Some n are indicated, (b): superfluid and vacancy-SS forL = 16 andJ/V — 0.1. 

of long-range interactions is 



<l,j> Kj ij 



(89) 



The phase diagram of Eq. ([89]) has been recently computed for 2D triangular'^ and square'^ 
lattices using exact QMC methods and no cutoff in the range of interactions. The zero-temperature 



phase diagram on a square lattice is shown in Fig. 26 For large enough J/V, the low-energy 
phase is superfluid, for all jl. For finite J, three main solid Mott lobes emerge with densities 
n = 1/2, 1/3, and 1/4. These are the checkerboard (CB), stripe, and star solids, respectively, 
the corresponding groundstate configurations being sketched in panels (b-d). Analogous to the 



standard Mott insulating phases of Fig. 25, these lattice solids at fractional filling factor exist in 
some low-tunneling region of the /i vs. 7 phase diagram. 

Interestingly, it was found that a supersolid phase can be stabilized in a broad range of pa- 
rameters by doping the checkerboard and star Mott solids either with vacancies (removing par- 
ticles) or interstitials (adding extra particles), in accordance with Andreev-Lifshitz's scenario of 
defect-induced supersolidity.'^For example, a vacancy SS is present for densities 0.5 > n> 0.43, 
roughly independent of the interaction strength. In contrast to, e.g., case (iii) above, the long-range 
interactions prevent phase separation from occurring below filling n = 1/2 in this microscopic 
model. This offers interesting prospects for the first observation of this exotic phase using polar 
molecules. 
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Figure 28: (color online) Spatial density profile in 2D for A'^ ~ 1000 particles in a harmonic poten- 
tial. Phases are indicated (CB, SR, ST stand for checkerboard, star, and stripe solids respectively), 
(a-b) V /J = 15, li/J = 55, Q./J = 0.05 and T /J = 0.0377, with temperature annealing performed 
in panel (a).; (c) V/J = 5, /i/7 = 19, a/J = 0.01 and T/J = 0.1; (d) V/J = 20, ii/J = 51, 
^/J = 0.04 and T/J = 0.25. [Adapted from Ref.^i^ 

Both the solid/supersolid and supersolid/superfluid quantum phase transitions are second-order. 
By increasing T the SS melts into a featureless normal fluid via a two-step transition, the interme- 
diate phase being a normal fluid with finite density modulations, similar to a liquid crystal.^^ 



Devil's staircase and metastable many-body states: The large Mott lobes with n = 1/2, 1/3, 
and 1/4 are robust in the presence of a confining harmonic potential and moderate finite r/7 ~ 1, 
and thus are relevant to experiments. Interestingly, for small-enough hopping 7/y ^ 0.1 the low- 
energy phase is found to be incompressible (dn/djl = 0) for most values of /i. This parameter 



region is labeled as DS in Fig. 26 and it corresponds to a finite-7 version of the classical Devil's 
staircase. First discussed in the context of atomic monolayers adsorbed on solids,!^ the latter 
consists of a succession of incompressible ground states, dense in the interval < n < 1, with a 
spatial structure commensurate with the lattice for all rational fillings. 1^^^^^ This behavior, which 
has no analogue for shorter range interactions, is shown in Fig. [27ta,b), where the particle density 
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n is plotted as a function of the chemical potential jl. In the figure, a continuous increase of n as 
a function of /i signals a compressible phase, while a solid phase is characterized by a constant 



n. The main plateaux in panel (a) correspond to the Mott lobes of Fig. 26 while the other steps 
are incompressible phases, with a fixed, integer, number of particles, indicating the Devil's-like 
staircase. We will come back to this point in the discussion of one-dimensional models. 

Determining the exact groundstate geometry and periodicity of the solid for a given set of 
Hamiltonian parameters in the DS region is however a computationally daunting task, since i) for 



many rational fillings [e.g. n = 7/24 in Fig. 28 'a)] it would require to consider system sizes much 
larger than those accessible with reasonable numerical resources, and ii) the long-range interac- 
tions determine the presence of numerous low-energy metastable states, which for finite T can 
result in the presence of defects or in disordered structures. The stability of low-energy metastable 
states has been thoroughly studied in Refs.,'^^^'^^ where it is found that especially in larger lat- 
tices, two metastable configurations might differ by the occupation of just a few lattice sites. Thus, 
because of the presence of these metastable states, in an experiment it will be very difficult to 
reach the ground state or a given metastable configuration. This is directly reflected in the numer- 



ics: Figure 28 shows snapshots of the spatial density distribution in the lattice in the presence of a 
realistic harmonic confinement (shown is a single quadrant). Each circle corresponds to a different 
site, and its radius is proportional to the local density. In panels a) and b), }i is chosen such that 
particles at the trap center are in the CB phase, with very small T . The density profile shows a 
wedding-cake structure, with concentric Mott-lobes with density n=l/2 and 1/4, analogous to 
the shells with contact interactions. I^^SllQl] However, while the system parameters are the same in 
both figures, panel (a) shows regular CB and star patterns, while in panel (b) extended defects are 
present in the CB phase and the star is barely visible. This is due to the different preparation of 
the states in panels (a) and (b). In fact, in panel (a) temperature annealing of the system prior to 
taking the snapshot was performed in order to eliminate defects, while this was not done in panel 
(b). Defects in (b) reflect the existence of a large number of low-energy metastable states, which 
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are a direct consequence of the long-range nature of the interactions. 
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Figure 29: Ground state phase diagram for the Hamiltionian Eq. ([89]) on a triangular lattice in 2D, 
around n= 1/3. The phases are a superfluid 'SF', supersolid 'SS' and a commensurate solid at 
density n = 1/3. A transition region of the Spivak-Kivelson bubble type (emulsions) in indicated 
with the double line, gradually going over to a region of incommensurate, floating solids with in- 
creasing interaction strength. For large interaction strength, and starting around half filling, the 
supersolid phase is suppressed by emerging solid ordering (stripes at half filling and incommensu- 
rate, floating solids at other fillings). [Adapted from Ref. ^^^J 



Additional exotic phases: Interestingly, the possibilities offered by polarized dipoles in a 2D 
are not exhausted by the phases above. In Ref. first numerical evidence has been proposed for a 



bosonic lattice version of the "microemulsion" phase^^ introduced in Sect 6.1 above, for a system 



of polarized hard-core dipoles on a triangular lattice. This is shown in Fig. 29, where an emulsion 
phase is indicated as separating the Mott lobe at n = 1/3 from a low-density superfluid, suggesting 
a qualitative difference between the behavior of this complex strongly-interacting system above 
and below the Mott lobe. In view of future experiments with ultracold atomic gases with compar- 
atively large magnetic dipolar moments, such as Cr, Dy and Er atoms, it is an exciting prospect to 
investigate how finite onsite interactions, which usually favor supefluid and supersolid behavior, 
will modify this picture. 



Fermi gases: A conclusive determination of the phase diagram of 2D fermionic dipoles has 
been so-far hampered by the lack of exact theoretical tools for the study of the strongly interact- 
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ing regime. The experiments with dipolar gases have motivated several theoretical works, which 
use different approximate, often complementary, techniques. In analogy with Coulomb systems, 
for dipolar gases polarized perpendicular to the 2D plane (e.g., with isotropic in-plane interac- 
tions as discussed before) the existence of both crystal-like'^^^^^^ and quantum liquid crystal 
p]^^ggg|2Toi2iiJ (y^[^]^ interactions cut-off at NNN) has been recently discussed. In particular, using 
a mean-field-theory approach Ref.'^has focused on determining the complex structure of phases 
and phase-transitions of charge-density waves with different lattice unit cells. 

ReferencesEnmi] 

have discussed quantum liquid crystal phases, which may be obtained for anisotropic interactions. 
Analogous to the case of classical liquid crystals,'^ these phases are classified as being nematic 
and smectic according to their symmetry breaking associated with the deformation of the Fermi 
surface as compared to the isotropic case. In the nematic phase, the rotational symmetry is bro- 
ken so the typical Fermi surface has a cigar-like shape, i.e., it is stretched in one direction and 
shrunk in other directions. In the smectic phase the system is effectively in a reduced dimension, 
accordingly the Fermi surface is divided into disconnected pieces. Experiments with atomic and 
molecular dipoles will offer an enticing opportunity to test these predictions. 

6.2.2 Polarized dipoles on a bilayer optical lattice 

As discussed in Sect. |3] above, polarized dipoles trapped in a bilayer configuration can form a 
paired superfiuid phase (PSF), which originates from the interlayer attraction due to dipole-dipole 
interactions, when tunneling between the two layers is suppressed. The addition of an in-plane op- 
tical lattice in the two layers makes this situation even more rich from a physical point of view: in 
Ref-I^it was shown that the presence of strong enough intra-layer repulsion and inter-layer attrac- 
tion for bosonic particles allows for the realization of a novel phase, named a pair-supersolid phase 
(PSS), which is defined as a supersolid phase of composite particles. Using a model Hamiltonian 
with finite on-site intra-layer interactions U and NN intra-layer interactions and within a mean- 
field analysis, it was shown that the existence of the PSS phase relies on second order tunneling of 
(composite) dipoles in parameter regimes close to insulating Mott phases of (composite) dipoles. 
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similar to the discussion of Sect. |6.2.T] above. One example is given in Fig. I30| The existence of 



this exotic phase had been previously discussed for anisotropic t — J models, but no evidence 
was found. It remains an open question whether exact (e.g., QMC) calculations will confirm this 
prediction for realistic models with long-range interactions and, e.g., onsite hard-core constraint. 

-0.375 
|i/U 

-0.425 



-0.475 



0.08 




Figure 30: Phase diagram of the effective Hamiltonian of Ref.^^ describing bosonic dipoles 
trapped in a 2D two-layer configuration. The model Hamiltonian has an extended Bose-Hubbard 
form with f/.^J, = 0.025t/^''', [7^7^' = Here, o and o' label the two layers. The 

1 ,1 ~] 1 lilt II 

black full lines are semi-analytic solutions indicating the boundaries of the insulating lobes for the 
checkerboard (n — 1/2) and the doubly occupied checkerboard (n = 1). The shaded area is the PSS 
phase predicted within a Gutzwiller mean-field approach. PSS indicates a pair-supersolid-phase, 
and PSF a pair-superfluid-phase, see text. [Adapted from Ref.^^^] 



6.3 Advanced Hamiltonian design with polar molecules 

The quantum phases of Sect. |6.2.T] are based on interactions of the dipole-dipole type, see Eq. ([T|), 



and are thus present for any dipolar gas. For polar molecules, the techniques described in Sect. |2.2| 
for modifying the shape of interaction potentials allow for an advanced engineering of microscopic 
Hamiltonians. This offers new opportunities to explore exotic many-body phases in these systems. 
In the following we review recent progress by discussing situations where pure three-body effective 
interactions can dominate the dynamics (Sect. [63TT ), and the use of internal degrees of freedom for 



each molecule in addition to rotational ones can provide toolboxes for the simulation of condensed 



matter models of interest, including exotic lattice spin models (Sect. 6.3.2). 
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6.3.1 Three-body interactions 

In Sect. |2.2| techniques were introduced for the modification of the shape of the effective interaction 
potentials using external DC electric and AC microwave fields coupling rotational excitations of the 
molecules. In particular, for a circularly polarized AC field the attractive time-averaged interaction 
due to the rotating dipole moments of the molecules has been shown to allow a strong reduction, 
and even a cancelation, of the total dipole-dipole interaction. In a dense ensemble of molecules this 
can lead to the realization of systems where the effective three-body interaction V7^'-*(r,,rj,rfc) of 
Eq. ([5]) dominates over the two-body term V^^{ri — rj) and determines the properties of the system 
in the ground state. We note that, as always, direct particle-particle involve two particles only, while 
few-body interactions emerge as an effective low-energy interaction after the high-energy degrees 
of freedom are traced out. 

Model Hamiltonians with strong three-body and many-body interactions are strong candidates 
for exhibiting exotic ground state properties. This is exemplified by the fractional quantum Hall 
states described by Pfaffian wave functions which appear as ground states of a Hamiltonian with 
three-body interactions.'^^^*^^ These topological phases support anyonic excitations with non- 
abelian braiding statistic. The possibility of realizing a Hamiltonian where the two-body interac- 
tion can be manipulated independently of the three-body term has been studied in Ref.*^^ There, it 
is shown that a stable system of particles interacting via purely repulsive three-body potentials can 
be realized by combining the setup above with a tight optical confinement provided by an optical 
lattice. The latter ensures coUisional stability of the setup and defines a characteristic length scale 
(the lattice spacing) where an exact cancelation of the two-body term can occur. The resulting 
extended Hubbard-like Hamiltonian has the form 

^ = --^E^?^7 + E^"^«;+ E -^n.n^nk, (90) 

{ij) ¥j ¥.i^k " 

where Wijktiinjnj^ is an offsite three-body term. The latter is tunable independently of the two- 
body term UijUitij, to the extent that it can be made to dominate the dynamics and determine 
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Figure 31: Central plot: Phase diagram of the minimal model obtained within the semiclassical 
approach as a function of the chemical potential jl and the nearest-neighbor hopping amplitude 
J. Grey (white) regions are superfluid (phase separated), and light (dark) blue denotes supersolids 
(solids). Surrounding plots: schematic representation of the nature of some of the solid and su- 
persolid phases. The greyscale of the circles represents the filling (white = empty, black = full), 
while the length and the direction of the red arrows denotes the amplitude and the phase of the 
superfluid component. The blue lines highlight the unit cell of the different structures. [Adapted 
from Ref.^ 

the properties of the system in the ground state. This is in contrast to the common approach to 
derive effective many -body terms from Hubbard models involving two-body interactions, which 
are obtained in the J <^U perturbation limit, and are thus necessarily small.!^^ 

The phase diagram for bosonic particles on a 2D lattice with three-body interactions has been 
recently investigated in Ref.,^^^ where a rich variety of superfluid, solid, supersolid and phase 
separated phases have been found. In particular, several solid phases at fractional filling factor are 
found to evolve upon doping into supersolid phases with complex spatial structures. For example, 
the checkerboard supersolid at filling factor 1/2, which is unstable for hardcore bosons with nearest 
neighbor two-body interaction, is found to be again stable in a wide range of tunneling parameter, 
similar to the case of pure dipolar interactions discussed above. 

The phase diagram for bosonic particles on a ID lattice has been studied in Ref.,'^the corre- 
sponding phases being discussed in Sect. |7] below. 
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6.3.2 Lattice Spin models and quantum magnetism 

Lattice spin models are ubiquitous in condensed matter physics where they are used to describe the 
characteristic behavior of complicated interacting physical systems. Recent works have focused 
on realizing effective spin Hamiltonian for the simulation of lattice spin models of fundamental 
interest in condensed-matter, obtained by considering several internal states of each molecule. 

Reference'^ has shown that cold gases of polar molecules can be used to construct in a natural 
way a complete toolbox for any permutation symmetric two spin-1 /2 (qubit) interaction, based on 
techniques of interaction engineering which are extensions of those discussed in Sect. |Z2 The 



main ingredient of this (and related) proposal is the dipole-dipole interaction: it couples strongly 
the rotational motion of the molecules, it can be designed by means of microwave fields, and it can 
be made spin-dependent, exploiting the spin-rotation splitting of the molecular rotational levels. 

The basic building block in Ref.l^^is a system of two polar molecules strongly trapped at given 
sites of an optical lattice, where the spin- 1/2 (or qubit) is represented by a single unpaired electron 
outside a closed shell of a heteronuclear molecule in its rotational ground state, as provided 
e.g. by alkaline-earth monohalogenides. Rotational excitations are described by the Hamiltonian 

//in = 5N2-frNS, (91) 



with N the rotational angular momentum of the nuclei, and S the dimensionless electronic spin 
(assumed to be 5 = 1/2 in the following). Here B denotes the rotational constant and y is the 
spin-rotation coupling constant. The typical values of B are a few tens of GHz, and y is usually 
in the hundred MHz regime. The interaction describing the internal degrees of freedom is //jn = 
Vdd + LLi where //dd is the dipole-dipole interaction. 

The molecules are assumed to be trapped in the optical lattice with a separation Az ~ ry = 
(2J^/7)^/^, where the dipole dipole interaction is cf' /r^ = y/1. In this regime the rotation of the 
molecules is strongly coupled to the spin and the excited states are described by Hunds case (c) 
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Figure 32: (a) Potentials for a pair of molecules as a function of their separation r. The symmetries 
\Y\^ of the excited manifolds are indicated, as are the asymptotic manifolds {Ni,Ji;Nj,Jj). Here 
the quantum numbers are: Y = Mm +M5 where = M^^ +Mm^ and Ms = Mg^ +Ms^ are the total 
rotational and spin projections along the intermolecular axis; parity eigenvalues a = ±1 denoted as 
g{u) for ^erade(Mngerade); reflection symmetry i? = ± of all electronic and rotational coordinates 
through a plane containing the intermolecular axis, (b) Example anisotropic spin models that can 
be simulated with polar molecules trapped in optical lattices: Square lattice in 2D with nearest 
neighbor orientation dependent Ising interactions along x and z- Effective interactions between the 
spins Si and S2 of the molecules in their rovibrational ground states are generated with a microwave 
field E(?) inducing dipole-dipole interactions between the molecules with dipole moments Di and 
D2, respectively. [Adapted from Ref.!^^ 
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states in analogy to the dipole-dipole coupled excited electronic states of two atoms with fine- 
structure. Thus, while the two-particle ground states are essentially spin independent, effective 
spin-dependent interactions in the groundstates can be obtained by dynamically mixing dipole- 
dipole coupled excited states into the groundstates using a microwave field E(x,?) with properly- 
chosen polarization, frequency COf and Rabi-frequency Q., which is tuned near resonance with the 
= -> = 1 transition. 
The effective interactions in the lowest-energy states are obtained by diagonalizing the //bo = 



Hm — L(=i.2diE potential, as described in Sect. 2.2 In second order perturbation theory they read 



where {\gi), \gf)} are ground states for two molecules with Ni = = and {|A(r))} are excited 
eigenstates of //;„ with A^i -f A^2 = 1 and with excitation energies {E{X{r))}. The reduced interac- 
tion in the subspace of the spin degrees of freedom is then obtained by tracing over the motional 
degrees of freedom. For molecules trapped in the ground motional states of isotropic harmonic 
wells with rms width zo the wave function is separable in center of mass and relative coordinates, 
and the effective spin- spin Hamiltonian is //spin = (//eff(''))rei- 



The effective Hamiltonian Eq. (92 1 can in general be rewritten as' 



^effW = ^ £ <A„.^(r)af, (93) 
^ ce,jS=0 

where {a"}^^Q = {l^o^^o^^o^} and A is a real symmetric tensor. Equation ( [93] ) describes a 
generic permutation symmetric two qubit Hamiltonian. The components Aq,.? describe a pseudo 
magnetic field which acts locally on each spin and the components A^ ^ describe two qubit coupling. 

For a given field polarization, tuning the frequency near an excited state induces a particular 
spin pattern on the ground states. These patterns change as the frequency is tuned though multiple 
resonances at a fixed intermolecular separation. For example, the anisotropic spin model Hxrz = 
XxO^o^ + XyO^o^ -\- X^o^o^ can be simulated using three fields: one polarized along z tuned to 

108 



0+ (3/2) , one polarized along y tuned to 0^ (3/2) and one polarized along y tuned to 0+ ( 1 /2) . The 
strengths can be tuned by adjusting the Rabi frequencies and detunings of the three fields. 

Of particular interest is the possibility of realizing highly anisotropic spin models such as the 
following one 

<!n =LLJ(^1,j^Ij+i+^^'^^Ij^Uj^ (94) 

which was first introduced by Dou90t et al^^ in the context of Josephson junction arrays. This 
model (for ^ ^ ±71 /2) admits a 2- fold degenerate ground subspace that is immune to local noise 
up to £-th order and hence is a good candidate for storing a topologically protected qubit for appli- 
cations in quantum computing. Other Hamiltonians that can be realized include the famous Kitaev 
model. 



References'223E27|j^^yg recently focused on the realization of tunable Heisenberg-type models 
for unit filling of the optical lattice. Effective spin degrees of freedom are encoded in internal, 
e.g. rotational, states of the molecules, and molecular interactions are tuned by a combination of 

In Refs^^ it is found that 



2.2 



DC and AC fields, using techniques similar to those of Sect, 
doping the lattice with vacancies (that is, removing particles from a few sites) leads to a tunable 
generalization of the so-called t — J model, which is relevant in the context of high-temperature 
superconductivity .'^^ Typical Hamiltonian terms which are found in addition to Eq. ([87]) are of the 
form 



Hi 



1 

2' 



J 



+ VniUj + WniS] /(r,; 



(95) 
(96) 



where the operators 5^- = {nj^(j — Hj^yi)/!, = c'j^Cja', and S- = {S^V spin- 1/2 angular 
momentum operators on site j describing a two-level dressed rotor degree of freedom o (e.g.. 
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the lowest two = states of the molecule in the presence of a DC electric field along z), and 
satisfying = ^'^f- By intuitively thinking of the ground state as a (classical) dipole /ip = 

/ioz oriented along the DC field (i.e. /io > 0) and of the excited state as a dipole /ij = /iiz oriented 
against the DC electric field (i.e. ^Ui < 0), the terms = {jlo - V = {i1q + /ii)^/4, and W = 
(Hq — jU^)/2 derive from re- writing the direct dipole-dipole interaction, while the term J± = 2jUq^ 
comes from the transition dipole moment /ioi, describing the exchange of a microwave excitation 
between the molecules. One difference with the derivation of the effective spin models of Ref.!^ 
described above is that dipole-dipole interactions in Eq. ( |96l ) are used in first order, rather than 
second order, which can allow for stronger interactions. 



We note that the design of extended Hamiltonians with terms as in Eq. ( [96| ) is in fact a hot 
topic of research in atomic and molecular systems: for example, Refs .12221230] j^^yg discussed the 
derivation of similar terms for trapped neutral atoms and molecules at unit filling, while Refs.'^SniS] 
have focused on ionic particles. The 7^ = case is also studied in the context of molecules in 
Ref-f^Slin Ref.,l22^a similar Hamiltonian describes the dynamics of exciton-like interactions with 
impurities in gases of polar molecules. There, the authors considered first a clean system with 
one particle per site, subject to a very deep optical lattice potential in order to prevent particle 
tunneling: in this regime, dipole-dipole interactions induce collective excitonic modes, in sharp 
contrast with individual rotational excitations of short-range interacting atoms. After substituting 
a small number of initial components with a different kind of molecule still strongly confined by 
the optical potential, the system displays an exciton-impurity interactions between the original 
excitonic modes and the new particles, which may be treated as effective impurities, leading to a 
disordered background potential with both diagonal and off-diagonal contributions. Such a setup 
represents an ideal platform to investigate the effect of disorder and interaction induced localization 
in excitonic gases in both 1, 2 and 3D. 

Reference'^Sl has considered an ensemble of molecules in the rotationally ground state 
trapped on an optical lattice and have shown that collective spin excitations can be controlled using 
external electric and magnetic fields, in the context of the formation of Frenkel excitons-^ ^^^ I ^^^ ^This 



110 



system may be used for the quantum simulation of spin excitation transfer in many-body crystals 
without phonons. In Ref.^^^ a similar Hamiltonian has been studied in the context of molecular 
self- assembled crystals for quantum memory applications. 

The inclusion of more internal states for the particles, such as hyperfine levels, offers exten- 
sions to spin systems with larger spin.l^^ ^^^ l ^^^ l For example, the design of a large class of spin-1 
interactions for polar molecules has been shown in Ref.,!^ which allows e.g. for the realization 
of a generalized Haldane model in ID, using a strategy of interaction-engineering similar to those 
described above. Recent work in Ref.*^ has started exploring the non-equilibrium dynamics of 
molecules trapped in optical lattices in the presence of several internal states, including, e.g., the 
rich hyperfine structure often present in the ground state rotational manifolds. In particular, it was 
shown that the number of effective degrees of freedom participating to the dynamics may be in- 
deed dynamically changed for each molecule using external electric and magnetic fields, similar to 
the discussion above. This raises interesting prospects for the study of non-equilibrium dynamics 
in these complex systems. In section Sec. [7] below, we discuss some of the intriguing scenarios 
opened up in the many-body context by the accurate fine-tuning of interactions presented here, for 
the particular case of polar molecules trapped in ID. 

6.3.3 Hubbard models in self-assembled dipolar lattices 




Figure 33: Floating lattices of dipoles: A self-assembled crystal of polar molecules with dipole 
moment dc provides a 2D periodic honeycomb lattice Vcp (darker shading corresponds to deeper 
potentials) for extra molecules with dipole <^ dc giving rise to a lattice model with hopping / 
and long-range interactions Vij. 



Ill 



In Hubbard models with cold atoms or molecules in optical lattices there are no phonon degrees 
of freedom corresponding to an intrinsic dynamics of the lattice, since the back action on the 
optical potentials is typically negligible. Thus, the simulation of polaronic materialsl^^ where the 
presence of crystal phonons affects strongly the Hubbard-like dynamics of particles remains largely 
a challenge. 

ReferenceslSlEia 

have shown that a self-assembled floating lattice of molecules as 



discussed in Sect. 6. 1 can provide a periodic trapping potential for extra atoms or molecules, whose 
dynamics can be described in terms of a Hubbard model with phonons, the latter corresponding to 
vibrations of the dipolar lattice. 

The Hamiltonian for extra atoms or molecules in a self- assembled dipolar lattice is 

= cjcj + \ £ Vif\c]cjC, + Y,hC0qalaq + Y,Mqe''^'^°^c]cj{ag + al^) . (97) 

<ij> iJ 9 qj 

Here, the first and second terms define a Hubbard-like Hamiltonian for the extra-particles of the 



form of Eq. ( [87] ), where the operators c/ (c- ) are destruction (creation) operators of the extra- 
particles. However, the third and fourth terms describe the acoustic phonons of the crystal [aq 
destroys a phonon of quansimomentum q in the mode A] and the coupling of the extra-particles 
to the crystal phonons, respectively. Tracing over these phonon degrees of freedom in a strong 
coupling limit provides effective Hubbard models for the extra-particles dressed by the crystal 
phonons 

H = -Jl^ clcj + ^l^V,jclc]cjCi. (98) 

<'J> iJ 

The hopping of a dressed extra-particle between the minima of the periodic potential occurs at a 
rate /, which is exponentially suppressed due to the co-propagation of the lattice distortion, while 
offsite particle -particle interactions Vij are now a combination of direct particle-particle interac- 
tions and interactions mediated by the coupling to phonons, as given by polaronic dynamics.'^ 



Figure 33 shows two implementations of this idea in two- and quasi-one-dimensional configura- 



tions, [panels (a) and (b-c), respectively]. The distinguishing features of this realization of lattice 
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models are the tunability of interactions among crystal dipoles, and of the particle-phonon cou- 
plings. In addition, the lattice spacings are tunable with external control fields, ranging from a 
jim down to the hundred nm regime, i.e. potentially smaller than for optical lattices. Compared 
with optical lattices, for example, a small scale lattice yields significantly enhanced hopping am- 
plitudes, which set the relevant energy scale for the Hubbard model, and thus also the temperature 
requirements for realizing strongly correlated quantum phases. Phonon- mediated interactions can 
be quite long-ranged, decaying with distance as 1/r^ in ID. 



7 Dipolar gases in one- and quasi-one dimensional geometries 

In reduced dimensionality, the effects of quantum fluctuations are so relevant that the standard 
Fermi-liquid picture breaks down due to the emergence of several strongly correlated states of 
matter. The addition of interaction with long-range, anisotropic tail such as the dipolar one leads 
to the stabilization of very rich phase diagram in both purely ID systems and coupled ones such 
as ladders and planar arrays. In the following, we will first review the basic physics of the single 
tube configuration and then illustrate recent results on two- leg and multi-leg ladders. 

7.1 Dipolar gases in a single tube 

We now consider a system of dipolar particles confined in a one dimensional (ID) geometry by a 
sufficiently deep optical lattice with frequency 0)±, so that their dynamics is purely ID. In order 
to ensure coUisional stability, we will consider all dipole moments aligned perpendicularly to the 
wire direction by an external field in such a way that the inter-particle interaction is always repul- 
sive. Defining m the mass of the particles and C3 the strength of the dipolar interaction, a proper 
description for, e.g., a gas of polar molecules in the long-distance regime (see Sec.|2]) is encoded 
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in the following model Hamiltonian: 

// = [dxxir^ (x) \- — d^] v/(x) + — [dx dx' i//'^ (jc) i/a'^ (x') , \ir(x') v/(x) (99) 

where V'^(jc), V'^^(-'c) are annihilation/creation operators with bosonic or fermionic statistics. For a 
polar molecule gas, C3 = J^/Cq, where £q is the vacuum permeability and d the dipole moment 
induced by a DC electric field. After a proper rescaling of all quantities, one can see that the 
only relevant parameter is the ration between the dipolar length R^^ = 2nh^ / (niC^) and the linear 
density po and, as such, we expect all thermodynamic properties to depend only on P0R3. While 
statistical properties are expected not to play a major role in this setup since bosons are subject to 
an effective hard-core condition due to the infinite short-distance repulsion, some properties of the 
system are still affected by statistics; thus, we will at first consider a bosonic gas, and then comment 



on the fermionic case at the end of the section. Qualitative features of Eq. ( |99| ) are obtained using 
the bosonization technique, which allows to map the original interacting problem into a free one 
via the identities tl^SHllT] 

I ^5 q7 ^~ odd I even 

V/]/^(x)^^Po + ^^e''P«x £ exp[/m(t^W + w)] (100) 



n{x) 



Po + 



£ exp[i2m(t?(.x:) -|- TTpox)] 



(101) 



Here, the field ^ and -& represent density and phase fluctuations of the original field, and satisfy 
bosonic commutation relations \d^'d-{x) ,(p{y)\ = illb{x — y). As long as the long-range tail of 
the interaction potential decays faster than Coulomb-like interactions J248 | 249 ] as in the 



dipolar case, the Hamiltonian in Eq. ( |99| ) is mapped into the so-called Tomonaga-Luttinger liquid 
(TLL):l23a 



(102) 
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which is a purely quadratic Hamiltonian with a linear dispersion relation at small momenta, Co{k) ~ 
v^; the excitations in the system are sound waves of the density with correspondent sound velocity 
V, whereas the factor K, known as Tomonaga-Luttinger liquid parameter, is related to the compress- 
ibility of the system ^ via K = vn^, and in Galilean invariant systems satisfies K = Vf /v,Vf = 
hpQn/m. The TLL picture captures the entire low-energy physics of the original model, and rep- 
resents the substitute of the Landau-Fermi liquid scenario in ID; analogously, the parameters v,K 
play a similar role to the Landau parameters in higher dimensions. In addition, the long-distance 
decay of correlations functions is completely captured by K, so that after the bosonization mapping 
one gets: 

B{x) = ivHxMy)) (103) 



\x-y 



for the one-body density matrix, also known as superfluid (SF) correlation, and 

DW = („(.)„(>■» -PI^^A (104) 

for the correlated part of density-density, or charge-density-wave (CDW), correlation function. 
The possibility of qualitatively describing the asymptotic decay of correlations is remarkable. 
Since spontaneous symmetry breaking is not allowed in many ID models because of the Mermin- 
Wagner-Hohenberg theorem,'^^^'^^ correlation functions encode the necessary information to dis- 
tinguish between different phases in low dimensional geometries; the slowest decaying correlation 
is usually referred to as dominant order, corresponding to the more diverging susceptibility in the 
system.l^^^'^^From Eq. ( |103| ) and ( |104[ ), we thus expect a transition from dominant SF correlations 



to CDW ones at a precise value of the TLL parameter, K = 1/2, where the two decay exponents 
coincide. 

While in principle all information regarding the quantum phases in ID is encoded in the LL 
parameters, the bosonization techniques sketched above do not in allow in general to establish a 



relation between the microscopic quantities in Eq. (99) and K, as the latter has to be regarded as 
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a phenomenological parameter to be determined by comparisons with experiments or numerical 
results. Before discussing the general case, let us comment on some relevant situations. In the 
weakly interacting limit, P0R3 <^ 1, the dipolar interaction is relevant only at very short distances, 
and thus the system behaves very similarly to a Tonks-Girardeau gas,E52l253] ^^j. ^^hich Kjg = 1 ; 
as such, we expect K^^'-'^ ~ 1 in the weakly interacting regime. In the strongly interacting case, 
P0R3 ^ 1, it has been noticed in Ref .'^^^'^^ that, approximating the system as a classical crystal 
where particles are quasi-localized at equally spaced distances, one has an energy per particle eo = 
^(3) [P0R3 1 (2m)] (where C, is the Riemann zeta function), and deriving from it the compressibility 
the resulting TLL parameter in the strong coupling regime reads: 

K{pqR3 > 1) = (0.73po7?3)"^/^ (105) 



Away from these two limits, the mapping of K into the microscopic parameters of Eq. ( [991 ) 
to be performed using non-perturbative techniques. Numerical QMC simulations based on both 
reptation, diffusion and Worm algorithm have been performedj ^^^ ^ ^^^l^ -^ allowing to estimate the 
TLL parameter from the so called static structure factor 

S{k)= / Jxe"'^^(n(x)n(0)) (106) 
Jo 



and from the long-distance decay of B{x) (see Fig. 36). Moreover, QMC also allows to quantita- 
tively estimate the momentum distribution 

1 

n{k) = - / dxEixy^"" (107) 



which is directly measurable in experiments, and, according to TLL theory, displays a clear sig- 
nature of the SF/CDW transition as the typical peak of quasi-condensation at A: = disappears 
just across the phase boundary between the SF and CDW phases. Typical QMC results for n(^) 
are shown in Fig.|35]just across the phase boundary, as can be argued from the sharp decrease of 
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Figure 34: Panel (aj: TLL parameter ^ in a dipolar wire as a function of the dimensionless interac- 



tion strength poRj (here, pp = n). T he blue line is Eq. ( |108[ ), while red points and black squares are 
QMC results from Ref anS^ respectively The inset shows a similar prediction with inverse- 
square interaction (blue line) compared with the exact result (red dashed line). Panel (b): sound 
velocity to Fermi velocity ratio v/vf (vq = vp) as a function of PqR^; blue line and red p oints are 
analytical and QMC results from Ref.'^and'^^ respectively. Image taken from Ref.!^ 



n{0) when the interaction is increased. Moreover, analytical results obtained through approximate 
methods lead to the following expression;!^^ 



K 



1 + 



6C(3) 



P0R3 



-1/2 



(108) 



which compares favorably with numerical simulations, as can be seen in Fig. 36 and recovers the 
limit in Eq. ( |105| ) at strong coupling. 

Let us briefly discuss the relationship between a dipolar interacting gas and a bosonic one with 
a contact repulsive interaction, the Lieb-Liniger model.!^ Besides sharing the same low-energy 
universality class (TLL) which allows to describe thermodynamical properties of both systems on 
the same footing, there are quantitative and qualitative differences between the two cases. The 
former is related to the very different domain spanned by the TLL parameter, that is, < ^ < 1 
for dipoles, whereas ^ > 1 for contact interactions. In addition to the momentum distribution, 
such a feature can in principle be detected by investigating the so called breathing mode, which, 
as shown in Ref.,'^^ displays a remarkably different behavior. Qualitative differences are not only 
captured by to the fact that non-local interactions allow to cross the SF-CDW transition; another 
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n(k) 




Figure 35: Momentum distribution as a function of the interaction strength in a single tube as 
obtained from QMC simulation; the quasi-condensation peak at k = smoothens out with in- 
creasing repulsion according to TLL theory. Here, n = po, = i?3 fix the notations, and / is the 
system size. Image taken from Ref.l2^ 

Figure 36: Panel (a): TLL parameter ^ in a dipolar wire as a function of the dimensionless interac- 



tion strength poR^ (here, Po = n). The blue line is Eq. ( 108), while red points and black squares are 



QMC results from Ref respectively The inset shows a similar prediction with inverse- 
square interaction (blue line) compared with the exact result (red dashed line). Panel (b): sound 
velocity to Fermi velocity ratio v/v/r (vq = vf) as a function of P0R3; blue line and red p oints are 
analytical and QMC results from Ref.'^^and'^^ respectively. Image taken from Ref.l258| 



very important feature is the response to a periodic potential, which changes qualitatively when 
considering non-local repulsion, as we will see in the next section. As a final remark, we resume 
the main differences between bosonic and fermionic dipolar wires: while the TLL picture remains 
quantitatively the same between the two, single particle correlation functions are slightly different 
from Eq. ( |103| ), and, since no superfluid instability may occur, the dominant correlation in the 
fermionic case is always a CDW. 



7.1.1 Bosonic and fermionic gases in an shallow optical lattice 

The effects of dipolar interactions become even more relevant in the presence of external pertur- 
bations such as an underlying periodic potential, as realized, e.g., by an optical lattice along the 
direction of the tube.-^ In this configuration, long-range interactions compete with both kinetic 
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energy and the periodic lattice potential, which in ID has the form:'^ 

HoL = U J dx \if\x) ^m^{2Ttx/X)\if{x) (109) 

where A is the laser wavelength, Ui = U /Er the dimensionless depth of the lattice potential, and 
Er = h^/{2mX^) the recoil energy. In the limit of a very shallow lattice, Ul <^ I, the periodic 



potential may be considered as a perturbation on the TLL Hamiltonian, Eq. ( |102[ ), and its effects 
may be investigated using a sine-Gordon description.l^^^I^SlIIeTI jj^g phase diagram, and thus the 
influence of the underlying potential, sharply depends on the ratio between the mean inter-particle 
distance l/po and the lattice spacing A/2; when such ratio is an integer, 2/(Apo) = p eN, and 
for large enough interactions, the system undergoes a Berezinskii-Kosterlitz-Thouless (BKT) tran- 
gj^^Qj^liMini from a gapless (e.g., SF or CDW) to a crystalline phase, characterized by a gap in 
the excitation spectrum and broken translational symmetry with one particle pinned every p lat- 
tice sites. On the contrary, if 2/(Apo) ^ N, the system is always gapless regardless of the dipolar 
interaction strength. This remarkable behavior suggests the following argument: when there is an 
allowed configuration where particles can, at the same time, sit on the minima of the underlying 
potential and maintain a constant interparticle distance, the dipolar interactions and the lattice po- 
tential can pin the gas, as both lattice and dipolar potential energy contributions may be minimized 
by the crystalline configuration. In the opposite case, the frustration between the long-range repul- 
sion, which tries to keep interparticle distance constant, and the periodic potential does not allow 
for any stabilization of crystalline structure. 



The sine-Gordon model, combined with Eq. ( |108[ ), allows to make quantitative predictions 



on the BKT transition in the regime of a shallow lattice potential where lattice bands are formed 
(which usually corresponds to U/Er < 2) .233 Typical phase diagrams for bosonic gases are plot- 



ted in Fig. 37, the fermionic case is obtained by substituting the superfluid phase with a CDW. 
A Luttinger staircase, that is, a series of insulating states with commensurate density satisfying 
2/ (Apo) = /> G N, appears as a function of the interaction strength at a fixed Ul^^ 
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Figure 37: Panel (a): commensurate phase diagram for bosons with dipolar interactions in a shal- 
low lattice with depth Ui = OA. Physical configurations correspond to commensurate fillings 
poA/2 = l/p, with /) G N (horizontal lines are guides to the eye for p < 10). Quantum phase 
transitions from a TLL to an insulating Mott insulator (MI) occur for each \/p at the position of 
the dots on the continuous line, while red and blue dots on dashed lines indicate crossovers. Mil 
and MI2 indicate MI with different low-energy spectra (see Ref.l^^for further details). Panel (b): 
phase diagram at commensurate filling l/p = 1/3 in the Ui vs pRj plane. Continuous line: quan- 
tum phase transition between a TLL and a MI. The phase diagram for fermions is identical to the 
one for bosons, except the TLL is always a CDW. In both panels, n = po; image taken from Ref.-^ 



7.1.2 Bosonic gases in a deep optical lattice: extended Bose-Hubbard models 

As discussed in the previous section, when all energy scales involved in the system dynamics are 
much smaller that the lattice bandwidth, both bosonic and fermionic gases are properly described 
in the context of Hubbard-like models. In the bosonic, single species case, the resulting extended 
Bose Hubbard Hamiltonian (EBHH) 



db 



i=\ 



u 



ni{ni-\)+ llTii 
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contains kinetic energy (7), local {U) and dipolar (V) interaction terms. The corresponding phase 
diagram has been studied in the grand-canonical ensemble in Ref.,!^^ where a Devil's staircase 
structure has been shown to appear both in the hard-core (e.g., infinitely-large U) and soft-core 
(e.g., finite U) limits. Moreover, close to rational fillings, supersolid behavior has been suggested 
via both numerical and analytical methods, in particular close to the half-filled case, for both 
nearest-neighbor ^-^ l ^^^ l and dipolar interactions,!^^ whereas low and incommensurate filling frac- 
tion are still described by TLLs with both CDW and SF dominant orders'^ 
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Hidden order. In the special case of unit filling, that is, when the number of particles is equal 
to the number of lattice sites, an additional instability may occur. When dipolar and local inter- 
actions are much larger than the hopping rate J <^U,V, the competition between U and V gives 
rise to two different states of matter: a Mott insulator phase when U >2V, where double occu- 
pancies are suppressed, and a density wave (DW) when 2V > U, where particles rearrange in a 
periodic structure in order to minimize the non-local repulsion (see Fig. |38]). While these phases 
are separated by a first order phase transition line at strong coupling, in the intermediate interaction 
regime it has been shown that an additional state of matter, the so called Haldane insulator (HI), 
can emerge between the two .'^^ The HI is gapped and characterized in real space by a non-regular 
density pattern, where doubly occupied and empty sites are spatially separated by strings of singly 
occupied sites of uneven length (see the corresponding cartoon in Fig. [38]). The occurrence of such 
magnetic-like order (between two doubly occupied sites there is always one and only one empty 
site) together with positional disorder (the length of the strings is not constant) is encoded into a 
string order parameter. 

^stringU - k) = {5nj exp[i7t ^'^/] (111) 

j<l<k 

which approaches a constant value at large distances, |?' — j| ^ 1; here 5nj = nj—l. In addition, the 
density wave order parameter, &Dw{i — k) = (dnjSnk), which is constant in the DW phase, decays 
exponentially in the HI. Numerical simulations based on the density-matrix-renormalization-group 
(DMRG)'^^^I^ algorithm have quantitatively determined how this phase persists up to relatively 
large values of the interaction strength U ^ S close to the U ^ 2V line.'^^ In such a regime, 
where local occupations with nj > 3 are strongly suppressed, the bosonic Hamiltonian can be 
mapped into a spin-1 problem where the spin component along the z-direction satisfies 5^ = 1 — 
Hj, and the J—,V— and i7— terms play the role of spin-exchange, dipolar interaction along the 
z-axis and on-site anisotropy. The system thus resemble the so called Heisenberg model with 
single ion-anisotropy, or A — Z) model, the only quahtative difference being the presence of a full 
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Figure 38: Left panel (image taken from Ref.l^^: numerical phase diagram of dipolar bosons 
at unit filling: the Haldane insulator is stable close to the U ^ 2V line up to relatively strong 
magnitude of U. Inset: string (solid line) and density wave (dashed line) correlation functions in 
the HI phase as a function of i — j. Right panels: cartoons of the magnetic phases at filling one, in 
both bosonic and spin-1 representation. 

dipolar tail instead of nearest-neighbor interaction, establishing a strong connection between the 
HI and the Haldane phase extensively studied in the context of magnetic sy stems .'^^^'^^ An exact 
correspondence between bosonic and spin problem can be obtained in case of strong three-body 
losses,!^ where the on-site Hilbert space is truncated up to doubly occupied states and the HI 
phase appears at smaller values of the dipolar interaction.!^^ The existence of such HI has been 
established even at higher integer fillings by mapping the problem to integer spin chains with spin 
5 > 1.12211 

Three-body interactions. Additional interesting states of matter can be found in ID version of 
Hamiltonians already treated in 2D in the previous section. As an example,it has been shown hot 
to engineer three-body interactions with polar molecules, leading to the effective Hamiltonian in 
Eq. ([90]). In ID, the three-body term is responsible for the stabilization of an insulating phase at 
filling 2/3, which is driven by the fact that, above a certain ratio ofW/J, particles will minimize 
their potential energy by sitting in a periodic patter at the expenses of the kinetic one. Such phase 
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transition belongs to the same universality class of the pinning one described in the shallow lattice 
case, and the critical value of the TLL parameter = 2/9 at the tip of the lobe is obtained from 
general considerations.'^ The addition of a weak nearest neighbor repulsion V /J ^ I can also 
stabilize a insulating state at filling 1/2 when W^/7 ^ 1 ; remarkably, the further addition of a next- 
nearest neighbor repulsion induces a competition between density wave and bond order wave, as 
the BOW phase broadens thanks to the three -body repulsion with respect to the W = Q cdiS&W^ 

7.2 Two-species mixtures 

The high degree of control over internal states of ultracold dipolar gasesl^^ has opened the way 
towards the theoretical investigation of multispecies Hubbard-like models with additional dipolar 
interactions. While such specific models were not considered in standard condensed matter lit- 
erature, Hamiltonians with nearest-neighbor interactions (which is, in many respects, similar to 
a dipolar one close to half-filling) such as the extended Hubbard mode^^ have been shown to 
present a richer plethora of phases with respect to system with contact interactions only. The theo- 
retical advantage of multispecies ID systems with long-range interactions stems from the fact that 
one can understand their basic physical properties employing the same low-energy formalism used 
for models with contact interactions only.l^SlSlS] jj^g following, we will present an overview 
of theoretical results on two-species models with dipolar interactions, with both Fermionic and 
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Figure 39: Phase diagram of hard-core bosons with three-body interactions in the grand canonical 
ensemble. In the insulating phase at J/W < 0.3, CDW and bond order coexist. Image taken from 
Ref.^ 
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Fermi-Bose statistics. 



7.2.1 Fermi-Fermi mixtures 

Contrary to their two-dimensional counterparts, fermionic dipolar mixtures in ID may be the- 
oretically investigated with the same accuracy as bosonic ones. Defining as cJ^,c/,<7 the cre- 
ation/annihilation operators at the site i of the species o (we take here a =tii), Fermi-Fermi 
mixtures in a deep optical lattice are usually described by the following Hamiltonian: 

Hff = -tY^{cl„Ci+i,a + h.c.)+UY,ni,^ni^i + -Y,l^(^'^i,(^ + \ L i -^^n ^'.cr^j.g (112) 

where the first two terms are hopping and on-site interaction, the third one is a species-dependent 
chemical potential, and the last one represent the dipolar contribution. When dealing with strong 
onsite two-body losses, the effective Hamiltonian has to be projected onto a constrained Hilbert 
space without double occupancies: 

H^^ff) ^ ^Hff^ (113) 



where ^ = Yiii^ ~ "i,t'^;,|)- Such a picture, which reduces the on-site Hilbert space to just three 
states, is connected with the spin-1 chains (once 5? = cJ|C,,| — cj^c,- 1 is defined) already discussed 
in the context of the Haldane insulator phase for single species dipolar bosons; however, differ- 
ently with respect to spin models, the key conserved quantity is as can be inferred by 
performing a rigorous spin mapping.^^Ii The phase diagram of Eq. ( |113[ ) has been investigated at 



the mean field level and via numerical simulations based on the infinite time-evolving block dec- 
imation (iTEBD)^^^ algorithm in Ref.!^ considering equal chemical potentials and intraspecies 
interactions, /i| = ll^.V^^ = V^^ = V, and varying the interspecies one = —Vcos{x)- This 
setup can be realized by, e.g., applying a pair of ac microwave fields close to the resonance of 
both internal states.!^^ In the repulsive case, V > 0, the ground state is expected to show a dilute 
antiferromagnetic order of the form .. 1 000 4- 00 00 1 where, between molecules of the same 
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species, there is always one and only one molecule of the other species, in addition to two strings 
of empty sites. Such configurations display long-range order in the string correlation function 



-(5^exp[/;r £ S^]S}) 
k<i<j 



(114) 



and can be distinguished into two different cases: when V co^{x) > ? — /i/2, all sites are occupied, 
and the ground state is an antiferromagnetic insulator (AFM); otherwise, the string order is dilute, 
the ground state is compressible and the so called Haldane liquid (HL) phase is stabilized. The 
mean field phase diagram presenting such competition is shown in Fig.|40} and has been supported 
by numerical iTEBD results, which show how, while both the AFM and the HL have long-range 
order in ^^tring^ other quantities such as standard spin-spin correlation functions have remarkably 
different behaviors ,1^^ and in the HL phase may be quahtatively described by mapping the low- 
energy physics to a spinless TLL. 

Using a combination of both DC and microwave external fields, as discussed in Sec. [8} it is also 
possible to engineer ? — 7-like Hamiltonian with dipolar spin-exchange interactions, as described 
by the t — J — V — W Hamiltonian: 



Htjvw = -tY,{claCi+i,<y + h.c.) + Y,\Ri - Rj 



-3 
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Figure 40: Mean field phase diagram for Eq. ( |113[ ) with /i| = V-^ = = V 
Image taken from.l^^ 
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Figure 41: DMRG phase diagrams of the t — J (upper panel) and tJVW (for V = W = 1^ = , 
lower panel) models: blue dashed lines denotes Kp = \, and green thick line separate spin gapped 
(SS+SG) and gapless (SDW or SS/TS) regions. In the lower panel, the blue shaded region displays 
a charge Luttinger parameter Kp = 1 ± 0. 15 within numerical accuracy, and, as such, may contain 
a gapless superfluid phase TS/SS. Image taken from.^^ 

where = (n, ^ — n,- 1) and 5^ = cLc/ 1. A complete investigation of the phase diagram of 



Eq. ( |1 15[ ) is still lacking; nevertheless, the special case with V = W = = 0, one of the sim- 
plest experimental realization, has been numerically investigated in Ref.'^^ by means of DMRG 
calculations; in this case, the main difference with respect to the standard t — J model resides in 
the long-range nature of the spin-exchange interactions, which are limited nearest-neighbor in the 
t — J itself. At small values of J, the ground state is a spin-density wave (SDW) at all densities, 
resembling the ground state physics of the strongly repulsive Hubbard model. Then, increasing 
J, one enters first a region of singlet superfluidity (SS) with a finite spin gap (SG), and finally phase 
separation. In the SS region, inter-particle attraction is encoded into the TLL parameter of the den- 
sity (or charg^% sector Kp > I. The main difference with the phase diagram of the t — J model. 



which is reported in the upper panel of Fig. 41 is the presence of a larger gapped region; besides. 



the absence of a superfluid region with no spin gap (denoted as TS/SS) cannot be completely ruled 
out, as the broad blue region may indeed present such instability. 

The realization of the tJVW model with polar molecules presents two main advantages for 
the purpose of exploring the physics of the f — 7 model with ultracold gases: first, the accessible 
regions are not limited to small exchange couplings J <^t, and second the spin gap in the superfluid 
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Figure 42: Density distribution ng (triangles), np (filled circles) and UB + rif (squares) for various 
phases as realized with dipolar Bose-Fermi mixtures: (a) BF liquid, (b) BF Mott Insulator, (c) and 
(f) BF solids with different bosonic fillings, (d) density- wave BF Mott and (e) density-wave BF 
Uqu id. R esults have been obtained by means of numerical iTEBD simulations. Image taken from 
Ref.l^ 



region is twice larger, making such phase more stable against thermal fluctuations.'^ 



7.2.2 Bose-Fermi mixtures 

Multispecies gases with different statistics, that is Bose-Fermi (BF) mixtures, have been the subject 
of numerous theoretical and experimental efforts in the field of cold atoms, and are currently being 
investigated even in presence of dipolar interactions. In ID, the effective Hamiltonian of such 
systems when confined in a deep optical lattice is: 



Hbf 



,(J^B/F 



+ 



i,B 



a,5i<j' J' 



(116) 



where aj^ are creation operators obeying bosonic/fermionic statistics with o = B/F respectively, 
the first line contains both hopping and chemical potential terms and the last two lines onsite 
and off-site interactions. The insulating instabilities of this Hamiltonian occurring close to half 
filling have been investigated in Ref.^^^ by means of numerical simulations based on the iTEBD 
algorithm by truncating the dipolar interactions up to nearest-neighbors and considering several 
combination in the broad parameter space of {Ja,U(j§,Vfy§}. In addition to gapless liquid and 
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density- wave phases (a BF liquid, Fig.|42](a), and a BP density wave, Fig.|42](e)), BF Mott phases 



with and without density oscillations are present (a BF Mott insulator. Fig. 42 (b), and a BF Mott 
density wave. Fig. |42](d)). Finally, in the special case of half filling for both species, a BF solid, an 



incompressible phase with an alternating density pattern (Fig. 42 (c) and (f)), is also stable in the 
strongly interacting regime 1/(^5 ^ Ja- Such phase displays true-long range order and, differently 
from the Neel antiferromagnetic case realized with purely contact interactions, is stable even when 
considering equal tunneling rates Jb = Jp', moreover, it displays very different melting processes 
as the interaction are increased depending on the relative value of Vbb/Vff^^ 



7.3 Quasi- ID physics: coupling between tubes 

The long-range nature of dipolar interactions is well suited to create hybrid systems with purely 
ID dynamics but 2D interactions, such as arrays of coupled tubes divided by very deep optical 
lattices preventing inter-tube tunneling together with a proper tuning of the external DC electric 
fields in order to manipulate the inter-tube interaction, as seen in Secj2} Such systems establish 
a deep connection between the physics of several condensed matter systems such as, e.g., spin 
compounds'^ and ultracold dipolar gases which cannot be explored with contact interactions only. 
In the following, we first present recent results on the two-tube case, or two-leg ladder, and then 
consider the 2D limit where a large number of tubes, a planar array, is taken into account. 

7.3.1 Two-leg ladders 

Dipolar ladders allow to investigate rich physical phenomena by matching the advantages of re- 
duced loss rates typical of ID confinement with the possibility of considering multispecies physics, 
where the physical species index is represented by the wire one |^ In this respect, the special case 
of two-coupled tubes has been widely investigated. The effective Hamiltonian of such systems, in 

"The case of multispecies gases confined in quasi- ID geometries, which allows to treat problems with even more 
degrees of freedom, has attracted little attention so far. 
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case an optical lattice in the tube direction is applied, is very similar to a two-species gas, 



a,i a,i 

+ mi'^a,l5iiJ)m,anj,ii (118) 
a,l5 ij 

with the remarkable difference that here the anisotropic nature of the dipole-dipole interaction 
plays a prominent role. Here, a = 1 , 2 is the wire-index, and the specific shape of (r) depends 
on both the inter-wire distance g and the angle between the dipole moments and the plane the 
ladder lies on. Moreover, for bosonic particles, an additional term Li.a f^'^i.a(";.a ~ 1) /2 denotes 
the on-site interaction. 

Bosonic ladders. As a first step, it is worth considering if and how phases typical of ID setups 
such as Mott insulators, superfluids and Haldane insulators are changed in ladder geometries. In 
case of interwire attraction, both MI and SF phases appear in the limit of small {J/U <ti l)and large 
(J/U ^ 1) intrawire tunneling respectively, and a pair superfluid phase (PSF), that is, a superfluid 
ground state of composite particles made of one boson on each tube, is stabilized at intermediate 
values of U /J.^^^The Mott lobes close to the PSF phase have remarkably different excitations, that 
is, particle -hole are substituted by creation/annihilation of composite particles, further increasing 



the re-entrant shape of the Mott lobes, as can be seen in Fig. 43 The evolution of the HI phase 
in ladder setups in indeed more tantalizing, as in the case of spin-1 chains it has been shown 
that the string order of the Haldane phase is unstable towards weak antiferromagnetic interchain 
exchange perturbations In the case of bosonics coupled tubes, combining TLL theory with 
renormalization group arguments, it has been shownl^^ that, while a small intertube repulsion 
^12(0)/-^ ^ 1 changes only quantitatively the shape of the phase diagram, interchain tunneling 
terms of the type: 

H^=J±Y.(cliC,,2+h.c.) (119) 
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Figure 43: Numerical phase diagram of bosonic two-leg ladders with inter- wire nearest-neighbor 
attraction 'fn{0)/U = —0.75: white, gray and black regions represent SF, PSF and MI phases 
respectively. Image taken from Ref.l^SQI 



open an additional gapless, superfluid phase between the two: a schematic view of such changes 
is presented in Fig. 44 and has been confirmed numerically via DMRG simulations.'^^^'^ 

Fermionic ladders. Fermionic ladders present more instabilities with respect to the single tube, 
single species case. In absence of a lattice, the weak-coupling phase diagram has been investigated 
using the bosonization 

formalismMlT] 

in various setups.'^^^'^ In case of short range interwire 
repulsion, the interplay between the longitudinal confinement length i?, the interwire distance d and 
the angle 9 (see Fig.|2]) is responsible for the stabilization of three main phases. In case of intrawire 
attraction, 9 < 0.3;r, pairing on the single wire is favored, giving rise to a triplet superfluid phase 
(TS); in the opposite case of intrawire repulsion, 9 >03it, a spin density wave with alternating 
order (SDW^) appears, while in the intermediate regime 9 ~ 0.3;r the dominant order is an in- 
plane spin-density wave (SDWjj,).'^In the case of a deep lattice, Eq. ( |118| ) becomes an extension 
of the so called anisotropic extended Hubbard model (AEHM)I^^ with both large anisotropy and 
long-range, dipolar terms; while all the previous phases are expected to appear in the low-density 
limit, the possibility of tuning independently inter- and intra-wire repulsion, making the former 
stronger than the latter, leads to a spontaneous breaking of the Z2 symmetry associated with the 
ladder geometry, stabilizing a fully polarized ground state where all particles stay on the same 
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Figure 44: Qualitative changes in the phase diagram of two-coupled bosonic chains at unit filling 
and fixed U /J close to the HI phase. Upper panel: at finite interchain tunneling and no interchain 
interaction, a superfiuid phase appears between the DW and the HI. Lower panel: at finite inter- 
chain interaction and no interchain tunneling, the phase diagram change s on ly quantitatively. Here, 
V^ = Yi2{f)),t = J and tx_ — J± fix the notations with respect to Ref. ,'^from which this image 
has been taken. 



tube to minimize the inter-tube repulsion. Moreover, the addition of inter-wire tunneling terms 
(Eq. (|119[)) may stabihze additional phases.!^ 



When inter-wire interactions are turned attractive, a different physical picture arises depending 
on the population ratio between the wires, n\/n2. When n\ = n2, pairing between the tubes is 
always favored, and the system behaves like a quantum liquid of composite bosonic particles made 
of one boson per wire; such a liquid can be both a superfiuid (that is, a PSF), or a CDW depending 
on both the interwire distance and dipole strength.!^ If a shallow optical lattice commensurate 
with the particle density is introduced along the tube direction, that is its lattice wavelength A 
satisfies 2/niX G N, the system can undergo a BKT transition to a composite crystal as a function 



of the lattice depth, in analogy with what happens in the single tube bosonic case (see Sec. 7.1.1 1: 



in a grand canonical ensemble, this picture evolves establishing a Luttinger staircase of dimer 
crystals.*^ In the deep lattice case, these predictions have been checked by means of DMRG 
simulations. 

Imbalanced ladders. A completely different phenomenon occurs in presence of a density im- 
balance between the tubes. For interwire distances g larger than the inverse particle density, the 
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Figure 45: Panels (a), (b): phase diagram of fermionic two-leg ladders with inter-wire nearest- 
neighbor attraction repulsion as a function of the interwire distance-to- wire width ratio d/R and 
the angle 6 (see Fig. |2]); kfR = 0.1,0.2 in the left and right panel, respectively, where kf is the 
Fermi momentum of the single tube. Panel (c), (d), (e): cartoons of the SDW^, TS and SDW^j 
phases discussed in the text; horizontal ellipses denote pairing between particles in the same tube, 
while intertube correlation is indicated by vertical ellipses. Image taken from Ref.l^MI 



long-range nature of the dipolar interaction creates an effective fixed range interwire attraction, 
thus favoring pairing not only between two, but also between many particles. In the special case 
of commensurate densities between the wires, ni/n2 = p/q E Q, a gas of composite particles, or 
multimers, composed by p{q) particles in the first (second) wire respectively can be stabilized for 
sufficiently low densities.'^ Such a multimer liquid picture has a description in terms of low- 
energy field theory ,1^^ and allows to study generalized pairing mechanism beyond two-particle 
ones. In the special case of p = 1,^ = 2, the stability of trimer liquids has been confirmed nu- 
merically for sufficiently large tunneling and dipole imbalance between the tubes. When also a 
trapping potential and a non-ideal population ratio is considered, such trimer liquids (which are 
not stable in homogeneous setups with unmatched densities) appear in the center of the system, 
and are surrounded by both dimer liquids or singles species ones depending on the population ra- 
tion, forming a pancake-like structure of composite liquids (see Fig. [46]) as shown in recent DMRG 
simulations.!^^ 

The possibility of further tuning the numerous parameters of the ladder Hamiltonian Eq. ( |1 18| ) 
may indeed allow to identify additional exotic state of matter, as the complete phase diagram 
of these system is still largely unexplored. Moreover, before turning our attention to systems 
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Figure 46: Panel (a): typical configuration of two-leg ladder with attractive short-distance interac- 
tion between the wires. Panel (b-d): cartoons of a dimer crystal {(b)), dimer liquid (or PSF, (c)) 
and a trimer liquid {(d)). Panel (e-f): typical density distribution of strongly imbalanced wires with 
2ni < ^2 (see Ref.^^for technical details). In both cases, the density match n2 = 2n\ in the central 
part of the system indicate the presence of a trimer liquid. Image taken from Ref.l2S5J 

composed by a large number of tubes, it is worth noticing that ladder systems with more than two 
tubes and mixtures in quasi- ID geometries are in principle good candidates to study new model 
Hamiltonians in connection to many-body problems with many degrees of freedom. 

7.3.2 Planar array of tubes 

Planar arrays of tubes of dipolar gases constitute a hybrid setup where one can investigate 2D states 
of matter with physical properties typical of ID systems such as, e.g., quasi long range order of 
correlation functions. Even in this case, the effective Hamiltonian: 

^^j ^ ct,P,j,l 

has many controllable parameters: here, i/a^, Xffa creation/annihilation operator on the tube a 
(sums over Greek indices denote sums over the tube index, while Latin indices denote particles 
along each tube), pa = WaWa, the first term represents the kinetic energy over each tube and the 
second one the dipolar interactions, which depend on both strength and orientation of the dipole 
moments. In the zero-density limit, where only one bosonic particle per wire is considered and 
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Figure 47: Phase diagram of a planar array of dipolar bosons as a function of the single wire 
Luttinger parameter K and the ration between the dipolar interaction and the scattering length 
7= AVqK/ (a^U) (see text). Here, the angle (a according to notations in Ref.^^^) is fixed at k/2. 
Image taken from Ref.l2§2l 



inter-wire interaction are attractive, a quantum rough chain forms with the off-diagonal long range 
order in the M-body correlation function, where M is the number of wires. 1^ 

Li the finite density regime, the many body physics displays very different effects. In the 
bosonic case with intrawire dipolar repulsion, the competition between local interactions and long- 
range ones is encoded into a dimensionless parameter 7 = 4VQK/{ua^), where u and K are the 
sound velocity and the TLL parameter of a single tube, Vq is the strength of the dipolar interactions 
and a the intertube spacing. By varying both K and 7, one can span a broad parameter region with 
respect to the dipolar interaction, the density of particles and the 5-wave scattering length awP^ 
In addition to a series of density wave states driven by the strong dipolar repulsion, which may 
lead to both a stripe and checkerboard order depending on the relative strength between inter- and 



intra-tube repulsion, a gapless phase (denoted as (i) in Fig. 47 1 with dominant superfluid correla- 
tions along the tubes is present for large values of K (that is, for weak intratube interactions) .12^ 
Such a phase, which despite being effectively two-dimensional still preserves an algebraic decay 
of correlation functions typical of the quasi-long range order of ID systems, is known as slid- 
ing Luttinger liquid (SLL), and has been investigated even in several fermionic models related to 
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Figure 48: Left panel: schematic phase diagram of a planar array of fermionic dipolar tubes as a 
function of the angles 6 and 0; all phases are discussed in the text, and their cartoons are presented 
in the right panel, where solid ellipses denote (intra- or inter-tube) pairing and shaded ones indicate 
inter- tube coherence. Image taken from Ref.!^^ 



stripe physics and high-Tf superconductivity Moreover, the SLL phase may also have domi- 
nant density-density correlation along the tube in a tiny region at intermediate 7 (denoted as (ii)). 
In the strongly interacting regime, a collapsed regime take place, whose precise shape depends on 
the sign of the intertube interactions. The complete phase diagram for a fixed value of the angle 



is shown in Fig. 47 



The planar array of fermionic dipoles displays a different phase diagram, albeit some features 
such as the density wave patterns emerge independently on the statistics. The phase diagram of 
this system as a function of the angles 6 and cp (see Fig.|2]) have been derived in the framework of a 
generalized TLL theoryl^^by comparing the long-distance decay of several correlation functions, 
as reported in Fig. |48} When both inter- and intra-tube interactions are strongly repulsive, the 
ground state exhibits CDW (checkerboard) order, similarly to the bosonic case. Close to the magic 
angle 6 = 6c = 0.3n, where the intrawire interaction turns from repulsive to attractive, an intertube 
superfluid (5-SF) phase together with a density wave with broken particle number conservation 
along the tubes (or gauge-phase density wave, GPDW) are present. Finally, for attractive intratube 
interactions, a ;?-wave superfluid (;?-SF) is the precursor of an unstable phase toward intratube 
collapse. 
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8 Conclusions and outlook 



The many-body systems discussed in this review are examples of the variety of physical proper- 
ties which originate from the anisotropy and long-range character of dipole-dipole interactions, in 
combination with their tunability with external fields. In many physical situations quantum dipo- 
lar gases behave qualitatively different when compared to atomic gases with short-range van der 
Waals interactions, and provide us with a large number of novel quantum systems with unique 
physical properties. Understanding the many -body behavior of these systems is a very interesting 
and challenging problem with several potentially significant consequences for fundamental science 
and practical applications. 

In the present review we have focused on dipolar quantum gases represented by polar molecules 
in the rovibrational ground state. An alternative realization of a dipolar quantum gas, although in a 
completely different regime, is a gas of laser excited Rydberg atoms'^ or molecules .'^^^'^ Highly 
excited Rydberg atoms and molecules interact via remarkably strong electric dipole moments or 
Van der Waals interactions. This leads to the phenomenon of a dipole blockade and formation 
of superatoms where within a given blockade radius only a single Rydberg atom can be excited, 
and, for example, crystals of these superatoms can be formed.l^Sll in view of the finite life time 
of Rydberg states these many body phases will only exist for a comparatively short time in the 
so-called /rozen gas regime where there is no atomic motion. However, as discussed in Ref.,UiO| 
the large Rydberg dipoles can also be admixed weakly to the ground state by off-resonant laser 
excitation which provides a situation loosely reminiscent of the polar molecule case, although 
decoherence due to spontaneous emission remains always an issue. 

Finally, we remark that the tools for manipulating interactions in dipolar systems, as described 
in the present review, also provide promising ingredients for controlled entanglement between 
polar molecules, and thus possible new scenarios for quantum computing.'225H222l Loading exactly 
one atom or molecules per lattice site via a Mott insulator transition in an optical lattice provides 
us with an array of qubits. In the case of polar molecules, the qubits can he represented by long- 
lived rotational or spin degrees of freedom. Single site addressing, as developed in Refs.!^ for 
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atoms, allows to both manipulate as well as read out the single qubit. Entanglement of qubits can 
be achieved via the strong and long range dipolar interactions between molecules .I^^SHIstIqj. 
atomic case via dipolar Rydberg interactions. I^^H 
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